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ABSTPACT 

A review of the ZJ-region ionization measurements and its solar 
zenith angle variation reveals that a unified model of the D region, 
incorporating both the netural chemistry and the ion chemistry, is 
required for a proper understanding of this region of the ionosphere. 

Model calculations are carried out with a view to interpreting the 
solar zenith angle variation of ZJ-region ionization as measured on 
July 24, 1968 at Wallops Island. All input data are taken cor- 
responding to this day. 

The model developed for the neutral chemistry includes the trans- 
port terms relating to molecular and eddy diffusion. It describes the 
diurnal behavior of the minor neutral constituents formed in an oxygen- 
hydrogen-nitrogen atmosphere, in the height interval between 30 and 
120 km. Computations carried out for two cases of the eddy diffusion 
coefficient models indicate that the constituents which are important 
for the Z?-region positive-ion chemistry do not show a significant 
variation with zenith angle for values up to 75° over the Z)-region 
heights . 

In the ion chemistry model, ion-pair production rates are calculated 
for solar X-rays between 1 A and 100 A, EUV radiations from 100 A up to 
the Lyman-a line, precipitating electrons and galactic cosmic rays. Two 
cases of 0 ^ production rates due to different fluxes of precipitating 
electrons and X-rays below 3 A have been considered. The model describes 
the solar zenith angle variation of the posirive-ion composition, negative- 
ion composition and the electron densities up to 75° zenith angle, in the 
height interval between 60 and 100 km. 
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A comparison o£ the computed electron- density profiles with the 
rocket measured profiles reveals that good overall agreement is obtained 
when the neutral chemistry model associated with the low eddy diffusion 
coefficient values is used. The computed solar zenith angle variation* 
however, is only a fraction of the measured variation. This discrepancy 
is probably due to short-term variations in the mesopause temperatures 
and X-ray fluxes below 3 A, 
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1. INTRODUCTION 


I . I The Ionospheric D Region 

Solar and other extra-terrestrial radiations incident on the upper 
atmosphere of a planet ionize its neutral molecules producing pairs of free 
electrons and positive ions. The region where such free electrons occur in 
significant amounts is known as the ionosphere. The single parameter that 
characterizes the ionosphere is the electron density, which is a function 
of the height, latitude and time of the day, among other factors. 

The propagation of a high frequency radio wave through the ionosphere 
is influenced by the, free electrons present in it. These free electrons 
give rise to a complex refractive index for the propagating radio wave, 
resulting in its deviation from the original path, and also causing its 
attenuation. It was in fact, the experiments conducted to study the long- 
distance propagation of radio waves in the 1920 's that led to the discovery 
of the Earth's ionosphere. 

In these early experiments, different layers of the ionosphere which 
cause reflection of radio waves of different frequency bands were identified 
and were designated D , E, and F. Later studies on the ionospheric structure 
revealed that these layers could be identified with different regions, simi- 
larly designated, where the electron-density profile exhibits a plateau or a 
peak, as illustrated in Figure 1.1. It is sometimes possible to identify 
more than one layer within the same region, depending on the time of the day. 
For example, during the sunrise period a C layer, having a distinct peak 
develops in the lower part of the D region. During nighttime, the F region 
exhibits two distinct peaks known. as the FI and F2 layers . 


T 
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The present investigation is concerned with the D region, or the 
region below about 90 km of the ionosphere. This region exhibits several 
distinct features not found in the other regions of the ionosphere. The 
main reason for this distinction is the fact that the neutral atmosphere 
over this region has sufficiently high density so as to cause strong 
coupling bett^een itself and the D region. The features referred to above 
are: (1) high absorption of radio waves, (2) complexity of the ion chemis- 

try and (3) the occurrence of negative ions. 

The absorption of radio waves in the ionosphere at a given height 
depends on the product of the electron density and the collision frequency 
between the electrons and the neutral particles . The collision frequency 
is directly proportional to the atmospheric pressure, and hence decreases 
urith the height. The electron density, on the other hand, generally 
increases with height, at least in the V and E regions. The product of 
these two quantities, therefore, reaches a maximum, which happens to occur 
at 0- region heights. Thus, the absorption of high frequency radio waves 
propagating in the ionosphere takes place mostly in the D region. The re- 
ception of radio signals from long-distance transmitters is limited by this 
ionosphe"i>' absorption. Therefore, a knowledge of the behavior of the D 
region is important for the proper operation of high-frequency radio 
networks. 

The second feature mentioned above is really a consequence of the neutral 
chemistry of the region. The composition of the major constituents of the 
atmosphere up to the top of the D region is similar to that near the ground. 

In addition, a variety of active minor neutral constituents is produced in 
this part of the atmosphere as a result of the photodissociation of the 
major constituents as well as others that are transported from the ground. 


The primary ionization process produces electron-ion pairs. The posi- 
tive ion will either undergo dissociative recombination with the electron. 
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or react with other neutral species in a charge transfer process to yield 
a large number of new positive ions . Some of these positive ions can undergo 
three-body reactions producing complex cluster ions. These reactions become 
important only when the third-body density is high; the necessary density 
being found at ZJ-region heights and below. The D - region positive- ion 
chemistry is, therefore, much more complex than that of the higher regions. 

Finally, the high efficiency of the three-body reactions at these 
heights is also revealed in the formation of the negative ions. These are 
formed when a free electron gets attached to a neutral molecule through a 
three-body reaction. These ions also undergo a series of charge transfer 
and attachment reactions producing a variety of negative molecular and 
complex cluster ions. In the D region during daytime, negative ions become 
important below about 75 km, the negative-ion and electron concentrations 
becoming equal around 70 km. At heights below, the negative-ion chemistry 
controls the electron concentration. At nighttime, negative ions dominate 
most of the 0-region heights. 

1.2 Methods for D-Begion Investigations 

The methods available for investigating O-region ionospheric processes 
fall into two main categories: ground-based and rocket-borne. Until about 

a decade ago, when the rocket techniques were developed, studies on the 
ionosphere were totally dependent on ground-based radio sounding experi- 
ments. The ground-based methods can also be classified into two groups; 
the conventional methods which provide only an integrated picture of the 
D region, and the more recent methods which directly provide the electron- 
density values as a function of the height. The rocket-borne experiments 
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measure the electron densities as well as the positive-ion densities 
in situ 3 and are considered to give the most reliable picture of the D 

region at present. In this section, a brief review of these methods will 

* 

be given. 

1.2.1 Ground-based experiments. In these experiments, a radio 
wave is employed as a tool for probing the ionosphere. The propagation of 
a radio wave in the ionosphere is governed by its refractive index, which 
is a function of both the electron density and the wave frequency. When 
the electron density increases with height, the refractive index continues 
to decrease, causing the wave normal to deviate more and more from its 
original path. When the refractive index becomes zero, which happens when 
the plasma frequency of the medium is equal to the wave frequency, the 
wave begins to return to the ground. 

According to the simple ray theory of propagation, the frequency of a 
vertically incident wave that is reflected from a point having an electron 

_3 

density, N cm , is given by 

/ = 8.98>¥ kHz. (1.1) 

For a wave incident at an angle 0, this is multiplied by sec 9. The 
maximum frequency that gets reflected from a peak or a plateau in the 
electron-density profile is called the critical frequency of the corre- 
sponding layer. The electron density at the top of the D region is about 
10 4 electrons cm " , and hence the maximum frequency that is reflected 
from the D region is of the order of 1 MHz, when incident vertically on 
the ionosphere. 
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Low frequency and very low frequency propagation experiments were among 
the early experiments carried out to investigate the D region. Frequencies 
in the range 15-100 kHz are reflected from the low and middle portions of 
the daytime D region. Signals from transmitters generally used for navi- 
gational purposes were monitored at various distances, ranging from 100’ s 
to 1000* s of kilometers. The amplitude of the skywave as well as its phase 
relative to the phase of the ground wave was recorded. From these re- 
cordings, the diurnal behavior of the region near the reflecting level, 
and also effects such as the sunrise and sunset phenomena were studied. 

In the vertical incidence absorption measurements, the signals from a 
pulse transmitter, after reflection from the ionosphere, are received in a 
receiver located near the transmitter. Usually, frequencies ranging from 
about 1.8 to about 6 MHz are used. The received echoes are displayed on a 
cathode ray tube, and the amplitude and the delay are recorded. At night- 
time when the ^-region ionization vanishes and the absorption is low, 
multiple echoes are observed. The amplitude of these yield the calibra- 
tion constant of the system, or the unattenuated amplitude of the signal. 
Hence, the absorption during daytime is determined by deducting the observed 
amplitude from this calibration constant. The time delay gives the virtual 
height of the reflecting layer. 

Generally, the lowest frequencies, i.e., those between 1. 8-2.0 MHz 
return from heights in the range 95-100 km, well above the D region. 
However, most of the absorption takes place at P-region heights, between 
about 80 and 90 km, as mentioned before. Any changes in the observed 
absorption could, therefore, be attributed to changes in the electron- 
density values at these heights. However, the absorption, being an 
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integrated effect, cannot yield information about the fine structure of 
the D region. 

The experiments developed more recently to study particularly the 
D region are: (1) Partial- Reflection Method, and (2) Cross-Modulation 

Method. A third method, developed for studying the upper ionosphere and 
applied recently for lower ionosphere studies, is the Incoherent- Scatter 
Method. 

In the partial-reflection method, weak signals reflected from irregu- 
larities in the electron-density profile at O-region heights are recorded 
in both the ordinary and extraordinary modes [Gardner and Pawsey , 1953] . 

From the ratio of these amplitudes, the electron-density values are 
derived over the height range 65-80 km. In view of the very low intensity 
of these echoes, high-power transmitters and low-noise receiving sites 
have to be used. The finite pulse width, the errors in reading the echo 
amplitudes, and the long sampling periods are some of the factors that 
limit the accuracy of the electron densities derived from this method. 

Coyne and Betrose [1973] estimated the accuracy in the absolute values 
of the electron densities at any height to be about 50%, while the 
accuracy of any changes was estimated to be within ± 20%. 

The cross modulation of two radio waves in the ionosphere was long 
known, and was commonly referred to as the Luxembourg Effect. This method, 
first applied to ionospheric studies by Feger [1955], essentially consists 
of measuring the difference between the absorption of two downcoming waves, 
one passing through the normal ionospheric region, and the other passing 
through a heated or disturbed region. The heating or the disturbing of 
the ionosphere is done by a high -power transmitter whose signal is directed 
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to meet the downcoming radio wave at a required height. The high power 
beam raises the electron temperature of the plasma, which in turn, modifies 
the collision frequency. The resulting change in absorption, which is of 
the order of 1-2%, is measured. From these data, both the electron 
density as well as the collision frequency could be determined as a 
function of the height [Throne et ai. 3 1968; Lee and Ferraro, 1969]. 

Incoherent-scatter facilities have been in existence for a number 
of years to study the ionospheric E and F regions. In view of the 
variety of information this method provides, it has proved to be a power- 
ful tool for investigating the middle and upper regions of the ionosphere 
[Bowles, 1958]. In the original systems, the presence of int erf erring 
ground clutter signals has made it impracticable to obtain useful signals 
from the lower ionosphere, particularly below about 150 km. Recently, 
Armistead et al . [1972] have shown that, with certain system modifica- 
tions, it is possible to bring down the heights of measurement to about 
80 km. However, this technique has not been perfected yet to provide 
D - region profiles with sufficient accuracies. At present, results could 
be considered reliable down to about 90 km only. 

1.2,2 Rocket experiments. In situ measurements of ionospheric 
parameters by rocket experiments have been widely used in recent years. 

Of these, the experiments designed for 0-region studies can be divided 
into two groups: (1) radio propagation experiments and (2) ion-probe 

experiments. In radio experiments, the generalized magnetoionic theory 
is employed to deduce the electron densities, while in ion-probe experi- 
ments, the total probe current gives the electron densities. 

The radio experiment essentially consists of sending two signals, one 
in the ordinary mode and the other in the extraordinary mode, on the same 
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frequency from ground-based transmitters, and receiving them on a rocket- 
borne receiver. The two signals are recorded separately and telemetered 
to the ground \Mechtly et at . , 1967]. The Faraday rotation or the 
rotation of plane of polarization of the wave containing both the ordinary 
and the extraordinary modes is also measured and telemetered to the ground. 
From the two independent signals, the differential absorption is measured. 
By combining these two measurements, the collision frequency is determined. 
This value, in turn, is used with either of the two measurements to yield 
electron densities. The Faraday rotation method gives reliable electron 
densities at the upper portion of the D region, while the differential 
absorption gives reliable results in the lower half of the D region. Both 
methods give overlapping electron densities in the central part. In view 
of the approximations made in the formulations, iterative techniques are 
employed to arrive at the final results. Recently, the accuracy and the 
height resolution of the measurements have been improved by employing two 
independent propagation experiments operating on two different frequencies 
[Meahtlyj 1974] . 

Propagation experiments by themselves do not yield high resolution 
data to reveal the fine structure of the electron- density profiles. 

Usually, in the measurement of absorption below 90 km, large fluctuations 
occur in the telemetered data. The differential absorption is therefore 
determined by employing a sampling technique, and averaging over a selected 
range of points. This results in poor height resolution in the final 
electron densities. 

Ion probes which respond to the electron or ion concentrations in the 
immediate neighborhood of the probe are widely used in rocket experiments. 
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This method has the advantage over the propagation experiments in that 
it reveals the fine structure of the electron-density profiles, but certain 
types lack the capability of yielding absolute values of electron concen- 
trations. Though several types of in situ probes are available for 
ionospheric investigations [Sayers A 1970], only the following used 
specifically for D region studies will be discussed here: (1) Langmuir 

probe, (2) Gerdien condenser, and (3) electrostatic probe. 

The Langmuir probe technique was developed for ionospheric studies 
by Smith [1964]. This consists of a metallic electrode mounted axially 
on the nose tip of a rocket. The probe is biased at a fixed potential 
relative to the rocket body and it collects either the positive or the 
negative ions depending on the polarity of the bias potential. The probe 
current is directly proportional to the concentration of the ions or 
electrons collected, the constant of proportionality being a function of 
the probe shape and size. In view of the uncertainties involved in deter- 
mining this parameter, direct measurement of the ion or electron concentra- 
tion by this method is not possible. Therefore, it is customary to convert 
the collector current into electron concentrations with the help of radio 
propagation experiments or ionosonde data. The most reliable D-region 
electron densities have been obtained by calibrating the Langmuir probe 
current profiles with the electron densities deduced from differential 
absorption and Faraday rotation experiments described previously. 

The Gerdian condenser has been generally used to measure the positive- 
ion concentrations. In the upper part of the D region where the abundance 
of negative ions is negligible, the electron density could be obtained 
from this measurement. This is essentially a pair of two concentric 
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metallic cylinders located axially at the nose tip of the rocket. The 
passage of the rocket through the plasma maintains a flow of ionized air 
through the two cylinders. An electric field applied across the annular 
region between the two cylinders causes a flow of electric current through 
the condenser. When the applied electric field exceeds a certain value, 
this current reaches saturation. At this saturation mode, the current 
collected is equal to the product of the charge density of the plasma, 
aperture area, electronic charge, and the effective velocity of ion flow 
into the condenser. The velocity is deduced from the rocket speed and its 
inclination to the trajectory. Thus, from the current measurements, charge 
density can be determined directly. At region heights where the mean 
free path of the particles is small compared to the dimensions of the con- 
denser, reflection of particles take place at supersonic speeds. At 
altitudes above 70 km, this effect can cause errors up to about 25% [ Kane, 
1972], which can become serious at lower altitudes. Further, at altitudes 
above about 85 km, vehicle potential can cause errors up to a factor of two 
in the measurement of positive- ion densities. Thus, the Gerdien condenser 
is best suited for measurements in the middle part of the D region. 

Sagalyn and Smiddy [1964] developed a technique for measuring inde- 
pendently the positive and negative charge densities of a plasma. Their 
instrument consists of two spherical electrostatic probes, one for the 
collection of positively-charged particles, and the other for the collection 
of negatively- charged particles, enclosed in outer spherical grids. The 
collector is maintained at a fixed bias voltage relative to the outer grid, 
which is maintained at the vehicle potential for the positive- ion probe, 
and at a slightly positive value for the negative-ion probe. Under normal 
operation, the fixed bias voltage collects all the charged particles 
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entering the outer grid, and the resulting current flow is measured. To 
evaluate the corrections due to vehicle potentials, the outer grid voltage 
is swept linearly and the resultant currents are measured. From a knowledge 
of the dimensions of the probes and pre- determined parameters, the total 
ion densities are evaluated in terms of the collector currents. 

In addition to the above probes which measure the total ion or electron 
densities, mass spectrometers are also flam in rockets to measure the dis- 
tribution of different ionic species in the ionosphere. With the develop- 
ment of this technique, a complete new picture of the D region has emerged. 
The presence of complex cluster ions in the lower part of the D region was 
first detected by Nca>oisi and Bailey [1965] who flew a quadrupole mass 
spectrometer in a rocket. 

This instrument has four rods mounted axially at the four corners of a 
square. The ions are dram in through an aperture at one end, while the 
ion detector is located at the opposite end. Suitable RF and DC potentials 
are applied across these rods which make it possible for an ion of only a 
certain charge-to-mass ratio to reach the ion detector. Sweeping of the 
RF voltage results in the collection of ions of different charge-to-mass 
ratios sequentially. In order to make the trajectory time smaller than 
the collision period, the spectrometer chamber is evacuated continuously 
and maintained at liquid nitrogen temperatures. The collector currents 
are proportional to the respective ion densities. The absolute values are 
obtained by conducting an independent experiment which measures the total 
ion density, using a different method such as the electrostatic probes 
described before. Inaccuracies in reading the current peaks, drifts in 
the DC amplifier could result in overall errors up to about 50%. 
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1.3 Evidence for Diurnal Variation 

Both the LF and VLF recordings and the absorption measurements pro- 
vided the first evidence of the strong solar control of '•‘he D region. 
[Bracewell et al . 3 1951; Appleton and Piggott 3 1954], Figure 1.2 illus- 
trates a typical recording of the phase variations obtained more recently 
by Belvose [1963], The phase variations could be attributed to variations 
in the apparent height of reflection. For a fixed path and frequency, the 
apparent height corresponds to a level having a given electron density. 
Hence, these curves in effect show the displacement of the electron-density 
profile near the level of reflection. For the shorter path, which corre- 
sponds to a higher level of reflection, the variation between sunrise and 
sunset is rather smooth. For the longer path, the variation is abrupt 
near the sunrise and sunset, and in between, there is no significant varia- 
tion. These results indicate that the upper portion of the D region has 
a strong solar control, while the lower region has virtually no solar 
control during the daytime hours. 

In summer, the daytime variations over short paths are so regular and 
constant, that the observed variations can be represented by the expression 

h = h * A log (sec x] » (1*2] 

X & & 

where h and are the apparent heights when the solar zenith angle is 
X and o, respectively, and A is a constant for the day. A closer look 
at the results show that the curves are not quite symmetrical about the 
local noon. This becomes apparent when h values are plotted against log 
(sec x) * The straight lines passing through the morning values do not fall 
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on the straight line connecting the afternoon values. There is a delay of 
about 20-30 minutes in the appearance of the minimum points in these curves. 

Vertical-incidence absorption measurements have been carried out in 
the past over long periods of time, both at mid-latitudes and low latitudes 
[Appleton and Piggott 3 1954; Gnanalingam 3 1974] . These measurements show 
that the absorption varies linearly with the sunspot number, has a diurnal 
variation, seasonal variation, and also a day-to-day variation. The diurnal 
variation of absorption is usually expressed in the form 

L = L o cos n x> Cl -3) 

where L and L . are the absorption values at solar zenith angles x and o 3 
respectively; and n is a constant. The simple ray theory predicts that 
the index n for the diurnal variation of the total nondeviative absorp- 
tion for waves traveling in a Chapman layer is 1.5. However, observations 
show that n takes the value 0.75 for mid-latitude stations [Beynon and 
Davies^ 1955]. Using the method of oblique incidence absorption, Sehwentek 
[1966] also found that n = 0.75 on most of the quiet days. On the other 
hand, Gnanalingam [1974] found that, for a low latitude station, n takes 
the values 0.8 and 0.9 for the equinoctial months at solar maximum and 
solar minimum, respectively. These results have important implications’ 
on the diurnal behavior of the Z)-region ionization. 

The plots of absorption against log ( cos also exhibit the diurnal 
asymmetry. The morning points lie on a straight line having a smaller 
gradient than that joining the afternoon points. This asymmetry is again 
of the order of 20-30 minutes. 
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Coyne and Belrose [1972] used the partial-reflection system to measure 

both the diurnal and seasonal variations of the electron densities between 

65 and 81 km. They found a strong solar variability in these measurements, 

at least up to a solar zenith angle of about 75°. Here again, the variation 

was found to have a slight asymmetry about local noon. 

Using both the partial -reflection and cross-modulation techniques, 

Throne et at. [1968] measured the diurnal variation of the electron 

densities between 65 and 85 km at a mid-latitude station during quiet solar 

conditions. Throne [1969] attempted to interpret these measurements using 

a ZJ-region model where the electron losses are due to a recombination 

. 2 

process, for which the production rate q is proportional to N . 

For a Chapman-type ionosphere, consisting of a single ionizable con- 
stituent ionized by monochromatic radiation, which is attenuated by a 
single absorbing gas, the rate of production q is related to x according 
to 


log g <7(3 jX) s A - Biz) sec x, Cl-4) 

where A is a constant, and 

Biz) is a function of the height s only. 

From the measured values of N, the production rates were deduced using 
known electron loss rates . The log q values were then plotted against 
sec x for different heights. It was observed that the gradients of these 
plots were much steeper than the values of Bis) calculated for the cor- 
responding heights. According to Thrane, this discrepancy could be 
explained if the ionizable constituent or the electron loss rate has a 
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solar zenith angle variation, or if the dominant ionization is not produced 
by a single gas. 

Hang and Throne [1970] have re-examined these measurements in the 
light of a new D - region model suggested by Haug and Landmark [1970], In 
this model, electron losses take place through an attachment- like process. 
Consequently, the production rate q is proportional to N under quasi- 
equilibrium conditions. Haug and Thrane compared the plots of log „V versus 
sec x made on the basis of this theory, with the experimental plots, and 
found that the discrepancy is smaller with this model than with that used 
by Thrane. In view of the uncertainties involved in the measurements, as 
well as in the other parameters used, this small discrepancy does not 
appear to be significant. 

The diurnal study of the ionosphere by rocket experiments is a diffi- 
cult task due to the high costs and various operational problems. Hence 
rocket experiments are not generally employed to study diurnal variations 
of the ionosphere. Nevertheless, Meehtly and Smith [1970] were successful 
in making four sets of reliable measurements of the electron- density 
profiles at four zenith angles on the same day (Figure 1.3). The extent of 
the solar control of the D -region ionization becomes evident from this set 
of measurements. 

The two large angle profiles were measured during the sunrise hours. 
The x = 18° profile corresponds to the noon. To study the variability of 
the D-region ionization during the daytime, only the two profiles at 
X = 18° and 60° are of interest. 

1.4 Some Outstanding Problems of the D Region 

Much progress in our understanding of the D region has been made 
during the past few years, though some problems still remain unsolved. The 
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ELECTRON DENSITY (cnr 3 ) 

Figure 1.3 Rocket measured D region electron-density profiles at solar zenith angles 90° 
18°, and 60° on July 24, 1968 [Mechtly and Smith, 1970] . 
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present status of the P region has been the subject of two recent reviews 
by Sechrist [1972] and Thomas [1974]. Some salient features of the P 
region, its production and loss processes will be given here in order to 
identify some of these outstanding problems. 

A large number of P-region electron-density profiles have been 
measured during the daytime at Wallops Island (latitude = 40°N) using 
rocket-borne propagation and Langmuir-probe experiments [Meehtly et al. 3 
1972a]. All of these profiles were measured at the same zenith angle, viz 
60°, and during different seasons and solar conditions. Yet, almost all 
of them exhibit some consistent characteristic features. 

Between about 82 and 88 km, the electron density increases very 

rapidly with altitude. In most cases, this increase is more than an order 

of magnitude and takes place within a height interval of about 2 km. 

Above 90 km, the increase is rather gradual, varying between 10 4 and 10 5 
-3 

electrons cm in the height range 90-105 km. A second ledge appears 
around 60 km. This is more prominent in the active-sun profiles. In the 
quiet-sun profiles, this ledge is absent, but they show a marked plateau 
region between 60 and 80 km, where the increase in the electron density is 
very slow. 

These features of the electron- density profiles are closely associated 
with some features present in the ion distributions. Direct measurements 
of the P-region positive-ion composition obtained from rocket -borne mass- 
spectroscopic experiments reveal (Figure 1.4) that heavy cluster ions of 
the type . (n = 1, 2, ...) become dominant below about 82-85 km 

level, and that molecular ions NO - and C> 2 + dominate the altitudes above 
[Narcisi and Bailey 3 1965; Goldberg and Aikin 3 1971; Krankousky et al . , 
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1972]. This sudden drop in the hydrated cluster ion concentration and the 
sudden increase in the electron density at the same height range has been 
interpreted as due to hydrated cluster ions having a higher recombination 
coefficient than that of the molecular ions [Sechrisb, 1970], Though 
laboratory measurements of the recombination coefficients of both the 
molecular ions and the cluster ions are available, uncertainties in the 
distributions of such minor constituents as I^O, 0, NO, H, etc. which are 
important in the reaction schemes yielding these cluster ions, have made 
it difficult to interpret the observations quantitatively. Attempts by 
several workers to reproduce theoretically the sudden disappearance of 
the cluster ions and the sharp ledge in the electron densities near 82- 
85 km have so far failed [Reid, 1970; Goldberg and Aikin, 1971; Hunt, 

1971a]. 

More than Ehe sharp drop of the cluster ions, the prediction of the 
actual type of the dominant cluster ion present in the D region by itself 
has remained an outstanding problem [Reid, 1971] . The known values of NO 
and the ionization rates would yield water clusters of N0 + around 80 km, 
rather than those of the type H + * (H_0) which have been observed to 
dominate these altitudes. Attempts to solve this problem have not been 
successful so far, due to the lack of proper knowledge of either the rate 
constants involved or the concentrations of the minor neutral species, or 
due to both. 

The other problem is concerned with the N0 + production rates in the 

D region. The strong solar Lyman-a line at 1216 A ionizes the neutral NO 

present in the mesosphere producing N0 + . Direct measurements of NO con- 

7 8 - 3 

centration show that it varies between 1.5 x 10 and 1 x 10 cm in the D 
region with a minimum near 85 km [Meira, 1971] . Near 80 km, or immediately 
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below the ledge in the electron-density profile, the major source of elec- 
tron production is this ionization of NO, Through various ion-molecular 
reactions , N0 + produced will yield a variety of cluster ions including 
the water clusters that are dominant at these altitudes. Under quasi- 
equilibrium conditions, the final electron density values will be deter- 
mined by the rate of recombination of these dominant ions and the 
electrons . 

The ion-pair production rates have been calculated using known values 
of fluxes, ionization cross sections. The rocket -measured electron densi- 
ties are then used to calculate the effective recombination coefficient, 
defined as 


“eff = Cl.S) 

Calculations show that near 80 km, « e ££ has to be greater than 
™S 3 -1 

1 x 10 cm sec in order to balance the production and loss rates 

[j Donahue, 1972] . Recent estimates by Mechtly et al. [1972b] using rocket 

observations of the changes in the electron-density profiles during 

eclipses give values about 5 x 10 cm sec for heights immediately 

below the ledge in the electron-density profiles. However, the labora- 

—6 3 

tory measurements of Leu et at. [1973] give values about 3 x 10 cm 
sec - ", almost an order of magnitude less than the values deduced from 

4 

ionospheric observations for the recombination coefficient of H 
which is the dominant ion near 80 km. 

This suggests that either the production rates used are too high, or 
the actual dominant ion below the ledge is something other than H + • (H^O ) 2 * 
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According to Donahue [1972] , a low production rate can be realized if a 
reduced ionization rate is used for NO above 75 km. On the other hand, 
the presence of a heavy cluster ion having a larger recombination coef- 
ficient than that of H + *(H 2 0) 2 could be accepted if it is assumed that 
these heavy clusters get fragmented during the process of rocket measure- 
ments, resulting in the recording of a lower mass number than the mass 
number of the ion actually present in the ambient D region. 

Another problem that has drawn relatively little attention is the 
production rates of 0 2 + between 70 and 90 km. Above 90 km, X-rays between 
30 and 100 A ionizing both 0 2 and N 2 , and the Lyman-B radiatiop ionizing 
0 2 , produce C> 2 + . Below 70 km, galactic cosmic rays ionizing both 0 2 and 
N 2 produce the necessary 0 2 + . In between, there are several sources pro- 
ducing ' a H of them being minor sources compared to NO ionization. 
X-rays between 2 and 10 A ionizing both N 2 and 0 2 > solar EUV radiation 

ionizing ) and precipitating electrons with initial energies greater 

* Q 

than 40 keV, all produce 0 2 + . Near equatorial latitudes, the precipitating 
electrons are not considered important. The theoretical calculations of 
Goldberg and Aikin [1971] show 0 2 + concentrations rapidly falling off 
below 82 km, while their measurements show 0 2 which remains almost 
constant between 75 and 82 km. Th? indicates that either the precipi- 
tating electrons are important near the equatorial latitudes where the 
measurements are taken, or some other source of 0^ + is present in this 
height range. Comparing the 0 2 + production rates due to EUV photo- 
ionization of 0-C^A ] with the results of Naroisi and Bailey [1965], 

Huffman et al. [1971] have also come to a similar conclusion, Horton and 
Beid [1972] have examined the possibility of such an additional source of 
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O 2 , without success. A source of 0 ^ below the ledge can also solve, at 
least partly, the problem of the formation of the water- cluster ions. 

It will be possible to identify the origin of the ionizing source 
if a diurnal study of the 0 ^ distribution could be made. A source such 
as precipitating electrons has no diurnal variation, whereas one with solar 
origin will show a diurnal variation, at least over the height range where 
the optical depth for the ionizing radiation is comparable to or larger 
than unity. 

Thus, a realistic model of the D region should be able to explain: 

1) the sharp ledge in the electron- density profile, 

2) the formation of the dominant water-cluster ions, 

3) the balance of production and loss rates, and 

4) the ©2 concentrations between 70 and 90 km. 

1,5 Objectives and the Outline 

The observation of the variation of ionization with the solar angle 
will give information as to the height at which sources of different 
origins, such as Lyman-a and precipitating electrons become important, 

Such a deduction is straightforward, however, only if the incident radia- 
tion alone is a function of the solar zenith angle. Departures from this 
behavior will show that others, such as the concentration of ionizable con- 
stituents and those responsible for the electron losses could also vary 
with the zenith angle. Thus, an understanding of the phenomena occurring 
in the D region requires an understanding of the behavior of the neutral 
constituents at different solar zenith angles, and also of the variations 
in the ionizing radiations. In this report, an attempt will be made to 
construct a unified diurnal model of the neutral atmosphere between 30 and 
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120 km, and the D region. The neutral atmosphere will contain the major 
constituents such as and 0 2 , as well as the minor constituents that arc 
either transported from the ground, or produced through photodissociation 
and photoionization. The final objective is to compare the calculated 
electron-density profiles with those measured by rocket experiments at 
X = 18° and 60° on July 24, 1968 at Wallops Island [ Medhtly and Smithy 
1970]. Therefore, the input data taken for this study will correspond 
to this date and location. The emphasis will be for the region where 
the negative ions are not important. 

In Chapter 2, the photochemistry of the minor neutral constituents, 
including their photodissociation rates and the production- loss mechanisms 
of different groups of constituents are discussed. The calculation of the 
initial values using the steady- state transport and photochemical equa- 
tions for use in the time -dependent transport model are presented in 
Chapter 3. The solar zenith angle variation of the minor neutral con- 
stituent concentrations is also investigated in this chapter. The ion- 
pair production rates responsible for daytime D-region ionization are 
calculated in Chapter 4. In Chapter 5, the production and the loss 
processes, as well as the measurements of the composition of both the 
positive and negative ions are discussed. The results of the model cal- 
culations describing the solar zenith angle variation of the positive- 
ion composition, negative-ion composition, and the electron densities are 
given in Chapter 6. Chapter 7 presents the summary and the conclusions. 
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2, PHOTOCHEMISTRY OF NEUTRAL SPECIES 

The ionosphere is formed as a result of the ionization of neutral 
molecules present in tl e upper atmosphere, by solar and other radiations. 
Therefore, a study of the ionosphere should necessarily include the study 
of the behavior of the neutral constituents. If one begins with a simple 
model of the atmosphere having only the two major constituents, molecular 
oxygen and nitrogen and subject them to the effects of solar radiation, 
the development of it into its present form can be followed up through 
different stages. 

This chapter begins with a discussion of the neutral atmosphere model 
used in the present study. The calculation of the photodissociation co- 
efficients of various constituents will be next given, followed up by a 
discussion of the photochemical production and loss processes of different 
groups of constituents. 

2.1 Neutral Atmosphere 

From direct measurements of the density of air it has been established 
that the neutral atmosphere composition remains unchanged from the ground 
level up to a height of about 100 fan. Below this level the atmosphere is 
considered to be well mixed maintaining constant mixing ratios for those 
constituents that are not chemically active. At higher altitudes, the dis- 
tribution of each constituent is governed by diffusive equilibrium. 

The ground level composition adopted here is taken from the model of 
Jacohia [1971], and is given In Table 2.1. 

According to Jacchia's model, molecular oxygen begins to photo- 
dissociate at 90 km and form atomic oxygen. The mean molecular weight 
is defined by a sixth order polynomial, from which the distribution 
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of the individual species is calculated. The hydrostatic equations using 
the scale heights corresponding to each constituent are used to calculate 
the number densities above 100 km. 

The temperature appears as an independent parameter in Jacchia's 
model. Up to 90 km, the temperature profile is constructed from three 
c^ts of direct measurements carried out on July 24, 1968 at Wallops Island 
[Smith et at., 1970]. The temperatures above 90 km are determined by 
empirical formulas whose coefficients depend on the 10.7 cm solar flux, 
latitude, longitude and the time of day. For the purpose of this study, 
a 10.7 cm flux of 150 units (10 ^ W m ^ Hz~*) and local noon (x = 17.8°] 
corresponding to Wallops Island geographic coordinates were used. The 
density and temperature profiles thus obtained are shown in Figure 2.1. 

The calculation of the distribution of minor neutral constituents 
involves the simultaneous solution of several continuity equations which 
are in the form of parabolic type partial differential equations. The 
solution of such equations call for the imposing of boundary conditions. 

For many constituents these boundary conditions are not 11 known, and 
in some cases, suitable values have to be determined by initial trial cal- 
culations. The departure of the solution from the exact situation due to 
the use of incorrect boundary values will be largest near the boundaries. 
Therefore, in order to minimize such errors within U-region heights, the 
boundaries are extended up to 120 km at the top, and down to 30 km at the 
bottom. 

2.2 Photodissociati on Coefficients 

The amount of radiation energy absorbed by a single molecule of a con- 
stituent (denoted i) at a given altitude over a wavelength interval AA at 


ALTITUDE (km) 


DENSITY (gm cm' 3 ) 
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Figure 2.1 The temperature and density profiles of the neutral atmosphere over Wallops 

Island on July 24, 1968. Values up to 90 km are based on rocket measurements 

by Smith et at. [1970] and above 90 km, on empirical formulas given by 
Jaachia [1971] . 
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a wavelength X is given by 

W, (X.X) = J(X,AX,s) c-CX), (2. ID 

where J(X,AXjs) is the intensity of solar radiation at the altitude z 
and wavelength X taken over the wavelength interval AX, and a. (X) is 
the absorption cross section of the ith constituent at the wavelength X. 

The intensity of solar radiation remaining at a given height depends 
on the amount absorbed by the column of air above that height measured 
along the ray path, and is given by 

JfX^AXjS) = J^CX.AX) expC-r (Xj 3 DD (2.2) 

where T^X, AX) is the solar radiation flux integrated over the wavelength 
interval AX incident on top of the atmosphere, and tCxj^D is the optical 
depth factor. 

The ability of the solar radiation to penetrate the atmosphere for 
different wavelengths is expressed in terms of the unit optical depth, 
defined as the altitude where the optical depth factor is unity. It 
follows from equation (2.2) that at this altitude, the fractional attenua- 
tion of the incident radiation is equal to 1/s. Figure 2.2 illustrates 
the penetration depths for radiation incident normally on the atmosphere, 
as a function of the wavelength up to 3200 A. It is seen from this curve 
that wavelengths above 1000 A could reach heights below 120 km. Between 
1800 and 2000 A, there is a sharp drop in the optical depth, from a height 
of 80 km to 40 km. This wavelength interval corresponds to the Schumann- 
Runge bands in the 0 2 absorption spectrum. Beyond 2000 A up to about 


PENETRATION DEPTH (km) 



Figure 2.2 The penetration depths of solar radiation incident normally on the atmosphere. 63% of 
the incident energy in a given wavelength interval is lost above the altitude shown. 
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3000 A, the penetration depth lies more or less constant in the 30-40 km 
interval. By comparing this curve with the absorption spectra of oxygen 
and ozone, it becomes apparent that the penetration of wavelengths below 
about 2000 A is controlled by absorption due to oxygen, while ozone 
absorption controls the penetration of longer wavelengths. 

In determining the optical depth factor, it is therefore necessary to 
consider only the absorption by oxygen and ozone. It is calculated using 
the formula 


(XjS) = Y* 0 (*) J [0 2 ]e?3 + a 0 (X)J [0 3 ]<23|secx 


(2.3) 


where [0 2 ] is the number density of oxygen, and 
[0 3 ) is the number density of ozone. 

The photodissociation coefficient, is obtained by summing up 
the energy absorbed over the entire spectrum effective in dissociating 
the given constituent, and is given by 


= I J 0 CX,AX) 0^0) exp(-x(x,3)) C 2 * 4 ) 

X 

The height dependence of J . is determined solely by the optical 
depth factor. For small values of t, it is almost independent of the 
height. In calculating the optical depth factor for wavelengths where 
absorption by ozone is important, it is necessary to assume initially the 
height distribution of this gas. The adopted distribution, which is based 
on direct observations, is given by 

[0 3 ] = 10 16 - 3/1 ° C2.S) 


for heights above 30 km. 
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2.2,1 Solar radiation in the 1000-4000 S region . The ranges of wave- 
length that are of interest in photodissociation studies of the neutral 
constituents lies between 1000 and 4000 A. This is mainly a continuous 
spectrum, being of thermal origin, except for the lines that appear towards 
the 1000 A end. Of these, the more intense lines are the Lyman-6 at 
1025.7 A and Lyman-a at 1215.7 A. However, for the purpose of calculating 
the photodissociation rates, only the latter is of importance. This line 
has been under investigation over several years, and it has been observed 
that it is a fairly stable line, with the intensity variation at most of 
the order of 1.5 between the solar maximum and minimum [Weeks, 1967]. The 
intensity of this line has also been measured using satellite-borne instru- 
ments in the recent past. Timothy and Timothy [1970] have measured the 
intensity variations of both the Hydrogen I Lyman-a [1215.7 A] and Helium 2 
Lyman-a (304 A] lines from 1967 to 1969. The intensity of 1215.7 A lines 

measured on July 24, 1968, the day that is of interest here, is 3.7 x 10 11 

-2 -1 11 
photons cm sec . The average for the whole period is 3.6 x 10 photons 

-2 -1 11 

cm sec . This value, however, is slightly higher than 3 x 10 photons 

-2 -1 

cm sec , the value reported by Vidal-Madjar et dl. [1973] as the average 
for the years 1969 and 1970, monitored on board the OSO 5. 

The flux values used in the region 1027-1310 A are those recom- 
mended by Hinteregger [1970] for medium solar activity. This region 
includes a few Si and 0 lines and a continuum which contains about 1/5 
the total energy (excluding the Lyman-a line] in the band. 

The fluxes between 1300 and 2000 A used in aeronomical calculations 
in the past were taken from the observations of Detwiler et at. [1961]. 
Recent measurements of the spectrum between 1400 and 1875 A using photo- 
electric instruments, carried out on 1968 September 24 (10.7 cm flux = 
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143 units) by Parkinson and Beeves [1969], however, give values much less, 
about 1/3 to 1/5 of those reported by Detwiler et al. Considering the 
reliability of the experimental technique employed by Parkinson and Reeves, 
Hinteregger [1970] suggests that the values quoted by Detwiler et al. 
for the above range of wavelengths should be scaled down by a factor of 3. 
In the present, calculations, therefore, this correction was applied over 
most of this region. Towards the two ends, this factor was made to reduce 
gradually down to unity in order to join smoothly with the rest of the 
spectrum. 

The region between 2000 and 4000 A was measured by Tousey [1963] as 
early as 1960, the intensities measured being averages over 50 A intervals 
from 2000 to 2600 A, and over 100 A intervals from 2600 to 4000 A. More 
recent measurements of Bonnet [1968] generally agree with Tousey's 
measurements. Slight disagreement over small wavelength intervals could 
be attributed to averaging effects by the measuring instruments. The 
spectral distribution expressed in photons cm sec A over the entire 
region 1000-4000 A, adopted in the present calculations , is shown in 
figure 2.3. 

2.2.2 Absorption cross sections. The absorption cross section of 
a constituent, which is a function of the wavelength, gives a measure of 
the radiation absorbed by the constituent resulting in its dissociation. 

In general, the absorption of energy need not necessarily give rise to dis- 
sociation; instead the molecule which absorbs the energy could be raised 
to a higher excitation level. However, in the present work, it is assumed 
that the absorption of energy by a molecule results only in its 


dissociation. 
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The solar spectrum between 1000 and 4000 
Detmiler et al . [1961] , Parkinson o' 
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As mentioned previously, wavelengths below 2000 A are absorbed mainly 
by oxygen while wavelengths above are absorbed by ozone. Hence, the 
absorption cross sections of these two constituents control the photo- 
dissociation coefficients of all constituents, and their values adopted 
in this study are illustrated in Figure 2.4. Of the wavelengths below 
1300 A, only the window at the Lyman-a wavelength is of importance. The 
presence of this atmospheric window allows the strong Lyman-a radiation 
to reach 0-region heights resulting in the photoionization of NO, which 
is the major source of 0-region ionization. This line is also important 
in the photodissociation of constituents such as H^O and CH^ which have 
large absorption cross sections around this wavelength. High resolution 
measurements of the absorption cross section of 0^ at the bottom of the 

window were carried out by Ogatia [1968], He finds that the minimum of 

- 20 2 

the transmission window, having a value 0.945 x 10~ cm , lies at 1216.0 A, 

-20 2 

while the value at the center of Lyman-a line is 1.05 x 10 cm . There 

-20 2 

fore, a value of 1.00 x 10 cm is taken as the mean absorption cross 
section of oxygen at the Lyman-a line. 

The cross sections for wavelengths in the Schumann- Rung e continuum 
(1950-1750 A) and up to 1900 A were measured by Watanabe et at. [1953]. 

This region shows a peak between 1400 and 1500 A, and a sharp continuous 
drop up to about 1750 A. Beyond this wavelength, up to about 2000 A, the 
spectrum consists of a large number of rotational bands (Schumann-Runge 
bands). The maxima of the rotational lines in each of the vibrational 
bands lying in the region 1760-1900 A have been measured by Ackerman 
et at. [1969]. The continuum (Herzberg continuum) in the 1850-2600 A 
interval, underlying the Schumann-Runge bands, has been measured by 
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Ogama [1971], and it is found that this continuum lies more than 100 
times below the peak values reported by Ackerman et al. In spite of these 
measurements, a proper knowledge regarding the detailed structure of these 
bands is still lacking, and this has caused much uncertainty is estimating 
the extent of penetration of solar radiation down to lower mesospheric 
heights. Even if such measurements are available, the calculation of the 
contributions due to radiation penetrating through each of these windows 
would be a difficult exercise in view of their fine structure. 

To get over this problem, Hudson et at. [1969] measured the trans- 
mission factor for radiation of different wavelength bands passing through 
a column of oxygen. These results can be conveniently used to calculate 
the photodissociation rates in the corresponding wavelength bands. Mole- 
cular oxygen actually undergoes predissociation at these wavelengths, and 
the photodissociation coefficients of 0 2 include the contributions from 
this effect. 

The absorption cross section of many of the constituents, with the 
exception of 0 3 , NO^, HC> 2 , ri 2 0 2 and HNO^ fall off below 2000 A. For these 
constituents the photodissociation coefficients are controlled solely hy 
0 2 absorption, and their values vanish off before reaching lower heights. 
Since the solar fluxes below 2000 A are much less than those in the higher 
wavelengths, the photodissociation coefficients of these constituents 
generally have small values even at zero optical depths. 

For the constituents mentioned above, however, the large flux above 
2000 A, and the smaller attenuation due to ozone absorption result in 
higher photodissociation coefficients. Further, they maintain almost 
constant values down to about 50-40 km, where they begin to get attenuated 
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due to the presence of the ozone layer. In fact, the second term in 
equation (2.3) is important only for these constituents. 

The ozone absorption spectrum is illustrated in Figure 2.5. Values 
betiireen 1100 and 2000 A are taken from Tanaka et at. [1953] , and for the 
region above 2000 A, values are taken from Inn and Tanaka [1959] . The 
ozone spectrum exhibits a peak in the Hartley continuum near 2600 A. 

Above 3000 A, the cross section decreases with the wavelength allowing 
the wavelengths in near ultraviolet and visible regions to reach the ground. 

The photodissociation of H 2 0 is most important in the chemistry of 0-H 
constituents, and also in estimating the residual H^O in the mesosphere. 

Tile cross sections used for H 2 0 are taken from Watanabe and Zelikoff [1953] 
and Thompson et at. [1963]. In the SR continuum, this curve reaches a peak 
near 1650 A, and drops off rapidly beyond about 1800 A. Since large values 
of q occur in the region of SR bands, any errors in estimating the pene- 
tration depth of the solar radiation at these wavelengths could cause a 
certain amount of uncertainty in the corresponding photodissociation co- 
efficients, Anderson [1971a] calculated the photodissociation coefficients 
of H 2 0 using the transmission factors determined by Hudson et at. [1969], 

The values obtained by him are about a factor of 10 higher than the values 
given by conventional data for similar solar zenith angles at 60 km. In 
the present work, the cross-section curve over the SR bands was modified 
to give H 2 0 dissociation rates close to those obtained by Anderson. 

Of the other 0-H constituents, both H0 2 and H 2 Q 2 can get photo- 
dissociated producing OH molecules. Due to non-availability of absorption 
data for H0 2 » it has been customary in the past to assume the same values 
for H0 2 cross sections as those measured for H 2 0 2 . However, recent measure- 
ments of the absorption cross sections of H0 2 by Paukert and Johnston 
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[1972] have enabled independent calculation of the photodissociation 
rates for this constituent. These new data give photodissociation rates 
about six times higher than those of H 2 0 2 zero °P^* ca ^ depth. The data 
used for H 2 0 2 are those quoted by Paukert and Johnston. 

Cross sections used for N0 2 are taken from the work of Nakayama et al. 
[1959] for wavelengths between 1080 and 2700 A, while the values of Hall 
and Blaeet [1952] are used for higher wavelengths. In the case of N 2 0, 
absorption cross sections measured over the region 1080-2100 A by Z elikoff 
et al. [1953], and values quoted by Bates and Hays [1967] for higher wave- 
lengths are adopted. The absorption spectrum of NC> 2 maintains a high value 
even above 2000 A, while the N 2 0 spectrum falls off rapidly after reaching 
its peak near 1800 A. Consequently, the photodissociation rates of N0 2 are 
generally much higher than those of NjO. Cross sections used for HN0 3 are 
taken from the recent work of Johnston and Graham [1973], 

Recently Strobel [1971b] suggested that 6(0,0), 6(1,0) and 6(1,0) 
bands, among others, could produce predissociation of NO in the mesosphere. 
The dissociation coefficients due to these three bands are also included in 
the present work, assuming an average value for the absorption cross section 
of 0 2 at the corresponding wavelengths. 

For C0 2 , cross sections measured by Nakata et al, [1965] and Thompson 
et al. [1963], and for CH^, values reported by Okabe and Becker [1963] are 
adopted. In both of these, cross section profiles are important at wave- 
lengths below 1800 A only, and hence, their photodissociation rates are 
respectively small. 

The complete set of absorption cross sections used in the present work 
are shown in Figures 2.4 and 2.5. 
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2.2.3 Solar variation of 'photodissociation coefficients. The entire 
photochemistry of the upper atmosphere rests on the photodissociation of 
its constituents . The rates of production* as well as the loss of certain 
constituents are solely governed by their photodissociation processes. 

The diurnal b*,»,>vior of these constituents could, therefore, be expected 
to have a strong dependence on the solar variation of the respective photo- 
dissociation coefficients . 

In certain constituents such as 0 2 and 0 3 , different regions of the 
spectrum yield different species. These have different chemical proper- 
ties, and therefore, to evaluate their yields separately the photo- 
dissociation coefficients over these different spectral regions have 
been determined separately. The complete list of the photodissociative 
reactions is given in Table 2.2. The calculated coefficients for both 
X = 18° and 60° are illustrated in Figures 2.6 and 2.7. 

For heights where the optical depth factor is small, the photo- 
dissociation coefficients are independent of both the height and the 
solar zenith angle. Towards the end of the penetration depth, the dis- 
sociation curves tail off as the exponential term in equation (2.1) takes 
large values. It is only at these heights that the effect of solar zenith 
angle and the 0 2 column density become prominent on the photodissociation 
values. As the zenith angle increases from 18° to 60°, the tail in the 
dissociation curves can get lifted by about 4-5 km, while the values de- 
crease by more than an order of magnitude. However, the role played by 
the dissociation coefficients become less important as their values drop 
off. Thus, it is really a limited height interval over which the solar 
variation begins to show up in the chemistry of the constituents. 

Exceptions to this behavior are the C> 2 and H 2 0 dissociation coef- 
ficients. These two curves are composed of contributions made by different 
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Figure 2.7 Photodissociation coefficients of neutral constituents having low 
dissociation energies, for solar zenith angles 18° and 60°. 
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regions of the respective spectra. Bach component has a different pene- 
tration depth and when one takes the sum, the height over which the solar 
variation occurs is extended. 

It should be noted that any variation of the 0 2 column density can 
cause the dissociation coefficients to vary in a manner similar to that 
caused by the zenith angle variation. Hence, in the time-dependent cal- 
culations, the dissociation coefficient at time t was computed by taking 
0 2 column densities evaluated at the previous time step, (£-A£) . At 
zenith angles greater than 75°, the sec x term in equation (2.2) was 
replaced by the Chapman function, Ch{x ) . 

2.3 Reactions of Oxygen Constituents 

The problem of the production of atomic oxygen in the upper atmo- 
sphere has been the topic of several investigators [Nicolet and Mange » 
1954; Colegvove et al. 1965; Shimazakiy 1967], The production of the 
atomic oxygen species, as well as ozone, is a consequence of the photo- 
dissociation of 0 2> The energy absorbed in the Schumann- Runge continuum 
results in the formation of one atom in the ground state 0( 3 P), and the 
other in the metastable state Of^D), out of one molecule of oxygen. Wave- 
lengths above 1750 A, produce both atoms in the ground state. The shorter 

wavelengths do not penetrate to lower heights, and therefore, the direct 
production of the metastable species from C> 2 is restricted to higher alti- 
tudes, generally above 80 km. The metastable atoms thus formed undergo 
quenching reactions with N 2 and making them to attain gTound- state 

0( 1 D) + N* 2 0( 3 P) + N 2 

OfS) + o 2 -> 0( 3 P) + o 2 


energies, according to 
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The oxygen atoms recombine with each other and with C >2 in the presence 

1 

of a third body to give metastable 0^( A^) and 0^, respectively, according 
to 

0 + 0 + M ■+ ) + M, 

2 g' 

o + o 2 + m o 3 + m 


The metastable 0 2 (^A^] gets converted to ground level 0 2 in the quenching 
reactions. 


° 2 c\) ♦ N 2 + 0 2 * N 

VV + °2 * °2 * 0 

and also through the spontaneous emission 


2 * 

2 , 


OjC 1 * } + o 2 + h\>. 

Ozone, through photodissociation, yields both ground-state and excited- 
state molecular and atomic oxygen (Table 2.2]. As the photodissociation 
coefficients of 0^ extend to lower heights, these processes become the 



and the stratosphere. 

The xist of the reactions used in the present calculations is given 
in Table 2.3, along with the adopted rate constants. Many of these re- 
actions have been reviewed by SehLff [1969, 1972]. In reactions involving 
atomic oxygen, quantitative measurements are generally difficult to make 
because of the errors that creep in due to wall effects. Most of the 
measurements of the rate constant of Reaction 1 have been made at high 
temperatures where the temperature dependence is found to be 2 1 " 1 . However, 


48 


Table 2.3 
Oxygen reactions. 


Reaction 

No. 

Reaction 

Rate Constant 

n -l* 

cm sec 

R 1 

C •*- 0 + M 0 ( 1 A ) + M 

^ 9 

+ 3,0(-33) (T/300)" 2,9 
Campbell and Thzmah [1967] 

[ 2 

0+0 2 + M + °3 + M 

6.57 (-33) expfl. 01/RT) 
Huie et al. [1972] 

R 3 

0 * °3 * °2 * °2 cl V 

1 . 05(-ll) exp(-4.31/RT) 
McCrumb and Kaufman [1972] 

R 4 

o/t g ) * n 2 - o 2 *k 2 

4.0C-19) 
Zipf [1969] 

R 5 

° 2 C 1 V + °2 + °2 + °2 

2.22C-18) CT/300) 0,78 
Findlay and Snelling [1971] 

R 6 

0 2 ( 1 A^ r ) -*■ 0 2 + fcv 

2.6C-4) 

Badger et al. [1965] 

R 7 

oc^) + o 2 -* o + o 2 

5 . 5 (-11) 
Sdhiff [1972] 

R 8 

O^D) + N 2 0 + N 2 

5.0C-11) 
Sdhiff [1972] 

R 9 

0( 1 D) + 0- ->■ 0 + o^c 1 * ) 

2 2 K g J 

1.0C-12) 

Vallance- Jones and 
Gattinger [1963] 

r io 

oc 1 d) to 3 o 2 + o 2 

3.0C-10) 
Sdhiff [1969] 

II 

for binary reactions, and 6 for 

3-body reactions 


Read a(-n) as a x 10 





at temperatures below 350 K, the temperature dependence was found to be 
nearly T . At D-region heights, therefore, this reaction is relatively 
unimportant. In his latest review, Sahiff [1972] has recommended the 
value 1.2 (-35) exp(2.1/RT) for the rate constant of Reaction 2 (Table 2.3). 
In this work, however, a more recent measurement [Huie et al. 3 1972] has 
been adopted. This reaction is different from many other 3-body reactions 
in that it has a negative activation energy. It also happens to be the 
main source of 0^ in the upper atmosphere, and therefore plays an impor- 
tant role in aeronomical studies. 

*7 

The binary reaction between 0( P) and 0^ (Reaction 3; Table 2.3) 
begins to dominate as a loss process for 0^, when losses due to direct 
photodissociation are absent . The value for its rate constant is taken 
from the recent work of McCrumb and Kaufman [1972]. This value agrees well 
with the results of Krazenski et al. [1971] within experimental errors. 

It is, however, smaller by nearly a factor of 2 than the values deduced by 
previous workers [ Sahiff 3 1969]. The rest of the reactions given in 
Table 2.3 are all concerned with the excited species, 0 2 (*A ) and 0(*D) 
which are produced in the photodissociation of 0 7 and 0^, and also through 
Reactions 1 and 3. Their loss mechanism is mainly through quenching by 
N 2 and C> 2 , and the corresponding rates have been adopted from the reviews 
of Zipf [1969] and Sahiff [1972]. 

2.4 u&y gen- Hydrogen Reactions 

Bates and Nioolet [19S0] first introduced the oxygen-hydrogen reac- 
tions into the upper atmosphere photochemical models. In their work, the 
presence of various 0-H products in the mesosphere was hypothesised in 
order to explain the OH airglow emissions observed to originate in the 
mesosphere. However, many of the parameters that went into the calculation 
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of the distribution of the species were not known at that time. Hence 
their results were rather of a qualitative nature. A more detailed 
calculation of the behavior of the O-H products and their influence on the 
ozone distribution was carried out by Hunt [1966] . He also included non- 
equi librium conditions in his calculations. Heestvedt [1968 j extended 
these studies by introducing the effects of vertical eddy transport on 
the distribution of the minor neutral species at mesospheric and lower 
thermospheric heights. More recent studies of this subject which include 
the solution of a system of time- dependent continuity equations having 
transport terms were carried out by Skimazaki and Laird [1970] , Hunt 
[1971b] and Thomas and Bowman [1972] . 

The basic process responsible for the production of 0-H constituents 
in the upper atmosphere is the photodissociation and oxidation of water 
vapor which is assumed to be constantly transported into the stratosphere 
from the troposphere. The exact mechanism of this transport through 
the tropopause is not well understood. It is expected that the tropopause 

with temperatures below the freezing point would act as a barrier for any 
transport of water vapor through it. The water vapor that is carried to 
the stratosphere and above is subject to photochemical loss processes; 
its dissociation by solar radiation, and oxidation by 0( 1 D], according to 

H 2 0 + hv H + OH 

H 2 0 + 0( X D) OH + OH. 

In the stratosphere where the photodissociation coefficients are small, the 
main loss process is the oxidation reaction. 

The by-products of the H 2 0 dissociation and oxidation reactions, viz. 
OH and H are highly chemically active species. In the presence of ozone 
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and atomic oxygen, these constituents participate in catalytic reactions, j 

giving rise to other 0-H products such as H0 2 and H 2 0 2 , and also to molecu- 
lar hydrogen, The two reactions 

C ?) + OH -> H + 0 2 
H + 0- -*• OH + 0 2 

-11 3 -1 

both of which have high rate constants greater than 10 cm sec 

[Kaufman, 1969], however, do not yield any new constituents. They merely 

<? ,1 

interchange OH and F while converting 0('P) and 0^ into molecular oxygen. j 

3 I 

Nevertheless, they are important as loss processes for 0( P) and 0^. j 

The main sources of H0„ are the reactions 1 

^ ! 

0( 3 P) + OH + M -> HO + M, ;j 

^ i 

3 ■ i 

H + 0 3 -V H0 2 + OC P). | 

■ 1 

Many workers in the past have also included the reaction 

OH + 0 3 + H0 2 + 0 2 

as a possible source of H0 2> However, Sahiff [1972] has cautioned against 
using this reaction as no direct evidence for its occurrence with the ground- 
state OH is available. The H0 2 thus formed undergoes fast reactions with 
0( 3 P) and H destroying it to produce OH according to 

H0 2 + 0( 3 P) + OH + 0 2 

H0 2 + H OH + OH, 

respectively. Another reaction that has been postulated as converting i 

H0 2 back to OH is 

H0 2 + 0 3 + OH + 20 2 
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According to Sdhiff [1972], this again has to be dismissed as no direct 
evidence is available for its occurrence. It is seen now that the produc- 
tion and loss of both OH and HO^ are closely linked with the chemistry 
of 0( 3 P) and Q y 

Reactions between two H0 2 molecules result in the production of a 
single H 2 0 2 molecule, according to 

ho 2 + ho 2 - h 2 o 2 *o 2 

3 

while the loss of H 2 0 2 is governed by its reactions with Of P) and OH: 

H 2 0 2 + 0( 3 P] OH + H0 2 
H 2 0 2 + OH -*■ H 2 0 + H0 2 

While breaking up H 2 0 2 , these reactions reproduce OH, H0 2 and also H 2 0. 

The molecular hydrogen, H 2 is produced as a result of the reduction 
of OH, H0 2 and H 2 0 2 by atomic hydrogen, according to 

H + OH H 2 + 0, 

H + H0 2 H 2 + 0 2 , 

H + H 2 0 2 ■+ H 2 + H0 2 . 

Several reactions that take place between these newly formed 0-H 
constituents result in the reformation of H 2 0, the original source of all 
the 0-H constituents. Some of these reactions are 


S3 

OH + OH -> H 2 0 + 0, 

OH + H0 2 -► H 2 0 + 0 2> 

OH + H 2 0 2 -*■ H 2 0 + H0 2 . 

In addition to the above reactions, a large number of other reactions 
involving the hydrogen species take place at relatively slow rates. These 
reactions, however, can become important at certain heights and hours of 
the day when the concentrations of the constituents in the above-mentioned 
reactions drop to insignificant values. Such instances occur during the 
night time when both the 0( P) and H concentrations below about 60 km 
decrease by several orders of magnitude from the values occurring in the 
daytime. 

The complete list of reactions involving 0-H species, along with their 
rate constants is given in Table 2.4. Many of these reactions and the 
available data on the rate constants has been reviewed by Kaufman 
[1964, 1969], and Schofield [1967], For some of the fast reactions in- 
cluding H0 2 (Reactions 21, 23, 24, and 25] only the lower limits of the 
rate constants have been estimated in these reviews, and these have been 
used by all investigators in the past. For two of these reactions 
(Reactions 21 and 24]-, improved estimates have been made recently by 
Eochanadel et at. [1972], and these have been adopted in this study. For 
the other two reactions, lower limits given by Kaufman [1964] have been 
used. However, for the production of OH, as well as for the loss of H 
and H0 2 , the role played by these reactions is of minor importance. Hence 
any errors in these rate constants are not expected to influence the results 
of the present computations to any great extent. 


Reaction 

No. 

Reaction 

Rate Constant 

n -1* 

cm sec 

R 11 

H + 0 2 + M -*• H0 2 + M 

+ 3.0(-32) (273/T) 1 * 3 
Schofield [1967] 

R 12 

H + H + M + H 2 + M 

1.2 (-32) (273/T) 0,7 
Kaufman [1969] 

R 13 

H + 0 3 OH + 0 2 

2.6C-11) 
Kaufman [1969] 

R 14 

H 2 + 0 ->■ OH + H 

7.0C-11) exp(-10.2/RT) 
Kaufman [1969] 

R 15 

H 2 + 0( 1 D) -> OH + H 

7.S(-11) 
de Move [1967] 

R 16 

H 2 0 + 0(*D) OH +OH 

3.5(-10) 
Schiff [1972] 

R 17 

OH + 0 ■+ H + 0 2 

S.O(-ll) 
Kaufman [1969] 

R 18 

OH + 0 3 ■+ H0 2 + 0 2 

See Text 

R 19 

OH + OH -v H 2 0 + 0 

2. 7 (-12) 

Eochanadel et al . [1972] 

R 20 

OH + H + M -*• H 2 0 + M 

2. 5 (-31) 
Kaufinan [1964] 

R 21 

H0 2 + 0 -*■ OH + 0 2 

6.6(-ll) 

Eochanadel et al. [1972] 

R 22 

H0 2 + 0 3 -> OH + 2 0 2 

See Text 

R 23 

H0 2 + OH * H 2 0 + 0 2 

2.0(-10) 

Eochanadel et al. [1972] 

R 24 

H0 2 + H H 2 + 0 2 

> 3.0(-12) 
Schofield [1967] 

R 2S 

H0 2 + H •+• OH + OH 

l.O(-ll) 

Schofield [1967] 

R 26 

ho 2 + ho 2 + h 2 o 2 + o 2 

6.0(-12) 

Eochanadel et al. [1972] 
Paukert and Johnston [1972] 

R 27 

H 2 + OH + H 2 0 + H 

6. 3(-ll) exp(-5.49/RT) 
Kaufman [1969] 

R 28 

H 2 0 2 + 0 -»• OH + HO^ 

1.0(-15) 

Fonev and Hudson [1962] 

R 29 

H 2°2 +H - H 2 + H °2 

1.2C-11) exp(-4.2/RT) 
Kijeuski and Tvoe [1971] 

R 30 

H 2 0 2 + OH -V H 2 0 + H0 2 

4.0(-13) 

Fonev and Eudeon [1962] 
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The rate constant for Reaction 26 (Table 2.4) has been estimated 
recently by Eochanadel et al. [1972], The value quoted is the average of 
the two measured values. The reaction between OH and 0^ (Reaction 18) has 
been observed to take place only when the OH is in the vibrational ly ex- 
cited state, and the rate constants have been measured for the states 
v = 2 to v = 9 [Coltharp et al. 3 1971], This reaction becomes the major 

3 

source of H0 2 at nighttime when both 0( P) and H vanish at night below 
80 km. Therefore, even a small value assumed for its Tate constant can 
make a significant contribution towards the production of HO^. Though 
such a value can be assigned to this reaction by extrapolating the measured 
values down to v * 0, this reaction has been left out from the calcula- 
tions in view of the absence of direct evidence for its occurrence. For 
similar reasons. Reaction 22 (Table 2.4) has also been left out of the 
present calculations. 

2 . 5 Nitrogen Reactions 

The photochemistry of N-0 reactions has been pursued in the past in 
view of the importance of NO in the formation of the D region [Nicolet, 
1965; Norton and Barth 3 1970; Strobel et at. 3 1970; Strobel 3 1971a, 1971b]. 
Nicolet, and also Norton and Barth, in their calculations included only 
the photochemical continuity equations. Identifying the weakness of these 
photochemical models, Shimazaki and Laird [1970], and Strobel et al. 

[1970] introduced diffusive transport terms into the continuity equations. 
Schimazaki and Laird's calculations are restricted to the chemistry of 
the neutral species only, while the work of Strobel et al. and Strobel is 
restricted only to the thermosphere and the mesosphere, respectively. 
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The chemically inert nature of N 2 make the photochemistry of N-0 
reactions somewhat different to that of 0-H reactions. The absorption 
of ultraviolet radiation by N 2 is very little. Its atomic species, 
the ground-state N( 4 S) and the excited N( 2 D) are produced from N 2 through 
ionization by high energy EUV radiation and X-rays. Therefore, the photo- 
chemistry of these constituents should also include the ionic species, 

N, + , N0 + , and 0 * Most of the EUV radiation and X-rays do not penetrate 

d d 

below about 90 km, and hence the production of these atomic species and 
the subsequent formation of NO is restricted to altitudes above this level. 

The ionization rates used in the calculation of production of ionic 
species N 2 *, 0 + , 0^ , and N0 + ivill be discussed in detail in Chapter 4. 

The important reactions which produce the neutral atomic nitrogen species 
from the above ions are: 

0 + + N 2 + N0 + + N( 4 S), 

N 2 + * 0 N0 + + N( 4 S), NC 2 D), 

N0 + + e -*■ 0 + N( 4 S), N( 2 D], 

N 2 + + e - N( 4 S) + NC 4 S). 


In the reaction of N 2 with 0 and the recombination of NO , the 

2 

formation of the excited state N( D) is energetically possible. The 

2 

yield of N( D) from these two reactions has important consequences m 
the production of NO. The conversion of N( 4 S) to NO is effected through 


the reactions ; 




57 


N( 4 S] + 0 2 NO + 0, 

N( 4 S) + 0 3 + NO + 0 2 , 

N( 4 S) + °2^V N0 + °> 

N( 4 S) + 0 + M ■* NO + M. 

However, the largest contribution in the thermosphere comes from the 
reaction 

N( 2 D) + 0 2 + NO + 0, 

-12 3 -1 

which has a rate constant of 6 x 10 cm sec [Lin and Kaufman 3 1971], 

-17 3-1 4 

compared to 5 x 10 cm sec for the reaction between N( S) and 0 2 - 

2 

The problem here is that the exact yield of N( D) in the above reactions 
is still not known accurately. Strobel [1971b] showed that a percentage 

4 * 4 * 

branching ratio of 100 for the N 2 +0 reaction, and 75 for the NO re- 
combination reaction could give results in agreement with the observa- 
tions of Meira [1971] . Strobel in his work has also included the 
2 

production of N( D) through the impact of fast electrons with N 2 - How- 
ever, the contribution of this process below 120 km is not ' r ery impor- 
tant, and hence it was neglected in this study. The predissociation of 
N 2 has also been suggested as a possible source of atomic nitrogen in 
the lower thermosphere. According to Hudson and Canter [1969], however, 
its yield of atomic nitrogen is small compared to other sources, and 
therefore, this process has not been included here. 

While NO is produced from N in the thermosphere, reactions between 
the two also form major mutual loss processes for the two constituents: 
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NO + N ■+■ N 2 + 0 . 

Therefore, the abundance of NO in the lower thermosphere is strongly 
coupled to that of N. 

In the lower regions the oxidation of N 2 0 by 0(*D) can produce NO 
molecules as well as nitrogen molecules according to 

N 2 0 + O^D) + N 2 + 0 2J 

-> NO + NO. 


Nicolet and Peetermans [1972] in th; ir recent study of the stratospheric 
N^O have assumed equal branching ratios for this reaction. The photo- 
dissociation of N 2 0 } having a dissociation energy of 38.5 kcal, is 
capable of producing N 2 , 0( 1 D) J 0(*S), N( 2 D) and NO, according to 


N 2 0 + ho N 2 + 0(D) , 


+ N 2 + 0( S) 


A < 3400 A 
A < 2100 A 


■* NO + N( S), 


A < 2515 A 


NO + N( D). 


A < 1695 A 


However, photolysis studies of N 2 0 carried out at wavelengths 1849, 2139, 
and 2288 A show that the excited oxygen atom produced is exclusively the 
OC^D) state, and that the yield of NO and N makes only a negligible con- 
tribution towards the photolysis of N 2 0 [Preston and Barr 3 1971]. Since 
the major contributions to the photodissociation of N 2 0 come from wave- 
lengths larger than 2100 A, the N 2 0 photodissociation is assumed to 
produce only N 2 and 0( 1 D). 
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The main loss process for NO is due to its reactions with 0 and 0 3 , 
which convert NO into NO2 : 

0 + NO + M + N0 2 + H, 

NO + 0 3 N0 2 + 0 2 , 

NO + 0 •+• N0 2 + Tzv . 

However, NO is reformed through the reaction of N0 7 with atomic oxygen: 

N0 2 + 0 -v NO + 0 2 . 

The photodissociation of N0 2 , which has a coefficient as high as 
- 3 - 1 

5 x 10 sec -..during the daytime, also reproduces NO: 

\ 

\ 

N0 2 + h\> ■+• NO + 0 

The chemistry of NO and N0 2 are, therefore, strongly coupled to 
each other through 0 and 0 3 , somewhat similar to the situation with H0 2 
and OH. During the daytime, N0 2 is lost rapidly through these reactions 
and photochemical equilibrium is attained. For NO, however, the loss 
processes indicated above are not so rapid, having the lifetime against 
chemical losses of the order of a day in the mesosphere. Though the 
pre-dissociation of NO as suggested by Stvobel [1971b] is a major sink 
for NO in the mesosphere, it is not sufficiently fast enough to make the 
lifetime of NO short. It is therefore necessary to include diffusive 
transport in the continuity equations for NO. For the purpose of cal- 
culating the initial values in the next chapter, however, only the 
photochemistry has been used. 

The reactions between nitrogen species and hydrogen species have 
been thought as the final sinks for the odd nitrogen species produced in 
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the thermosphere and the mesosphere. It is believed that species such 

as HNC >2 and HNO^ produced in the stratosphere diffuse downwards and 

finally get washed down with rain water [ Sahiff a 1972]. The presence 

of HN0 3 above the tropopause has already been detected by Murcray et al, 

[1969] and Harries [1973] using spectroscopic techniques, and found 

_ 3 

mixing ratios to be about 3 x 10 ppm. Crutzen [1971] who studied 
the production of 0^ in a 0-H-N atmosphere found that the data available 
at that time for the relevant reaction rates could not explain the 
observed concentrations of HNOj. 

The production of HN0 2 and HNO^ takes place according to the 3-body 
reactions 

OH + NO + M HN0 2 + M, 

OH + N0 ? + M ■+ HN0 3 + M 

The HN0 3 formed gets destroyed through the oxidation reaction 

HN0 3 + 0( 3 P) OH + N0 3 , 

and by photodissociation, both of which occur only during the daytime 
in the stratosphere. 

The reaction between H and N0 2 is a possible means of converting 
N0 2 into NO: 

H + N0 2 -*■ NO + OH. 

However, at heights where N0 2 concentrations are large, H concentrations 

3 

are much smaller than 0( P] concentrations, so that it cannot compete 

3 

with the reaction between N0 2 and 0( P) in the conversion of N0 2 and NO. 
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Strobel [1972b] has also included the reaction 

N( 4 S) + OH + NO + H 

in his study of the 2?-region nitric oxide abundance. As a means of 
converting N( 4 S) into NO, this reaction has to compete with the reaction 

N( 4 S) + o 3 + NO + o 2 , 

-13 3 -1 

which has a rate constant of about 7 x 10 cm sec . Since the con- 
centration of OH is much less than that of 0^ at all heights, the reac- 
tion between N( 4 S) and OH will be of no importance, unless its rate 
constant is several orders higher than that of the reaction with 0^. 

This reaction is therefore not included in the present study. 

The complete list of the neutral reactions between 0-N species and 
O-N-H species considered in this study are listed in Table 2.5. Some 
of these 0-N reactions and their rate constants have been discussed by 
Schiff [1969] in his review of reactions involving oxygen and nitrogen. 
Unlike the case of the 0-H species, many of these reactions are relatively 
slow, and the only fast reaction here involving ground-state species is 
the Reaction 39. The reaction between N and N0 2 may also be considered 
fast. This reaction produces N 2 and 0, and to a lesser extent NO 
(Reactions 42 and 43) . 

The more recent data adopted here are the rate constants of Reaction 37 
(Table 2.5) measured by Slanger and Black [1970], Reaction 34 measured 
by Lin and Kaufman [1971], Reaction 41 measured by Davis et al. [1973], 
and the oxidation reaction of N^O (Reactions 44 and 45) measured by 
Saott et al, [1971]. Among the N-0-H reactions, the rate constants of 
Reactions 46 and 47 are those measured recently by Westenberg and 


Table 2.5 

0-N and O-N-H reactions 


Reaction 

Reaction 

Rate Constant 

No. 


1* 

cm see" 


N + 0 2 

-t- NO + 

0 

N + 0 2 C 

V 

-y 

NO + I 

n + o 3 


NO + 

°2 

n( 2 d) + 

°2 

-*■ 

NO + 0 

N( 2 D) + 

0 

N 

+ 0 

NO + 0 

«* 

no 2 

+ hv 

NO + 0 ■ 

► M 


no 2 + i 

NO + 0 3 


N0 2 

+ °2 

NO + N 


n 2 + 

0 

NO + NO 

+ c 

*2 

1 N0 2 ' 

no 2 + 0 


NO 

+ ° 2 

N0 2 + N 

■+ 

n 2 o 

i + 0 

N0 2 + N 


NO 

+ NO 

n 2 o + 0 

C X D) 


n 2 + ' 


N 2 0 + 0(D) -»• NO + NO 


OH + NO + H HN0 2 +■ M 


OH + N0 2 + M -*■ HNOj + M 


HN0 3 + 0 -»■ OH + N0 3 


H + N0 2 -»• NO + OH 


OH + N('S) •+ NO + H 


6.5(-12) exp(-7.0/RT) 

Sahiff [1969] 

2.8C-1S] 

Clark and Wayne [1969] 

3.4C-11] exp(-2.4/RT) 

Phillipg and Sahiff [1962] 

6.0C-12) 

Din and Kaufman [1971] 

2.0C-13) 

Weill [1969] 

6.4C-17) 

Fontigin et al. [1964] 
6.8C-32) 

Slanger and Black [1970] 

9. SC-13) exp(-2.46/RT) 

Sahiff [1969] 

2.2C-11) 

"Phillips and Schiff [1962] 

1.0(-33) 

Sahiff [1969] 

9;12(-12) 

Davis et al. [1973] 

7. 3(-12) 

Phillips and Sahiff [1965] 

6.0 C— 12} 

Phillips and Schiff [1965] 
l.O(-lO) 

Scott et al. [1971] 

1.0C-10) 

Scott et al. [1971] 

9.2(-31)CT/300)' 4-0 
Westenberg and de Haas [1972] 

1.6(-30) (T/300)" 2,8 
Weetenberg and de Haas [1972] 

£ 2.0C-14) 

Morns and Niki [1971] 

1.2 (-9) exp (- 1 . 93/RT) 
Schofield [1967] 

See Text 


* 

n = 3 for binary reactions and 6 for 3-body reactions 
+ Read a(-n) as a x 10" n 
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de Haas [1972] . The value for the HNO^ oxidation reaction is only an 
upper limit. These values differ widely from those used by Crutzen 
[1971] in his calculations of HNO^ concentrations. 

2 . 6 Carbon Reactions 

The chemistry of carbon species has not drawn much attention in the 
past as compared to the chemistry of other species, in photochemical 
model calculations of the terrestrial atmosphere. The importance of 
including the C-O-H chemistry in the studies of the other minor neutral 
constituents was highlighted by Strobel [1972a] who suggested that CO 
could act as a sink for OH above 80 km. Further, CH^ is believed to 
yield H.>0 through oxidation processes increasing the mixing ratio of H^O 
in the upper stratosphere [Nicolet, 1970]. Both of these constituents, 
OH and ^0 play important roles in the ion chemistry of the D region. 
Hence it is important to investigate the effects of the chemistry of 
C-O-H species on the distribution of OH and H 2 0 in the mesosphere. 

The sources of C-O-H species in the upper atmosphere are the photo- 
dissociation and oxidation of C0 2 and CH^. Hays and Olivero [1970] 
studied the production of CO due to photodissociation of C0 2 only, while 
Wofsy et at. [1972] investigated the production of CO due to photo- 
dissociation and oxidation of CH^ in the troposphere and stratosphere, 
and its effect on the distribution of 0-H species. 

C0 2 is relatively inert and it does not pay any significant role 
in the chemistry of neutral constituents. However, it gets photo- 
dissociated by solar radiation below about 2000 A yielding CO and 
atomic oxygen. These products recombine to form C0 2 in a 3-body 

CO + 0 + M C0 2 + M 


reaction 
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While this reaction converts CO into CO 2 efficiently at higher altitudes, 
the reaction 


CO + OH ■* C0 2 + H 

becomes important in this oxidation process at lower altitudes. However, 
as a sink for OH, this reaction becomes more important at higher alti- 
tudes where the concentration of CO is expected to be large. According 
to Hays and Olivevo [1970], the CO profile obtained by considering only 
the photodissociation of C0 2 has a peak of about 30 ppm between 100 and 
120 km, the exact height depending on the eddy diffusion values and the 
rate constant of the recombination reaction used in the calculations. 

It is known that CH^ is present in the lower stratosphere in con- 
centrations slightly above 1 ppm around 24 km [Bairibridge and Heidt 3 
1966; Kyle et al. 3 1969]. In the stratosphere, CH^ is subject mainly 
to oxidation losses. Reactions with both 0 and OH result in the removal 
of one H atom from CH^ yielding CH^. The removed H atom forms an OH 
molecule with 0, and a H 2 0 molecule with OH, according to 

CH 4 + 0 CH 3 + OH, 

CH 4 + OH •* CH 3 + H 2 0. 

The CH 3 formed gets converted to CH 2 0 by its reactions with 0 and 

° 2 ' 

CH 3 + 0( 3 P) -»■ CH 2 0 + H 

CH 3 + 0 2 CH 2 0 + OH. 

% 

Here, the by-products are H and OH, respectively. Next, 0("P1 and OH 


I 


i 


oxidize CH 2 0 to produce HCO through the reactions 

CH 2 0 + 0( 3 P) -*■ HCO + OH 

CII 2 0 + OH •* HCO + H 2 0, | 

1 

' . } 

yielding OH and H,jO. In the final step, HCO is converted to CO j 

j 

according to j 

3 i 

HCO + 0( P) ^ CO + OH ; 

HCO + 0 2 CO + HO r ] 

! 

! i 

• i 

It is therefore seen that the four H atoms in the CH^ are removed 
by 0 and OH, forming in the process OH and H 2 0, respectively. 

Table 2.6 shows the complete set of reactions used in the calcula- 
tions. The rate constant for Reaction 51 ivas in dispute for several 

•••; j 

years, mainly in view of the uncertainty in its activation energy. The 

latest values for this rate constant were obtained by Slangev et at, j 

[1972], who found that it has a positive temperature dependence, with j 

CO as the third body. In this work, the value obtained at 296 K with 
N 2 as the third body was employed along with the temperature dependence 
obtained with CO. 

Established rate constants for the reactions with CO and CH^ 

[Reactions 52-55] have been adopted here. For the rest of the C-H-0 
species, most of the rate constants are not so well known. The value 
for Reaction 56 is only a lower limit, while that of Reaction 57 is an 
approximate figure. The values for Reactions 59 and 62 are lower 
limits. The rate constant for Reaction 61 is only an estimate. 





Table 2.6 
C-O-H reactions 



n * 3 for binary reactions and 6 for 3-body reactions 
Read a(-n) as a x 10~ n 






In view of the uncertainties in the rate constants of these reac- 
tions, the height profiles of these intermediate species cannot be con- 
sidered as accurate. However, the loss rates of CH^ as well as the 
production and loss rates of CO can be evaluated with less uncertainty 
as more reliable values for the rate constants of these processes 


become available. 
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3. MODEL CALCULATION OF NEUTRAL SPECIES DISTRIBUTION 

Most of the minor neutral constituents discussed in the last chapter 
are produced as a result of the photodissociation and oxidation processes of 
the major constituents 0 2> N 2 , C0 2 and those transported from the tropo- 
sphere to higher levels through various transport mechanisms. In this 
chapter, a time- dependent transport model will be developed to investigate 
the distribution and the diurnal behavior o ii these minor neutral consti- 
tuents in the height interval 30-120 km. 

It is necessary first to set up a model for the transport processes 
which are responsible for carrying the constituents such as H 2 0, N 2 0 and 
CH^ from the troposphere into the stratosphere and above. This is discussed 
in the following section. The method of solution of the continuity equa- 
tions which are in the form of second- order partial differential equations 
is discussed next. It is then possible to calculate the set of initial 
values required as input parameters for the solution of the time-dependent 
continuity equations. These equations are next solved for the case of a 
constant solar zenith angle which corresponds to noon conditions. From 
the results of these calculations it is possible to eliminate those con- 
stituents that need to be further investigated for their diurnal behavior. 
Finally the solar zenith angle variation is incorporated to simulate the 
diurnal behavior. These calculations are continued over several simulated 
days until 24-hour reproducibility is reached. 

3.1 Transport Processes in the Upper Atmosphere 

Several types of motions are responsible for the transport of minor 
constituents in the upper atmosphere. These include macroscopic motions 
such as mean motions, wave motions, and winds as well as transport mechanisms 


I 
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such as turbulence and molecular diffusion. Of these motions, only those 
which cause vertical transport of minor constituents are considered in the 
present study. 

Mean motions are associated with the large-scale general circulation 
of the atmosphere. The vertical components of these circulations con- 
stitute the mean-vertical motions. Thermodynamic considerations show that 
mean- downward motion results in the heating of the atmosphere while the 
upward motions cause cooling. The constraints on the atmosphere heat 
budget show that such mean-vertical motions cannot have velocities more 
than a few centimeters per second. According to Gudiksen et at. [1968], 
however, the mean-meridional circulation is found to play only a minor 
role in the transport processes at stratospheric heights. At higher alti- 
tudes its effect is considered to be further reduced. Hence, this motion 
has not been included in this study. Zonal and other prevailing winds, 
being associated with horizontal motions, are also excluded. 

Transport effects caused by such irregular motions as wave motions, 
tides, and turbulence are all described under a single mechanism, viz. 
eddy diffusion, in terms of a lumped parameter called eddy diffusion co- 
efficient. This quantity is analogous to the molecular diffusion coeffi- 
cient, except that the former deals with the diffusion of 'eddies', while 
the latter deals with the diffusion of molecules. 

3.1.1 Eddy diffusion coefficients. Transport by eddy diffusion is 
considered to be taking place through the exchange of 'eddies' having 
scale lengths ranging from several meters to several kilometers. These 
eddies are formed by such phenomena as wind shears, internal gravity waves, 
tides, and in the case of large-scale eddies, by planetary waves. 


[Sheppard^ 1959; Hines 3 1963; Hodges 3 1969; Geislev and Dickinson, 1968]. 
Tile main outcome of such exchange processes is local mixing. 

The behavior of air parcels under eddy diffusion has been evident in 
experiments conducted to study the dynamics of the upper atmosphere. In 
these experiments chemi luminous vapor trails are released from rockets 
and their rate of growth is observed visually or photographed from the 
ground. Such observations have shown that the trails initially grow under 
molecular diffusion and later disperse into erratic shapes under eddy dif- 
fusion below a certain level. Above this level, the growth is controlled 
only by molecular diffusion at all times [Blamont and de Jag:,v 3 1961; 
Zirmevman and Champion 3 1963; Justus 3 1969], This transition level where 
the turbulence ceases allowing molecular diffusion to dominate is generally 
known as the turbopause. Blamont and de Jager [1961] found that this 
turbopause level appeared distinctly on each occasion the observations 
were made, with its value lying between 102 and 110 km. According to 
Zirmevman and Champion [1963] the transverse growth of the turbulent clouds 

observed in the height range 60-100 km showed two sequential time depen- 

1/2 

dencies; an initial growth due to molecular diffusion with v * t and a 
subsequent growth due to eddy diffusion with v a t, where v is the Gaussian 
half-width of the clouds and t is the time Irpse. The diameter of these 
clouds generally indicates the order of magnitude of the scale length of 
the eddies. 

Lettau [1951] expressed the eddy diffusion coefficient, as 
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where X^ is a length fixing the scale of turbulent displacements (mixing 

length), and is a velocity fixing the speed of turbulent motions. 

Using this definition, Zinur* -wan and Champion [1963] calculated the value 

6 7 2 -I 

of D which was found to be between 2 x 10 and 1.5 x 10 cm sec" at 

3 

7 8 

75 km and between 1.5 x 10 and 10 at 100 km. More recent studies 

6 

by Keneshea and Zimmerman [1970] show D & having values between 2 x 10 
6 2—1 

and 8 x 10 cm sec" around 100 km. 

The observation of the vapor trails from ground level below can detect 
only the horizontal spread of the turbulent clouds, and as such the values 
obtained are in respect of the horizontal component of the eddy diffusion 
coefficient. In the lower thermosphere where the negative lapse rate of 
the temperature tends to inhibit any vertical movements, one cannot expect 
the turbulence to be isotropic. However, Keneshea and Zimmerman [1970], 
considering the scale lengths of turbulent sources near 100 km, assumed 
, vertical component of the eddy diffusion coefficient to be equal to 
the measured horizontal component. 

In another method, it is possible to estimate by studying the heat 

budget of the lower thermosphere and the mesosphere. Here, the diffusion 

coefficient for thermal conductivity of the atmosphere is calculated, and 

from the equality between this parameter and the eddy diffusion coefficient, 

the latter is estimated [Johnson and Gottlieb 3 1971; Olivero 3 1970]. 

This method really gives upper limits for the eddy diffusion coefficients, 

6 2 -1 

and shows that the D values decrease from about 3 x 10 cm sec near 

e 

5 2 -1 

100 km to about 10 cm sec near 60 km. 

In the stratosphere, the negative lapse rate of the temperature makes 
the atmosphere stable against vertical convections. This in turn, lowers 


i 
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the value of vertical eddy diffusion coefficient. The effective values 
for this quantity in the stratosphere have been determined by studying the 
distribution of the radioactive debris from upper atmosphere nuclear ex- 
plosions. Karol [1966] has deduced from aircraft measurements of radon 
daughters RaD and RaF, values in the range [0, 2-2.0) x 10^ cm 2 * 4 sec - * for 
the vertical eddy diffusion coefficient in the lower stratosphere between 
35°N and 34°N latitudes. In a different method, Reed and German [1965] 

attempted to explain the large-scale mixing processes using published heat 

3 3 2 

flux data, and arrived at values lying between 1 x 10 and 5 x 10 cm 
sec -1 for the vertical component of the lower stratosphere eddy diffusion 
coefficient for the months of July and September. 

Using a different approach, Skimazaki [1971] determined effective 
eddy diffusion coefficients which yield solutions to model calculations 
that would agree with the observed composition structure in the lower 

7 

thermosphere. He found that a D & profile having a peak value of 1 x 10 

2 -1 

cm sec in the height range 95-105 km, and a peak half-width of about 
8 km explains most of the rocket-obersved composition height variations. 
However, in view of the uncertainties in other variables such as solar 
flux, reaction rate constants, etc., it is difficult to arrive at unam- 
biguous results from such model calculations. 

From the results of previous workers [ Keneshea and Zimmermans 1970; 
Shimazaki and Laird, 1970] it is known that values of eddy diffusion coef- 
ficient play a crucial role in determining the height profiles of minor 
neutral species. Hence, in this work, two models for height profiles 

have been adopted, with the values fixed on the basis of the above dis- 

4 2 -1 

cussions. In both cases, is kept constant at 10 cm sec below 


.1 
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40 km. Between 40 and SO km, its values are given by (Figure 3.1). 


D = io 
e 


(3.2) 


where a and b are constants, and /As = (s-80) km. The values of a and b 

6 2 — 1 

are selected such that at 80 km, D = 2 x 10 cm sec in the 'high' 

5 2 -1 

model, and D = 4 x 10 cm sec in the 'low* model. Above 80 km, D is 

<3 0 

described by the following expression so that it peaks at 105 km. 


D _ yf + MAs) 2 + e(As) 3 


C3.3) 


where a 3 b s and a are constants, and As = (s-105) km. The 'high* model 

7 2 -1 

has a peak value of 1 x 10 cm sec and the 'low' model a peak value 
6 2 -1 

of 5 x 10 cm sec , both at 105 km. The above expression shows to 
drop rapidly near 110 km, so that molecular diffusion takes over above 
these heights. 

3.1.2 Molecular diffusion coefficients. From the kinetic theory of 
gases, the molecular diffusion coefficient, I> m can be derived as 

* ■ w iM ' 2 ■ 

where |Tj is the total number density, d is the mean molecular diameter, m* 
is the reduced mass of the molecules and k is the Boltzmann’s constant. In 
the upper atmosphere, the molecular diffusion coefficients of constituent 
gases is calculated for the motion of the gas molecules in the medium of 
an N 2 ~ atmosphere. Hence, for a constituent other than oxygen, [X] is 
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essentially equal to [N 2 ]. This quantity decreases exponentially in the 
atmosphere , and up to about 100 km, the variation of T does not exceed a 
factor of 1.5. Therefore, in this height interval, D » increases exponen- 
tially with the same scale height as that of the density variation. About 
100 km, the increase is further enhanced due to the increase in T which 
takes place in the thermosphere until it reaches a constant value. 

The values of D ^ for atmospheric gases at standard temperature and 

pressure lies between 0.14 and 0.21 cm sec , except for hydrogen and 

2 -1 

atomic oxygen whose values are 0.67 and 0.47 cm sec , respectively, 
[LettaUj 1951]. The variation of D m for oxygen with altitude is illus- 
trated in Figure 3.1. It is seen that is much less than Z? g below about 
110 km. Above 110 km, molecular diffusion dominates over eddy diffusion. 

Direct observations of the spread of rocket released vapor trails have 
established the existence of a transition level separating the molecular 
diffusion and eddy diffusion regions, as mentioned before. 

The observation of the spread of vapor trails above the turbopause 
yields the molecular diffusion coefficients. Values obtained by Golomb 
and MaoLeod [1966] and other workers are in general agreement with the 
values calculated from equation (3.4). 

3. 2 Solution of Continuity Equation with Transport Term 

Diffusion results in the loss of a constituent within a given volume. 
The rate of loss of its concentration is given by the divergence of its 


flux caused by diffusion. Since the transport effects are considered 
only in the vertical direction, this can be simply written as 


am 

U 


dif 


9 $ 

9a 


(3.5) 
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where M is the concentration of the constituent X, and <j> is the 
vertical flux. 

The continuity equation with this transport term then appears as 

^ - 9 - fffl - H (3.6) 

where Q is the production rate per unit volume and P is the loss co- 
efficient. The loss due to recombination of two like particles is 
neglected here. 

It is seen from equation (3.6] that the diffusion term is important 
only if the chemical loss rate is smaller than the loss rate due to dif- 
fusion. These rates are described in terms of respective time constants, 
which are given by (neglecting molecular diffusion) 

r(chem) = 1/P (3.7) 

x(dif) = Jp /D , (3.8) 

CL & 

where is the scale height of the mixed atmosphere. With the V fj 
profiles used, the time constant against eddy diffusion in the meso- 
sphere is of the order of a day. Therefore, for constituents having 
fast chemical loss rates, it is not necessary to include the transport 
term in their continuity equation. For such cases the continuity equa- 
tion is simply 


= Q - P M 


(3.9) 
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For constituents having chemical time constants more than about a 
day, the transport term is included. The total flux due to diffusion 
is given by 


d> = d> * A 
v 




(3. 10) 


where d> and d> axe the flux contributions due to molecular diffusion 
rrj e 

and eddy diffusion, respectively. The molecular diffusion flux is 
expressed in the form [ Chapman and CcwVlng 3 1952] 


m ' 


m 


/iia 

\ as 


+■ (1+Qt) 


m 

T 


I? + ml 

33 hJ 


C3.ll) 


where a is the thermal diffusion factor and E \ is the scale height of 
the constituent X.. In the same manner, Cotegvove et at. [1965] ex- 
pressed the eddy diffusion flux as 


<f> = 

y e 


- D 


/m + m . M 


e [ as 


t as ff f 

a 1 


(3.12) 


where E is the scale height of the mixed atmosphere. 


Combining equations (3.11) and (3.12) one gets 


+ = - {^(s)^||l + S( 3 ) [J]j , (3.13) 


4(s) = D + D (3.14) 


where 
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B(s) = { d & + (l+a)'0 m 


U 

T 


3 T jn 

3s + ff + H. 

a t- 


C3.15) 


Substituting equation (3.13) in the continuity equation for the £ 
constituent (equation (3.6)), 


th 




* 2 [X\i nx] i 


where 


a M 

b = si 


+ fl 


C 



(3.16) 


(3.17) 


(3.18) 


If the distributions of all the constituents are known at time 
t-At, it is possible to evaluate the loss term P and the* production 
term Q, and in turn, the parameters a and d appearing in equation (3.16) 
written for time t. Strictly, the values of F and Q deduced at time t 
should be used to evaluate these parameters. ‘This, however, is not 
possible, unless one does an iterative calculation for each time step. 
Nevertheless, if the variation of the distributions between successive 
time steps is small, the use of P and Q corresponding to time step 
(£-A£) is not expected to introduce any serious errors in the solution 
for [X\, at time t. From trial calculations, this assumption was found 
to be valid, except at sunrise. Since the emphasis of this work is to 
study the variations occurring at daytime, this approximation was used 
and any deviations occurring at sunrise were hot considered important. 


i - J 1 - J -.1 
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The concentrations of the constituents at height s and time t 
are therefore given by a system of non-linear equations, comprising 
a set of second-order partial-differential equations (equation (3-, 26)) 
and a set of single-order ordinary differential equations (equation 
(3.9)). In view of the above mentioned assumption it is possible to 
solve these equations taken individually. 

In order to solve these non-linear equations, one has to resort to 
numerical techniques. It is customary in these methods to replace the 
differentials with difference quantities; and to solve the resulting 
algebraic equations. . In converting the differential equation to a dif- 
ference equation, the space between the two boundaries z ^ and z^ is 
divided into N- intervals with spacing A z, so that 


z = Zq * nhz 3 n - 0, 1, 2, ...A? 


(3.19) 


Similarly, the time lapsed from an initial time t = 0 is expressed as 


t = mAt, m = 0, 1, 2, ... (3.20) 


The concentration of a given constituent whose distribution is governed 
by equation (3.16) can then be written as (see Appendix II for details 
of converting the differential equation to its difference form) 


a y 
n 3 


n+1 


jn~\ m m 

+ 6 y + y y , 
n 3 n ' n J n~L 


3 


6 m-l 

n 


(3.21) 


and the coefticient s a *. 0 , and y^ are functions o£ a« b+ As* and ht* 

n n ti 

3 m ~ , in addition, includes terms containing [X] .(3.,t-Ai) j ^ i 3 
n j 

w X 

while 6 ~ ± contains these terms as well as [XL* (z 3 t-At) . For sim- 

Tt *lr 

plicity omitting the prefix m and suffix n from the coefficients in 
equation (3,21), one gets 


ay,, + 3y + Y y , 

°n+ 1 a n a n~ 1 


For N intervals, one can write CX-2) such equations containing N 
unknowns. The two boundary conditions provide the two additional equa 
tions required for the solution of this set of simultaneous equations. 
The values of [Xl^ at t = 0 also have to be specified as initial 
conditions. 

This system of equations can be expressed in the form 


w m = Pi 


where [Aj is a tridiagonal matrix of order (ff-2) and [B] is a column 
matrix. The most direct method of solving this system of equations is 
to invert [A] and express [T] in the form 


p 


[B] 


The finite-difference scheme is only an approximate method and it 
is expected that the solution obtained by this method is also an approxi- 
mation to the exact solution of the differential equation. If the solu- 
tion is stable, it Will converge to the exact solution as A z •*- o. 

However, in the computer solution of these equations, the accumulation 
of round-off errors results in the solution diverging more from the exact 
solution as the number of steps is increased by making As smaller. This 
problem was overcome by adopting a more efficient method than the above 
one, developed mainly for the solution of a system of equations whose 
coefficients form a tridiagonal matrix [Richtmyer, 1957]. The details 
of this method are found in Appendix II. 

The boundary conditions required for the solution of the above set 
of equations are obtained from known physical quantities at the two 
boundaries. These can be either the number density the flux across 
the boundary. With the number density, one specifies either y q or j/ ff , 
while with the flux as a boundary condition, one specifies the gradient 
of y n at the boundary. 

In order to get the necessary information regarding the number densi- 
ties or the flux at the boundaries, one has to depend on direct measurements 
of these quantities or estimates of them by other indirect means. In the 
next section, sources of such information and the values adopted as 
boundary conditions are discussed. 

The transformation of the ordinary differential equation (equation 
(3.9]) into its difference form yields 




if 


1 + P(s, At 


(3. 25) 
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The use of this implicit form in the conversion from differential to 
difference equation results in better stability [Bhimasaki, 1967] , and 
does not impose restrictions on the choice of the height and time 
increments , 

From trial calculations it was found that height increments of 1 km 
were most convenient with regard to resolution and computer time. Time 
increments during the daytime were selected corresponding to x = 10° U P 
to x = 70° and thereafter x = 5° up to x = 100°. During nighttime one- 
hour time increments were used. 

3.3 Boundary Conditions 

The boundary conditions required for the solution of the continuity 
equations having the transport terms appear either as a number density 
or a flux across the boundary. The constituents for Which transport 
effects are important can be divided into tWo groups. Those originating 
in the troposphere and transported Upwards into the stratosphere, such 
as H 2 0, N 2 0, C0 2 , and CH^ fall into one group. The constituents in the 
other group, such as 0, NO, CO, H and H 2 are produced in the lower 
thermosphere and upper mesosphere and are transported either downwards 
or upwards. The boundary conditions that are applicable for these two 
groups are therefore quite distinct. 

3.3.1 Boundary conditions for lower atmosphere produced consti- 
tuents. In the stratosphere and mesosphere, these constituents are 
subject to losses due to either oxidation or photodissociation. Hence 
an upward flux across the lower boundary is required to maintain their 
concentrations at these high altitudes. The most suitable boundary con- 
dition for them is to specify this upward flux. Values of this quantity. 
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however, cannot he measured directly and can only be estimated from measure- 
ments of their number densities at the lower boundary. Such measurements 
have been carried out in the lower stratosphere using either air-borne or 
balloon-borne instruments, and these have been employed in specifying 
the lower boundary conditions for the above constituents. 

Though H 2 0 is one of the important minor constituents in the upper 
atmosphere its concentration above the tropopause in only little known. 

Early measurements carried out using frost-point hygrometers have shown 
wide disparity in the mixing ratios obtained above the tropopause, with 
the values varying from a few parts per million to about 40-50 ppm 
[Gutnick 3 1961]. More recent measurements carried out using improved 
techniques, however, have shown consistently lower mixing ratios, in the 
range 2-3 ppm between 100 and 40 mb levels or 16. and 22 km [Masteribrookj 
1968, 1971]. In view of the precautions taken in these measurements to 
minimize errors due to any contaminations and also, because of the con- 
sistency in the measurements which have been taken over a period of six 
years, the lower values are preferred to the high values reported pre- 
viously. Water vapor concentrations deduced from solar spectra recordings 
using air-borne instalments have also shown mixing ratios in the range 
2.4 - 2.6 ppm at 18 km in the mid- latitudes [McKinnon and Morei)ood 3 
1970]. The measurement camried out near the stratopause by Sdhotz et at. 
[1970], however has shown a high mixing ratio, about 6 ppm. This increase 
in the water vapor concentration above the lower stratosphere values has 
given rise to much speculation regarding production of water vapor through 
other means such as oxidation of methane, as mentioned before. For the 
purpose of fixing the boundary value, however, a mixing ratio of 3 ppm 
at 30 km has been adopted. From trial calculations it was found that an 
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upward flux of 1 x 10 cm sec with an eddy- diffusion coefficient of 
4 2 -1 

10 cm sec at 30 km gives rise to this mixing ratio. 

For nitrous oxide, which is the most abundant nitrogen compound in 
the lower atmosphere next to N 2 , direct measurements are available only up 
to about 24 km. These measurements carried out by Sohutz et at. [1970] 
show mixing ratios about 0.1 ppm near 24 km. Some recent spectral obser- 
vations have shown that the mixing ratio reduces from a tropospheric 
value of 0,27 ppm to a value 0.21 ppm taken 3 km above the tropopause 
[Harries 3 1973] . Theoretical calculations of Bates and Hays [1967] 

show that the mixing ratio at 30 km could lie in the range 0.02-0.15 

3 4 

ppm for stratospheric eddy diffusion coefficients between 10 and 10 
2 -1 

cm sec , corresponding to a tropospheric mixing ratio of 0.25 ppm. 

On the basis of these data, a mixing ratio of 0.1 ppm has been adopted 

7 -2 -1 

at 30 km, which corresponds to an upward flux of 1 x 10 cm sec at 
this altitude. 

In the case of C0 2 , the loss rate due to photodissociation is rather 

small, and hence it maintains a constant mixing ratio up to the mesopause 

level [Hays and Olivero 3 1970] . The mixing ratio in the troposphere has 

been found to be in the range 310-320 ppm [Seiler and Jmge 9 1970] and 

a similar value has been detected near the stratopause [Scholz et at. 9 

-4 

1970], To be compatible with a mixing ratio of 3.14 x 10 ppm, a flux 

_2 -1 

value of 6.4 x 10 cm sec has been assumed as the boundary condition 
at 30 km. 

Methane is a trace constituent of high abundance in the troposphere. 
According to measurements of Bainbridge and Heidt [1966], it has a 
constant mixing ratio of 1.6 ppm in the troposphere while around 24 km 
its mixing ratio drops to about 1.0 - 1.3 ppm. Kyle et at. [1969] found 


85 


that the total content of methane above 30 km is about 1/5 that above 
the tropopause. McConnell et al. [1971] calculated the vertical dis- 
tribution of CH^ below 20 km using eddy diffusion coefficients determined 

from trace studies, and arrived at a mixing ratio of 0.9 ppm and a flux 
9 -2 -1 

of 1 x 10 cm sec at 20 km. Their results show that the flux drops 

by about an order of magnitude between 10 and 20 km. Assuming the same 

8 “2 

loss rate holds between 20 and 30 km, the same flux value, 1 x 10 cm 
sec"* was adopted as the lower boundary value for CH^, and from trial 
calculations, this was found to yield a mixing ratio of 0.15 ppm at 30 km. 

In specifying the upper boundary conditions for these constituents 
H 2 0, N 2 0, C0 2 , and CH^, as well as for 0 2 , one can make use of the fact 
that all of these constituents are subject to photodissociation and are 
lost above the upper boundary. Upward flux across the upper boundary re- 
plenishes this loss. Near this altitude level the dissociation rates are 
independent of the height as their optical depths are nearly zero. In 
the case of 0 2 , however, the level of zero optical depth is near 160 km. 
Though its dissociation rate has a small variation between 120 and 160 km, 
for the purpose of calculating the upward flux its value has been assumed 
to be constant. Assuming each of these constituents has a distribution 
under its own scale height abovj 120 km, the flux can be determined from 

W = w* -w w C3 - 26) 

where J .(zj\ and 27 .( 3 ,.) are the dissociation rates and the scale heights 
z'- 2lr ir 

of the constituents evaluated at 120 km. 

3.3.2 Eonndczry conditions for upper atmosphere produced 
constituents . The constituents in the second group all move downwards 
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after being produced in the thermosphere, except in the case of H and 
Hj. In the case of 0 and CO, the downward flux across the upper boundary 
is determined by the rate at which 0 2 and C0 2 are carried upwards. The 
time average of the downward flux of 0 and CO has to be equal to that of 
the upward flux of 0^ and C0 2 , respectively, in order to conserve the 
total amounts of these species. Strictly speaking, this equality may 
not hold if one takes the instantaneous values of the fluxes. However, 
for the sake of simplicity the boundary values for these two species have 
been calculated with this assumption, and are given by 


" 2<f> o 2 ts i!p * 

(3.27) 

WV = “ *C0 o Cs ¥ 3 ■ 

(3.28) 


The estimation of the downward flux of NO at 120 km is more complex, 
and one has to go into the details of ion chemistry in the E and F 
regions to evaluate this quantity. Stvobel [1971a] has studied this 

8 

problem and found that the NO concentration at 120 km is about 1 x 10 
cm and that it has no diurnal variation below this height. This con- 
centration is also in agreement with the measurements of Meira [1971]. 

Since no estimate of the flux value at 120 km is given, a number density 
8-3 

of 1 x 10 cm - is used as the upper boundary value for NO. 

Atomic hydrogen, being light in weight escapes from the earth's 
atmosphere, and gives rise to the geocorona surrounding the earth. From 
measurements of the intensity of the geocorona radiation, and also from 
measurements of the absorption of the solar Ly-a line by geocoronal 
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atomic hydrogen the abundance of this species; in the thermosphere and 
above can be deduced. In these calculations it is customary to use a 
distribution profile such as the Ko charts -Nico let model which takes into 
account the upward escaping flux, and normalize the concentration at the 
base of the thermosphere to fit the observations. The atomic hydrogen 

• 7 

concentrations so determined near 100 km fall in the range (2.5 - ) x 10 
-3 

cm {Meier and Prim j 1970; Meier and Mange 3 1970; Vidal-Moddan et al. 3 

7 -3 

1973]. Based on these values, a number density of 1 x 10 cm was 
adopted as the upper boundary condition for H. 

Molecular hydrogen is produced mostly in the mesosphere through 
various chemical reactions, and it is possible to have both an upward 
flux above the mesosphere as well as a downward flux below the mesosphere. 
However, at the upper boundary level no evidence is available to estimate 
flux. Hence, it has been assumed that this flux is zero at the upper 
boundary. Since H ^ is rather inactive in the ion chemistry of the D 
region, any errors in the H ^ profile caused by improper upper boundary 
value is not of much significance. 

At the lower boundary, the above constituents which have their 
source in the upper altitudes are subject to chemical reactions having 
short time constants. A zero flux at the lower boundary is therefore a 
reasonably valid assumption for use in this study. 

3.4 Calculation of the Initial Values 

In the solution of equations (3.22) and (3.25) for the number density 
of the constituents at time t 3 a knowledge of these values at time {t-ht) 
was required. Therefore, these initial values have to be calculated first, 
before one can solve the above time dependent equations. This was done by 
employing a steady-state semi -transport model. 
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First, the distribution of the pure oxygen constituents, 0( P) , 

0[ 1 D), and 0^ are obtained using the Chapman's reaction scheme. Next, the 
height distribution of the constituents belonging to the first group men- 
tioned in the last section, vis. H 2 0, N 2 0, C0 2> and CH^ are calculated 
using steady-state continuity equations which included transport terms 
and loss terms due to photodissociation and oxidation subject to the 
boundary conditions mentioned in the last section. Finally, the height 
distributions of the rest of the constituents are obtained assuming photo- 
chemical equilibrium corresponding to noon conditions. 

By writing the continuity equations and assuming photochemical 
equilibrium, the following expressions for the oxygen species are readily 
obtained: 


3 t J 1 J 3 f0 2 ] 1 1/2 

[<K p,] = w/l * R 2 R 3 [o 2 n m! 


(3.29) 


EO^D)] 


CR 2 [0C 3 P)] [Ml + J la ) [0 2 ] 
R g [N 2 ] + (R 7 +R 9 ) [0 2 j 


(3.30) 


R [0 2 ] M [0( 3 P)J 

Ptf - - LJ -r 


(3.31) 


J 1 =J a* J b 


J Z = J 3a + b + J Z o 


where 
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The production of OC P) vanishes at nighttime, but its loss processes 
are maintained throughout the day by collisional reactions. Hence, at 
lower heights where the collisional rates are large due to the high density 
of the atmosphere, its concentration at night will vanish. At higher alti- 
tudes, where the collisional losses are small, the loss rates are con- 
trolled by diffusive processes having lifetimes of a few days. Therefore, 

3 

the 0( P) concentrations are maintained at these heights even at nighttime. 
The excited-state, 0(*D), abundance, on the other hand, vanishes at all 
heights during the night. In the case of 0^, its production is proportional 
to the square of the neutral particle density while its loss rate, mainly 
due to photodissociation, remains almost constant at all heights during the 
day and vanishes during the night. Therefore, 0^ is expected to have 
higher values at nighttime than at daytime. 

The photochemical loss term used in calculating the 0('P) concentra- 
tions as given in equation {'5.12'), yields unrealistically high values 
above about 90 km. It is known from the work of previous investigators 
[Shimazaki and Laird , 1970; Keneshea and Zimmerman , 1970] that the 0C P) 
profile has a peak around 90 km with the distribution above this peak 

3 

given approximately by the diffusive equilibrium conditions. Hence OC P) 
values above 95 km were adjusted to have an exponential drop with the 
scale height factor determined by its own molecular weight. 

In calculating the height distribution of the 0-H constituents, it 
Is convenient to start with an profile that has already been trans- 
ported into the high altitudes subject to photodissociation and oxidation. 

In this preliminary calculation, the recycling of H^O through reactions 
between constituents such as HC^ and OH was neglected. Also steady- 
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state conditions were employed. Therefore, the continuity equation for 
HO was written independent of the concentrations of other species, and 
the profile thus obtained depends only on the upward flux and the eddy 
diffusion coefficient used. 

The photodissociation products of 1^0, vis. OH, H0 2 , and H 9 0 2 have 
short time constants against chemical losses and under steady conditions 
their production rates become equal to their photochemical loss rates. 

This is not so in the case of H and H 2 . However, to obtain the initial 
distributions, photochemical equilibrium was assumed for these two species 
as well. The production and loss terms for the odd hydrogen constituents, 
H, OH, and H0 2 are written using the reaction scheme given in Table 2.4. 
These three continuity equations when added together yield the expression 


R 24 [h][ho 2 3 + r 23 [ho 2 ][oh] = J- 7 [H 2 o.] 


(3.32) 


Next, taking the major terms in the equations for OH and H0 2J 


[OH] - 


r 21 [ho 2 ][o] + R 13 [h][o 3 ] 

R^Io] 


(3.33) 


[ho 2 ] = 


R 11 [H][a 2 ][»] 

k 21 [o] 


Substituting these equations in equation (3.32) yield 

_ f J 7 [ H 2 0] , 1/2 

1 1 ' VcR 24 ♦ S *2j)/ 


(3.34) 
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( 3 . 35 ) 




a- « m/ * n m> 

ft <’ w n H/m ♦ * 

It Ut tlWt pOQtilhto to iMiitufa the of feet a of the Reactions 11 and 
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(3.36) 


Itllit 
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(3.37) 


h\ nvhuH en»\\Krt«H«n n qxmdratie equation, obtained after elimi- 
jx^| in aquation ($*$<0 using equation ($.57) was solved. The 
iinni valuer wow xxbtaixxeti aftox' iterate the whole set of equations 
VM-W wart x^aehed, Since the effects of fH) terms on [Of 3 ?)] 
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VoxiiOnvO whu x^aehod X'&fddly « The x-esults of these calculations are 
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necessary to calculate the distribution of the charged species which 

4 2 

appear in the continuity equation for N( S] and N( D). For this purpose, 
an electron-density profile given by 

[e] = 10 C S “ 5 0)/ 10 , s < ioo km 


and 


(3. 38) 


[e] = 10 5 100 <s < 120 km 

liras assumed. The region of the electron-density profile which is impor- 
tant for the production of N( D) , the principal source of NO, lies above 
the ledge where the cluster ions are absent. Hence, the distribution of 
the positive ions could be expressed in a straight forward manner in terms 
of the assumed electron densities. The calculation of the ion-pair pro- 
duction functions required for this purpose will be given in detail in 
the next chapter along with the relevant reaction rates. 

2 

The quasi-equilibrium photochemical continuity equations for N, N( D) , 
NO, and N0 + can be combined to yield the equation 


[N] [NO] + <2 [NO] = b ■' 

where 

a = CR 36 [0] + R 37 [0][Afj + R 38 [0 3 ])/2 R 39 

b = CR 45 [N 2 ol [OC^D] + R 101 [0 + ][N 2 ] + R 103 [N 2 + ][0] + N 2 [N 2 + ][ e ]]/R 39 


Eliminating [N] in the above equation using the expression for [N] obtained 
from the photochemical continuity equation one gets a quadratic equation 
in [NO], the solution of which gives the required initial values. 
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The [N0 2 ] and [Ih D)] profiles are next calculated using the photo- 
chemical equations, 


[W 2 ] - 


((R 36 *R J7 [«mo] * R 38 [0 3 ]) [NO] 
-r 4 * R 4i[°'l 


(3.40) 


[N( 2 D)] = 


Tl R 10 3f N 2 + lf°] * ^2 

R 34(°2] + R 3S C°I 


(3.41] 


where y^ and y 2 are the branching ratios of the reactions N 2 +0 and 
N0^ + e, which yield NC^D). 

The resulting profiles of the 0-N species are illustrate! in Figure 

3.4. 

The initial distributions of C0 2 and CH^ are first determined 
employing the steady-state transport equations having loss terms due to 
photodissociation and/or oxidation, subject to the boundary conditions 
given in Section 3.3. The behavior of CO is somewhat similar to that of 
0( P) above the mesopause. Its chemical loss term has a long time constant 
in this region and, therefore, the profiles calculated using photochemical 
equations alone would yield unrealistically high values. Hence, the 
height distribution of this constituent was also obtained by solving the 
steady-state transport equation. The profiles of [CH^] and [CO] thus 
obtained are shown in Figure 3.5. 

In addition to CO, the CH^ and CO^ dissociation and oxidation yield 
minor compounds such as 0^0, CH^, and HCO..- Of these only CH 2 0 reaches 
any significant level, and this is shorn in Figure 3.5 along with the 
profiles of [CH 4 ] and [CO]. 
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The calculated amounts of these two constituents, CH^ and CO were 
found to have negligible effect on the distribution of the 0-H species, 

OH, in particular [Ratnasiri and Sechnist* 1971] . Hence, these con- 
stituents are not incorporated in the diurnal model calculations. 

3.5 Constant Zenith Angle Calculations 

The height distributions of the minor neutral constituents obtained 
as initial values in the last section do not represent any real situation. 
The calculation of more realistic values involves the simulation of the 
diurnal variation of the incident solar flux, and repeating the calcula- 
tions over several diurnal cycles. In this study, however, an inter- 
mediate set of calculations has also been made in order to investigate 
the effects of transport mechanisms on the distribution of the long-lived 
constituents. This is done by keeping the solar zenith angle constant at 
18° corresponding to noon conditions. Only the high values have been 
used in these calculations which are continued with time increments of 
one-half hour until convergence in the distribution of most of the con- 
stituents is reached. 

3.5.1 Oxygen and oxygen-hydrogen constituents. The height dis- 

3 1 1 

tributions of the oxygen species, 0( P} , 0( D) , 0 2 ( A ) and 0 3 calculated 

C 

using the transport model with the solar zenith angle fixed at noon, are 
shown in Figure 3,6. Comparing these with the corresponding distributions 
obtained as initial values using the photochemical model, it is seen that 
significant differences between the two sets of results appear mostly 
above about 55 km. 

With the eddy diffusion coefficients adopted, the time constants 
against eddy diffusion in the stratosphere are found to be of the order 
of several months. As such, for any constituent having a small chemical 
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loss rate, the effects of diffusion do not show up in a time-dependent 
calculation where the simulation is carried out over a few days only. Any 
changes in the height distributions in this region, therefore, occur due 
to changes in the photochemical loss and production rates, rather than 
due to diffusive processes. In both sets of calculations discussed above, 
the photochemical production and loss rates in the stratosphere are more 
or less the same. This explains the similarity in the two sets of results 
which appear below 55 km. 

On the other hand, near the mesopause the time constants against 
eddy diffusion are of the order of a day. Hence eddy diffusion controls 
the distribution of the species whose chemical loss rates have long time 
constants. This becomes apparent in the case of Of P) , where the two sets 
of results show marked differences in the mesosphere and above. The 
[0( P)] profile reaches a peak between 90 and 95 km. Above this peak, the 
distribution is determined by eddy diffusion and the upward flux across 
the boundary. This upward flux causes the concentration in this altitude 
range to deplete below the values given by the photochemical model. Below 
the peak, both the photochemical and diffusive models show a sharp drop, 
the drop b- h ,; g more prominent in the latter than in the former. This is a 
consequence of the fact that the diffusive model is also time dependent 
while the photochemical values are obtained on the assumption of quasi- 
equilibrium conditions. When the oxygen species OC P) and 0 3 are allowed 
to remain with the hydrogen species H and OH, the resulting reactions cause 
the net conversion 

0C 3 P} + 0 3 — *■ 2 0 2 


Unless the eddy diffusion values are very high, the above conversion 

3 

proceeds with time, resulting in the sharp drop of OCP) concentration. 
This drop, however, is arrested below about 75 km when the H concentration 
drops causing one of the conversion reactions, 

H + 0 3 OH + 0 2 

to slow down. Further, in the initial set of calculations only this 

hydrogen reaction was included, and the rest of the reactions involving 

3 

both OH and Ho 2 were left out in obtaining the Ot P) and 0 3 values. 

The above reaction between H and 0 3 is also important as a loss 

process for 0 3 in the mesosphere. During the daytime, ozone is lost 

_2 

mainly through photodissociation which has a high rate of about 10 
sec - '*' down to the stratosphere. However, over a small height interval 
corresponding to the region where the peak in the [H] profile exists, 
the loss through this reaction becomes comparable to or even dominate 
the loss through photodissociation depending on the peak [H] value. 

This causes a cyclic chain of reactions to take place among the consti- 

3 3 

tuents 0( P), 0 3 , H and OH. The depletion of odd oxygen species 0( P) , 
and 0 3 caused by these chain reactions is responsible for the irregulari- 
ties present in both [OC^P)] and [0 V ] profiles around this height lange. 
When the eddy diffusion coefficients used are too low, this mechanism 
can even produce a sharp valley in the [Of P)] profile. 

In the stratosphere, the production of odd oxygen species is con- 
trolled mainly by the Herzberg band dissociation of 0 ^. Because of the 
high efficiency of the three-body reaction 
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0 + 0 2 + M — *■ 0 j + M 

in the stratosphere, the oxygen atoms will get converted to 0 3 molecules. 
Between 0 and 0 3 , therefore, the dominant species is 0 3 at these altitudes. 

3 

The time constants for the chemical loss of both 0( P) and 0 3 are very 
short, and the vertical diffusive processes have no significant effects on 
their concentrations . As such, both the photochemical and the diffusive 
models give more or less the same values for the [0( P)] and [Og] profiles 
in the stratosphere. 

1 1 

The distribution of the two metastable species 0( D) and 0 2 ( A ) are 

^ y 

closely related to the distributions of the other oxygen species. The 
changes in the profiles of these two species can therefore be attributed 
to the changes talcing place in the [0( P) ] and [0 3 ] profiles. 

In Figure 3.7, the height distributions of the hydrogen species, H, 
H 2 , H^O, OH, H0 2 , and H 2 C> 2 are shown. The initial values of these consti- 
tuents were determined first by calculating the distribution of H 2 0 using 
a steady-state continuity equation having the transport term, and then 
using the photochemical equations to determine the distribution of the rest 
of the constituents. In these calculations, the reformation of H 2 0 through 
reactions occurring among other species was not considered. This, however, 
was accounted for in the new set of calculations. However, the two I^O 
profiles do not show any significant changes except that in the diffusive 
model, the values are slightly higher towards the upper end. 

Among the other 0-H species, the constituents that were affected most 
as a result of diffusive processes are H and H 2 » Both of these have long 
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time constants against chemical losses, particularly in the lower thermo- 
sphere. Atomic hydrogen, in a manner similar to the atomic oxygen profile, 
exhibits a peak abundance near 85 km. Above this peak, the profile re- 
sembles one given by diffusive equilibrium. Uelow the peak level, the 
photochemical and the diffusive profiles are generally of the same shape, 
though the absolute values differ by about a factor of 10 or more in the 
stratosphere. In view of the approximations made in calculating the 
initial values such a difference could be expected. 

For the molecular hydrogen too, values are mostly affected in the 
lower thermosphere where the diffusive processes cause an increase in 
its concentration. Around 70 km, the profile has a broad peak, and 
below this the values are only slightly changed in the new profile. 

The remaining 0-H species, OH, HC^, and H,^ all have short time 
constants at all heights during the daytime. Hence their concentrations 
could be determined by photochemical equations. Their chemistry is 
closely linked with that of 0( P) and H. As mentioned before, both of 
these constituents have their concentrations reduced in the lower thermo- 
sphere due to the action of diffusive processes. As a consequence, the con- 
centration of the species OH, HO 2 , and ^Og also drop at a faster rate 
above 80 km in the diffusive model. Below this height the three profiles 
have almost the same features as those found in the photochemical model. 

The differences in the absolute values that are present in these profiles 

3 

could be attributed to changes in the [0( P)3 and [0^] profiles. 

3.5.2 Nitrogen and oxygen-nitrogen constituents. This series of 
calculations was done first by determining the height profile of ^0. In 
both the initial and the new set of calculations, continuity equations 
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with the diffusion term were used. Hence the profile of N 2 0 appearing 
in Figure 3.8 is similar to that shown in Figure 3.4. There is no signifi- 
cant production and its abundance there is due to the particles carried 

upwards from the troposphere. 

4 2 

The atomic species N( S) and N( D) are calculated in a manner similar 

to that discussed in the previous section. The branching ratios used for 

2 

the reactions producing NC D) were both 0.75. It is seen that while the 
2 

profiles of N( D) in the two sets of calculations are similar, the 

4 

profiles of N( S) differ widely, particularly in the mesosphere. This 
difference can easily be attributed to the difference in the profiles of 
NO. This constituent has a time constant of the order of several dsys 
in the mesosphere and therefore the photochemical equations, as used in 
the initial set of calculations, do not yield correct values. Conse- 
quently, the [NO] profile shows increased values with a peak near 75 km. 

In the transport model, however, diffusion causes the profile peak to 
move downwards as the NO, which is produced in the thermosphere, is 
allowed to diffuse downwards. This causes a minimum near 85 km in the 

new NO profile, and above this height the profile remains almost constant 

7 -3 

at a concentration about 3 x 10 cm . Its broad maximum now lies around 

9 -3 

40 km with a value nearly 2 x 10 cm . The reduced NO concentration in 

4 

the mesosphere causes the corresponding N( S} concentration to increase, 

4 

as the loss of N( S) is mainly due to its reaction with NO. 

The changes occurring in the [NO] profile as a result of introducing 
the diffusive transport term into its continuity equation are also reflected 
in the distribution of N0 2 . This constituent and NO are strongly coupled 
to each other, as discussed in the previous chapter. During daytime, 

N0 2 is lost mainly through its reaction with 0( P) down to about 50 km, 
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and through photodissociation below 50 km. The time constants vary from 

about one second in the lower thermosphere to about a few minutes in the 

stratosphere. Hence, in this new set of calculations too, photochemical 

equations were made use of for the determination of the N0 2 profile. The 

production of NC^, however, depends on the concentration of NO. For this 

reason, the values of N0 2 were also reduced several orders of magnitude 

below the levels obtained previously. Near 75 km, this value is about 

1 x 10' cm v , while above this height, it remains constant corresponding 

to the constant values of NO. In the stratospheric height range 40-55 km, 

8 -3 

the N0 2 concentration is of the order of 10 cm . The [N0 2 ]/[NQ] 
ratio given by these profiles is about 0.5 below 35 km, and drops through 
several orders of magnitude in the mesosphere. In the photochemical 
model, however, N0 2 has a value of 1 x 10 cm - at 75 km. a factor of 
1000 more than the corresponding transport model value. For values 
towards the two boundaries, such as at 40 km and 100 km, the two models 
give similar results. 

The above results show that eddy diffusion affects not only the con- 
stituents having long time constants, but also the constituents having 
short life times as well. This is a consequence of the dependence of 
the short lived constituents on the long lived types, whose height 
distributions are determined by the transport processes. 

3.6 Diurnal Variation of Neutral Species 

The simulation of the diurnal variation of the constituent concen- 
trations was next carried out using the results obtained in the previous 
section as initial values. The zenith angle was varied in steps of 10° 
between 20° and 70°, and in steps of 5° between 70° and 100°. Between 
noon and sunset, and between sunrise and noon, the time increments required 
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for the solution of the partial differential equations were obtained 
from the above x increments. Beyond 100° of the zenith angle, time steps 
of one hour were used. Since the scope of this work is to study the 
changes of the V region taking place during the daytime, the detailed 
variations occurring at sunrise and sunset were not included. In all the 
calculations, the height increment was maintained at 1 km. 

During daytime the photodissociation rates were determined at each 
time step using the 0 2 concentrations required for the calculation of 
the optical depth factor obtained from the previous time step. For x 
angles greater than 100°, the photodissociation rates were all made equal 
to zero. At the upper boundary, constituents such as 0 2 , H 2 0, and N 2 0 
were subject to a flux boundary condition expressed in terms of their 
dissociation rates. Since the dissociation rates at night were zero, 
the fluxes of these constituents across the upper boundary during night- 
time were automatically made to vanish. In reality this may not be true. 
However, the errors introduced by using this method are not expected to 
affect the values below about 100 km. 

The calculation of the diurnal variation was done for two eddy 
diffusion models, discussed in Section 3.1.1. The simulation was con- 
tinued until 24-hour reproducibility in values were obtained for most of 
the constituents. For certain constituents such as NO and H this condi- 
tion was not reached within a reasonable period. Since this variation was 
small, the computation was, however, terminated after 10 days of simulation. 
On the final day of simulation, the concentration values of all the neutral 
constituents were written onto a magnetic tape at every time step be- 
tween sunrise and sunset, for use with the ionic species calculations. 
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3.6.1 Otcygen constituents. The noon distributions of the oxygen 

species 0,, 0(^P), 0(*D), and CLC^A ) obtained on the 10th day of simu- 
o i g 

lation are shown in Figures 3.9 and 3.10, which correspond to the high 

3 

and low eddy-diffusion models, respectively. Comparing the 0( P) profile 

in Figure 3.9 with that in Figure 3.6, which was given for a fixed sun, 

it appears that the two profiles differ only above 70 km. The main 

feature in the new profile is the reduced gradient between 70 and 85 km 

11 -3 

and the reduced peak value, which is now about 1 x 10 cm . In view 

3 . 

of the shorter duration over which photodissociation production of 0( P) is 

3 

allowed to take place, this reduction in the amount of 0( P) in the region 

where diffusion occurs can be understood. With the low eddy diffusion 

3 

model (Figure 3,10), however, the rate of removal of the 0( P) into the 

thermosphere is reduced and consequently the peak concentrations occurring 

11 -3 

near 90 km increases to a new value of 3.2 x 10 cm . When the transport 
effects becomes reduced, the coupling between the reactions involving 
0( P), OH, 0 3 , and H become stronger, and as a result, irregularities in 
the profiles of these constituents can occur in the regions where such 
coupling dominates. This is shown in the low eddy diffusion profile 

3 

of 0( P) which has a minimum and a maximum near 80 km. 

3 

The diurnal variation of Of P) under both high and low eddy diffusion 
values are shown in Figures 3.11 and 3.12, respectively. The high eddy 
profile shows a significant variation between noon and predawn in the 
lower thermosphere, whereas the corresponding variation in the low eddy 
profile is not very significant. The behavior of the two profiles in the 
mesosphere at x = 60°, both in the morning and evening, is similar 
though the absolute values differ slightly. In both 6f these sets of 
curves, the higher values in the thermosphere are given by the afternoon 
profile while tile morning profile shows lower values. 
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Figure 3.9 The height distribution of oxygen constituents at noon 3 calculated using 
the diurnal model with high eddy diffusion coefficient profile. 
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Figure 3.10 The height distribution of oxygen constituents at noon, calculated using 
the diurnal model with low eddy diffusion coefficient profile. 
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An interesting feature in these profiles is the increase in the 
morning values over the noon and evening values between 65 and 75 km. 

3 

This could probably be due to the higher contribution to 0( P) produc- 
tion due to 0 3 dissociation which has higher values over this height range 
in the morning than during noon or evening. The secondary peak appearing 
around 60 km in both the noon profiles is also absent in the two x 3 60° 
profiles. At this height, the difference in the concentration between 
the noon and the evening x ® 60° concentration is about a factor of 4. 

In the stratosphere, all three profiles have the same slopes though the 
actual values are lower by about a factor of 2 in the afternoon than in 
the noon. The predawn profiles, in both sets of profiles, vanish below 
about 80 km. 

The two ozone profiles shown in Figure 3.9 and 3.10 corresponding to 
high and low eddy diffusion models have more or less the same values below 
about 75 km. The low diffusion profile shows a small irregularity around 
85 km as mentioned before, so that around this height the profile corre- 
sponding to the low eddy diffusion has higher values. Between 85 and 75 km, 
however, it has lower values than the profile corresponding to the high 
eddy diffusion model. This behavior of 0^ is directly governed by that of 
0( 3 P), whose profiles are influenced greatly by the transport effects. 

The variations occurring in the 0 3 profiles are also reflected in 
the profiles of O^A ) shown in Figures 3.13 and 3.14, for the two cases 
of high and low eddy diffusion coefficients, respectively. Below about 
65 km, the photodissociation of C> 3 is mostly responsible for the produc- 
tion of CL (A I, and therefore its concentration follows closely that 
of 0 3 . Consequently, both the high diffusion profile (Figure 3.13) and 
the low diffusion profile (Figure 3.14) show similar values below 65 km. 
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In these profiles, the evening values exceed the morning values between 
55 and 65 km, while immediately below and above, the opposite takes place. 
The exact height limits depend on the eddy diffusion coefficient used in 
each case. The morning profile with the low D » also shows an irregu- 
larity at 80 km not shown in the other profiles, and this could be directly 
related to the behavior of the C> 3 profile in the morning. Between 80 and 
95 km, both sets of profiles do not show any variation with the zenith 
angle. Above this height again, however, there is marked diurnal varia- 

tion taking place, caused by the dependence of 0-( A ) production on the 

z 9 

3 

0( P) concentrations. 

The diurnal variation of the oxygen constituents has been studied 
recently by Shimazaki and Laird [1970, 1972], Thomas and Bauman [1972] 
and Hunt [1971b j 1973]. Direct comparison of the results of the present 
study with the results of these workers is rather difficult in view of the 
different parameters used by them. . Nevertheless, the general features 
found in the results of these studies are present in the profiles of the 
present study. 

Both the solar flux in the 1300 - 2000 A band and the eddy diffusion 
coefficient used in this study agree closest with those used by Thomas 
and Bowman and by Hunt. Hence, a meaningful comparison of the results 
obtained here could be made with their results. 

3 

The Of P) peak concentration at noon at 90 km obtained by Thomas and 

11 -3 

Bowman is about 1.6 x 10 cm while the corresponding value obtained by 
11 -3 

Hunt is 3 x 10 cm . These agree well with the results obtained with 

11 -3 

the low diffusion model which has a peak value of 2.5 x 10 cm at 

10 -3 

91 km. The minimum present near 70 km which has a value of 1 x 10 cm 
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also agrees well with the results of these workers. The nighttime values 
in all these models vanish below about 80 km, while above this height the 
profiles do not show any significant diurnal variation, as found in the 
present study. The slight discrepancies present in these height distri- 
butions are probably due to the different values adopted for various 
parameters such as the eddy diffusion coefficients, recombination rate 
constants and solar flux. 

3 

In recent times experimental determination of the 0( P) distribu- 
tion in the lower thermosphere has been carried out using rocket-borne 

instruments. Henderson [1971] obtained values between 1 x 10^ and 
11 -3 

8 x 10 cm in the height range 89-94 km. The presence of the peak 

in the 0( P) profile near 100 km was demonstrated by Sehols and Offerman 

11 11 -3 

[1974], who obtained peak values of 2.5 x 10 and 6 x 10 cm in two 

flights using mass spectrometric measurements incorporating cryoion 

sources. Good and Goloirib [1973], using N0/0 chemiluminescence measure- 

11 11 -3 

ments obtained peak values lying in the range 4 x 10 and 8 x 10 cm 
under different solar conditions, with the average peak height near 99 
km. It appears that generally, both the theoretical and experimental 
peak values lie in the same range, while the experimental peak heights 
are a little higher than the corresponding theoretical values. In view 
of the present uncertainties in the various parameters mentioned earlier 
the agreement between the model calculations and the measured values of 

3 

0( P) can be considered as satisfactory. 

The ozone profiles obtained in these model calculations all show 
the irregularity around 75-85 1cm, which is obtained with the low dif- 
fusion model in the present study. In Hunt's profile this has a value 
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around 1 x 10 8 era -3 , while in Thomas and Bowman's profile the flat portion 

7 -3 7 -3 

has a value near 3 x 10 cm . rhe value of 3 x 10 cm obtained in the 

present study agrees with this latter value. These results are however 

slightly less than the values deduced by Evans and Llewellyn [1970] who 

obtained peak concentrations of (1-3) x 10 cm in the height range 

85-90 km. The noon concentrations obtained near 60 km in all their 

10 -3 

profiles have more or less the same value of 1 x 10 cm , in agreement 
with the experimental values [ Hilsenrath , 1971]. 

The observations by Evans a>id Llewellyn [1970] of the 1.27 y 

1 ] 
emissions gave 0 7 ( A ) distributions having peak values of 2.5 x 10 
2 9 

cm -3 between 50 and 60 km during the daytime. These values are about a 

factor of 3 higher than the calculated values. Such a discrepancy 

could easily arise in view of the many uncertainties in the production 

rate of CUf^A ), such as the reaction rate constants and solar flux 
2 g 

intensities responsible for the photolysis of 0 3 . The calculated values, 
however, are in closer agreement with the profiles obtained by Thomas 
and Bowman [1972] . 

3.6.2 Oxygen-hydrogen constituents . The noon profiles of the 
hydrogen species, obtained after diurnal simulation are shown in Figures 
3.15 and 3.16, corresponding to the high and low values of D . respec- 
tively. The two sets of curves are almost identical below about 60 km 
while above this height the effects of using different eddy diffusion 
values become evident. 

The [I^O] profile, under high eddy diffusion gives values higher 
by nearly an order of magnitude in the thermosphere than under low eddy 
diffusion. The loss of H 2 0 in the mesosphere and above due to photo- 
dissociation is replenished by the upward flux generated in the troposphere. 
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Figure 3.15 The height distributions of hydrogen and oxygen-hydrogen constituents at 
noon, calculated using the diurnal model with high eddy diffusion 
coefficient profile. 
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Figure 3.16 The height distributions of hydrogen and oxygen-hydrogen constituents 
at noon, calculated using the diurnal model with low eddy diffusion 
coefficient profile. 
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As such, the higher the eddy diffusion coefficient, the greater the con- 
centration of H 2 0 at high altitudes, subject to the limiting values given 
by complete mixing ratios. Because of the large time constants involved 
in the loss of H 2 0, it has no significant diurnal variation. 

The low D atomic-hydrogen profile, while showing reduced values 
between 70 and 80 km, has increased values above this height range. At 
heights below, the values remain unchanged. The increase of values above 
the peak for this low profile is somewhat similar to the behavior of 
atomic oxygen in this height range. The values at 120 km are the same 
for both curves, as constrained by the boundary condition. The profile of 
molecular hydrogen corresponding to the low model also shows increased 
values above 65 km. 

The more active constituents, OH and H0 2 exhibit only a little change, 
with the variation of the eddy diffusion coefficients, and even this was 
evident above 70 km only. Since both OH and H0 2 have short time constants 
in this part of the atmosphere, their concentrations are not directly 
affected by eddy diffusion. However, the variations in both [0] and [0 3 ] , 
or [H], easily produce significant changes in [OH] and [H0 2 ] . Both sets 
of curves show near constant values below 70 km, while a sudden drop in 
their concentrations become a characteristic feature above this height. 

The [OH] profile shows an irregularity near 75 km, which is caused by the 
corresponding variations in both [0] and [H] at 75 km. This irregularity 
causes [OH] to have higher values than [H0 2 ] over the height range 73 and 
78 lan in the low diffusion model. In the high diffusion model there is 
no such cross-over. Further, the latter set of curves have slightly 
higher gradients than the low diffusion curves above 80 km. The H 2 0 2 
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profile for both low and high diffusion profiles give the same values 

except over the region around 83 km where the low £> values are slightly 

6 

higher than the high D g values. 

The diurnal behavior of H is shown in Figures 3.17 and 3.18 for the 

3 

high and low eddy models, respectively. As in the case of 0( P) , the [H] 
profile reaches a peak at 85 km and vanishes in the nighttime below 
this height. The evening x = 60° curve and the noon curve do not differ 
very much, except that the values cross over near 70 km. In both sets of 
curves near the peak, the evening values are slightly higher than the 
noon values, while below the peak, the noon values dominate. The morning 
values in both cases, however, exhibit a large swing between 60 and 80 km. 
The enhancements of H between 60 and 70 km can be regarded as a conse- 
quence of the photodissociation processes. The photodissociation of H^O 
at these heights produce sufficient amounts of H causing the ledge to 
appear in the morning. Subsequent losses due to various chemical reac- 
tions later in the day bring down its concentration to the calculated 
value. The rate of production of H by photodissociation of 1^0 near 75 km 
is not adequate to yield concentrations as given by the noon profile, and 
hence the presence of the minimum in the profiles at these altitudes. 

The valley, however, is later filled in by the transport processes. 

These results are generally in agreement with the results of Thomas 

and Bowman [1972] and Hunt [1971b; 1973], The [Hj peak in Thomas and 

8 -3 

Bowman's profile occurs near 83 km with a value of about 5 x 10 cm , 

8 

while in Hunt's profile, it occurs near 82 km with a value about 2 x 10 
cm . In the present calculations with the low D model, the peak 

O' 

8 - 3 

ab indance of H obtained is 4 x 10 cm near 85 km, showing good agreement 
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Figure 3.17 The height distributions of H at noon, 60° zenith angles and pre-dawn, 
calculated using the diurnal model with high eddy diffusion coefficient 
profile. 
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with the results of the previous workers. However, the only experi- 
mental estimate of the H concentration in the mesosphere carried out using 

7 -3 

OH airglow data give values about 3 x 10 cm in the height range 80 and 
90 km [ Evans and Llewellyn s 1973]. This is about an order of magnitude 
lower than the corresponding theoretical values discussed above. 

The diurnal variations of OH and HO 2 are shown in Figures 3.19 and 

3.20, respectively, for the low D model only. The nature of the diurnal 

& 

m 

variation for the high diffusion model is also similar, and hence it is 
not shown separately. For OH in the region below 70 km, the variation 
between the two x = 60° profiles is not appreciable, at the most a factor 
of 1.5 except below 35 km. Above 70 km, up to about 80 km, where the OH 
curve has a moderate slope, the evening values exceed the morning values 
by a factor of about 4. Above 80 km, the diurnal variation is negligible. 

The more prominent feature here is the behavior of the pre-dawn 
profile of OH. The almost negligible loss rate in the night causes the 
OH concentration to build up by pre-dawn, below about 70 km. The reaction 
between OH and 0^, when taken into consideration is the major loss process 
of OH in the night. However, in this study, this reaction was not in- 
cluded, because of the lack of experimental evidence in favor of this 
reaction [Sdhiff 3 1972] . The reaction between OH and HO 2 and that between 
two OH molecules do not contribute significantly to the loss of OH in the 
night, in view of the long time constants of these reactions. 

A similar diurnal variation is also observed in the case of HO 2 
[Figure 3.20). In the two x - 60° profiles, the difference occurring 
below 68 km is rather small, while the noon profile gives the largest set 
of values. Between 68 and 80 km, the evening values dominate while the 
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Figure 3.19 The height distributions of OH at noon, 60° zenith angles and pre-dawn, 
calculated using the diurnal model with low eddy diffusion coefficient 
profile. 
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Figure 3.20 The height distributions of H0 2 at noon, 60° zenith angles and pre-dawn 
calculated using the diurnal model with low eddy difussior coefficient 
profile. 
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morning values get reduced below the noon values. In this region, the 
morning and evening values again differ by a factor of 4. Around 80 km, 
the morning values display an irregularity caused by similar variations 

3 

found in 0( P) and 0^ profiles. Above 85 km, there is no significant 

diurnal variation. The pre-dawn profile displays a much larger variation. 

However, unlike in the case of OH, the stratospheric values of HO 2 at 

nighttime do not exceed the noon values. In the upper stratosphere H0 2 

3 _3 

has a nighttime concentration of about 5 x 10 cm . It has a deep valley 
near 75 km, and reaches the daytime values above 80 km. 

In the previous calculations of 0-H species concentrations by other 
workers, the inadequate knowledge of some of the rate constants (Ryy and 
R 3q ) gave rise to a certian degree of uncertainty in their height distri- 
butions. However, recent determinations of more exact values for these 
rate constants have enabled the calculation of more accurate height dis- 
tributions, in particular for OH and H0 2 . Further, the assumption of 
zero -rate constants for the reactions of OH and H0 2 with 0^ causes 
enhancement of these constituents in the stratosphere during nighttime. 

Generally, the daytime profiles of both OH and H0 2 agree with those 
of the previous workers. The more or less constant values up to about 
75 km, and the rapid drop thereafter are some noted features present in 
these profiles. The calculated values of OH are also in agreement with 
the measured values of stratospheric OH by Anderson [1971a, 1971b]. No 
direct measurements are available for any of the other 0-H species at 
mesospheric heights. 

3.6.3 Nitrogen and oxygen-nitrogen constituents. The results of 
the previous set of calculations of the 0-N products were employed as 
initial values in the calculations of the diurnal distribution of these 
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species. The values obtained for the 10th day noon are given in Figures 
3,21 and 3.22 corresponding to the high D g and low D g values, 
respectively. 

As expected, the new values of [^0] in the mesosphere, are higher 
than the corresponding values obtained in the previous calculations. 

During nighttime, N 2 0 is not subject to any losses due to photodissocia- 
tion and hence the mixing ratio tends to reach the limiting values 
determined by the lower boundary conditions. The concentrations in the 
mesosphere and the lower thermosphere are controlled by eddy diffusion, 
and as a result, the profile under low diffusion gives values at 120 km 
ttfhich are higher by an order of magnitude compared to those given by the 
high D g model. 

In the previous calculations, it was seen that eddy diffusion plays 
an important role in the determination of the height distribution of NO 
in the mesosphere. This is because of the long time constant NO has 
against chemical losses, which is over two days between 70 and 90 knu 
Below 70 km, it comes down due to the increasing concentration of 0^, which 
oxidizes NO into NO,,. The two [NO] curves in Figures 3.21 and 3.22 demon- 
strate this. The two curves are identical below about 65 km, while above 
this height, the high diffusion model gives higher values than the low 
diffusion model values. The difference in the concentrations is within a 
factor of 3. The apparent coalescing of the values near the upper 
boundary is due to the fixed boundary value, and does not represent 
the real situation. 

The production of NO in the thermosphere depends on the branching 

f 4- 

ratios of the NO recombination reaction and the reaction between NO and 
0 2 . Different values between 0.5 and 1.0 were tried for these parameters 
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Figure 3.22 The height distributions of nitrogen and oxygen-nitrogen constituents at noon, 
calculated using the diurnal model with low eddy diffusion coefficient profile. 


and their effect was noticeable only in the height range between 95 and 
105 km. The differences obtained with these values of the branching 
ratios appear to be within a factor of 2 at 100 km. The profiles shown 
correspond to y = 0.75. 

As mentioned in Section 2.5 of the last chapter, NO and N0 2 are 
strongly coupled through the reactions 

R 38 N0 + °3 * N °2 + °2 

R 41 N0 2 + 0 * NO + 0 2 , 

which are actually effective as a catalytic process for the conversion of 
o7one and atomic oxygen into molecular oxygen. The reaction R 3g acts as 
the main loss process for NO, and also as the main production process for 
N0 2 below 70 km. The other reaction R 41 is the main loss process for N0 2 
over the entire height range during the daytime. This reaction also acts 
as the main production process for NO below about 95 km. In view of this 
strong coupling, when the individual time -dependent equations are solved 
for NO and N0 2 in sequence, any small errors introduced at the beginning 
could get magnified leading to an instability in their solution. This 
is particularly important in the height range below 70 km. 

During the daytime, N0 2 is in photochemical equilibrium with NO 

and, therefore, its continuity equation can be used to eliminate terms 
involving [0] and [0^] which appear in the continuity equation for NO. 

This procedure results in removing the instability mentioned above. The 
profiles shown in Figures 3.21 and 3.22 have been obtained after incor- 
porating the above procedure in the height range below 70 km. This 


simplification, however, was not introduced in the nighttime calculations 

as during this period the condition of equilibrium for NC^ becomes invalid. 

Strobel [1971b, 1972a] who studied the NO at D- region heights showed 

that its concentration in this height range is sensitive to the eddy 

diffusion coefficients used, and to the lower boundary conditions. Under 

high eddy diffusion he found the profiles to have broad minima between 65 

7 -3 

and 70 km with values over 10 cm , while the low eddy diffusion models 

6 3 

gave minimum values in the range (2-4] x 10 cm . More recently, Brasseur 

and Niaolet [1973] showed that the NO profiles have minima between 70 and 

5 6 -3 

80 km with the values extending from 10° up to 5 x 10 cm under dif- 
ferent values of eddy diffusion coefficients and nitrogen atom production 
rates. In the [NO] profile obtained by Hunt [1973], however, no such 
minimum is observed in the mesosphere. 

In the present study, the minima in the two noontime profiles of [NO] 
7 7-3 

have values 1.4 x 10 and 4 x 10 cm , both lying at 80 km. These pro- 
files, which show identical values below 65 km, have broad maxima between 

10 -3 

45 and 50 km, with values 8 x 10 cm However, in view of the uncertain 
lower boundary conditions, these latter values cannot be relied upon as 
representing the real situation. 

Though these profiles differ slightly from the theoretical profiles 
of the other workers, they are in good agreement with the available ex- 
perimental profiles. Meira's [1971] profiles of [NO] show minima around 
7 -3 

1.5 x 10 cm at 84 km. The more recent measurement by Tisone [1973] 

7 -3 

shows a height distribution having a broad maximum of 5.5 x 10 cm 
between 80 and 100 km. In view of the inherent weaknesses of the tech- 
niques employed in these measurements they are not, however, considered 
to be reliable for determining D - region NO concentrations [Thomas, 1974]. 
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Figure 3.23 shows the diurnal variation of [NO] for the low diffusion 
case. Above 70 km, the variation is negligible. This is expected because 
of the effect of the transport processes in this height range. Between 
the two x = 60° curves, the morning profile is nearly a factor of 10 below 
the noon profile at the stratopause, tfhile the drop in the afternoon 
profile is only very small. At pre-dawn, the [NO] profile vanishes below 
40 km. 

These results show that the diurnal variation of the NO concentra- 
tion appears only below the region where the Ly-a ionization of NO is 
important. Therefore, the variation of the [NO] does not contribute 
much to the diurnal variation of the ion-pair production rates in the 
D region . 

The values of [NO,,] obtained in the new set of profiles show marked 

changes from the constant x profile above 80 km as well as below about 

50 km. The changes taking place in the higher altitude range are due 

3 

to the changes in the [0( P] ] and [0^] values, whereas in the strato- 
sphere the diurnal cycle causes enhancements in the N0 2 concentrations. 

This is because at nighttime both loss processes of N0 2 in the strato- 
sphere, vis. reaction with 0( P) and the photodissociation, are absent. 

In the height range between 60 and 75 km, where NC> 2 plays an important 
role in the ^-region negative-ion chemistry, the diurnal cycle does not 
introduce any significant changes in the noon profiles. 

Unlike in the case of oxygen only or oxygen-hydrogen species, the 
calculation of nitrogen and oxygen-nitrogen species concentrations poses 
some additional problems. The odd nitrogen species, N and NO are pro- 
duced in the thermosphere, where the ionic chemistry plays a significant 


ALTITUDE (km) 


LOW D, MODEL 

NOON X s 18° 

0753 hrs X=60 

1647 hrs X= 60 

PRE-DAWN 


NUMBER DENSITY (cnr*) 

Figure 3.23 The height distributions of NO at noon, 60° zenith angles, and pre 
dawn, calculated using the diurnal model, with low eddy diffusion 
coefficient profile. 




137 


' 


■ ■ 
i- 

i- 

i: 

S'. 

t- 

i 

i; 


role. In view of the many uncertainties present in the ionic chemistry, 
the rate of production of those species cannot be determined accurately. 
Secondly, the sinks for the odd nitrogen species in the stratosphere have 
not been correctly estimated so far. Finally, the wide range in the time 
constants involved in the loss rates of odd nitrogen species makes the 
computer simulation of the diurnal behavior of this quantity a little 
more difficult task than in the case of the other constituents. 



138 


4. ION-PAIR PRODUCTION RATES 

Several sources contribute to the ion-pair production in the D region. 
These sources include radiations of solar origin as well as others such 
as galactic cosmic rays and precipitating electrons which come from the 
outer radiation belts. The solar radiations that are capable of ionizing 
any atmospheric gas and penetrating down to 100 km have wavelengths below 
100 A, and also between 1000 and 1300 A. However, as mentioned in the 
last chapter, wavelengths outside of these limits have also been included 
in order to extend the region of ionization up to 120 km for use in the 
calculation of neutral ntirogen species concentration. 

The solar radiations get attenuated in the atmosphere mainly due to 
absorption by oxygen. Molecular nitrogen also plays a significant role 
in the attenuation of incoming radiation, particularly the X-rays below 
100 A, in view of its high abundance in the atmosphere. The intensity of 
these solar photon radiations at D- region altitudes drops with the 
increase of solar zenith angle, as the thickness of the atmosphere 
traversed by these radiations increases. 

The calculation of the ion-pair production rates due to these solar 
EUV and X-rays and their variations with the zenith angle are presented 
in the next two sections. 

Section 3 presents the ionization rates due to sources of 
non-solar origin, viz. galactic cosmic rays and precipitating electrons. 
Unlike in the case of solar radiation, these do not have a solar zenith 
angle variation. In the last section, a discussion of the resulting ion- 
pair production rates is given. 
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4.1 Solar EUV Radiations 

The ionization thresholds of most of the gases present in the earth's 

atmosphere correspond to wavelengths below about 1330 ft. Table 4.1 lists 

these ionization potentials and their threshold wavelengths. Only 0 2 , 

1 

0 9 ( A ), N0 9 , and NO have sufficiently low ionization potentials so as 

to be ionizable by wavelengths greater than 1000 ft. The most significant 

contribution for the C> 2 ionization in this region comes from the solar 

Ly-B line at 1025.7 A, The several spectral windows appearing between 

1050 and 1300 ft in the oxygen absorption spectrum (Figure 2.4) makes the 

remaining three constituents potential sources of 27-region ionization. Of 

these, N0 2 has a very low abundance at Z7-region heights and therefore does 

not make a significant contribution as an ionization source. Of the 

other remaining two constituents, NO has been long recognized as the 

1 

major source of 27-region ionization, while 0„( A ) is considered impor- 

z g 

tant as an 0 2 + source in the mid -27 region. The contributions from these 

as well as from 0 2 are discussed in the next section. 

4.1.1 NO ionization. Even though NO is only a minor constituent 

at 27- region heights having only an abundance of a few parts per million 

or less, the combination of the high intensity of the solar Ly-a line at 

1215.7 ft and the low absorption of this line by atmospheric 0 2 makes this 

a major source of ionization in the 2? region. As discussed in Section 

2.2.1, the intensity of this line at the top of the atmosphere corre- 

11 -2 -1 

sponding to July 24, 1968 is 3.7 x 10 photons cm sec or 6.0 ergs 
cm sec . According to rocket measurements [Goldberg and Aikin , 1971], 
the unit optical depth of this line lies near 77 Ion for low zenith angles, 
and its intensity is completely attenuated by the time it reaches 70 km. 


Constituent 

Ionization 

potential 

CeV) 

Threshold 

wavelength 

CM 

N 2 

15.6 

790 

H 2 

15.4 

803 

N 

14.5 

852 

CO 

14.0 

882 

co 2 

13.8 

895 

0 

13.6 

910 

H 

13.6 

910 

OH 

13.3 

930 

n 2 o 

12.9 

960 

°3 

12.8 

970 ' 

h 2 o 

12.6 

980 

°2 

12.0 

1030 

°2 Cl V 

11.1 

1118 

NO, 

9.76 

1270 

NO 

9.27 

1330 
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The ion-pair production rate due to this line is given by 

<7(1216A) = [NOjo^CNO) 1,(1216 A) expC-rCx^)) (4.13 


where 

a. (NO) is the ionization cross section of NO at 1216 A, and 
I o (1216 A) is the unattenuated intensity of the solar Ly-a line. 

Since this line is absorbed mainly by molecular oxygen in the atmosphere, 
the optical depth factor is given by 

t(x, 2 ) = a a (° 2 ) sec x | [0 2 K 3 3 dz (4.2) 

3 

where 

o^(0 2 ) is the absorption cross section of 0 2 at 1216 A. 

-20 2 

As discussed in Section 2.2.2, a (0„) is taken to be 1.0 x 10 cm . 

- IS 2 

The value of o.(NO) at 1216 A adopted is 2.0 x 10 cm . 

The required [0 2 ] and [NO] values are taken from the diurnal model 
of the neutral atmosphere presented in the last chapter. Above about 
70 km, transport effects control the [NO] profile and therefore its varia- 
tion with x in this region is negligible. Hence, the variation of 
<gr (1216 A) with x is mainly due to the variation of the value of sec x 
appearing in the optical depth factor. 

The production rates obtained using the previously calculated [NO] 
profiles are shorn in Figure 4.1 for different x angles. These curves 
show that the zenith angle variation becomes significant for heights below 
about 90 km only. Above this height the optical depth factor is very 



ION-PAIR PRODUCTION RATE (cm' 3 sec' 1 ) 

Figure 4.1 The ion-pair production rates due to Lyman -a (1216 A) ionization o£ NO at zenith 

angles 18° (noon) and 60°. ^ 

•e* 
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small and the variation of the sec x factor has no influence on the 
intensity of the downcoming radiation. As the optical depth factor begins 
to grow the variations become very prominent, and even exceed an order of 
magnitude between x = 18° and x = 60° below 72 km. 

4.1.2 0 9 (*A ) ionization. The metastable 0_(*A ) was suggested 

as an important source of P-region ionization by Hiarben and MoElroy 

4 * 

[1968], Its importance is in the production of 0^ in the mid-P region 
rather than in the total ion-pair production. The formation of water 
cluster ions in the mesosphere requires the presence of an 0^ source, 
according to the original theories put forward to account for the presence 
of these cluster ions [Fehsenfeld and Fevguson 3 1969], Even with the 
current theories on the production of these cluster ions, it is easier 
to explain the presence of ions such as H^O* and if the 0 2 

production rate becomes enhanced as in the case of a PCA event [Naroisi 
et al.t 1972b]. 

+ 1 

The original values of 0 2 production rates due to 0 2 ( A ) ioniza- 

9 

tion were of the correct magnitude to account for the presence of the 
cluster ions around 80 km. However, a re-evaluation of this production 
rate by Huffman ei at. [1971] , taking into account the C0 2 absorption of 

j 

solar radiation between 1027 and 1118 A responsible for 0 2 ( A ) ioniza- 

9 

tion, showed that the 0 2 + production rates are reduced significantly 
from the previously calculated values. This modified ionization rate is 
calculated here using the formula given by Paulsen et at. [1972] 


q(0 2 ( 1 A^]) = O^A^) jo. 549 x 10~ 9 exp(-2.406 x 10“ 2 ° N(0 2 ]] 
+ 2.614 x 10" 9 expC-8.508 x 10" 20 N(0 2 ))j 


(4.3) 
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where 

N C°2-* column density of 0 2 molecules along the solar ray path. 

According to Paulsen et al, the combined uncertaintly due to any 
uncertainty in the [C^] and [C0 2 ] profiles is at most 15%, Larger un- 
certainties, however, occur in the solar flux values and in the 0 2 
and C0 2 absorption cross sections used to derive this formula. The rms 
error due to all these data has been estimated to be less than 10% with 
a maximum error less than 33%. 

In the above formula, values obtained from the diurnal model calcu- 

1 1 
lations are used for 0 ? ( A ). These calculations have shown that 0 7 { A ) 

*9 ^-9 

has a significant diurnal variation, with the values at 80 km increasing 

by as much as an order of magnitude from post-sunrise (x = 85°) to noon 

+ 1 

Cx = IS 3 . As a source of 0 9 in the mid-C region, 0 7 ( A ) has to 

z z g 

compete with 1-8 A X-rays and precipitating electrons. Of these, the 
X-ray ionization rates have a solar zenith angle variation caused by the 

optical depth variation only, while 0 2 (*A ) ionization rates have a solar 

9 

zenith angle variation caused by both optical depth variation and 

[O^C^A )] variation. On the other hand, ionization rates due to precipi- 
^ 9 

tating electrons do not have any solar variation at all. Hence, the con- 
tributions made by each of these sources vary at different rates during 
the day. Figure 4.2 illustrates the calculated 0 2 production rates due 
to O^C^A ) ionization at different solar zenith angles. 

4.1.3 0 2 ionisation. In the spectral region below 1030 ft, the 
threshold wavelength for 0 2 ionization, there are some wavelength bands, 
particularly in the Lyman continuum, for which N 2 offers only a small 
absorption cross-section. Consequently, these wavelengths penetrate to 
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heights below 100 km and contribute to lower E region ion-pair produc- 
tion. The biggest contribution comes from the Ly-B line at 1025.7 A, 
and to a lesser extent from the C III line at 977 A. The unattenuated 
intensity of this Ly-B line has been measured over a long period of time 
and was found to vary by about a factor of 1.5 during a solar cycle [ Hall 
&t al. s 1969]. However, no measurements are available for this line on 

July 24, 1968. The closest day on which a Ly-B measurement has been re- 

9 

ported is for November 21, 1968, the flux on this day being 3*7 x 10 
-2 “1 

photons cm sec [Hall et al* j 1969]* The Ly^a flux reported for this 
day by Timothy and Timothy [1970] is 3.5 x 10 photons cm” sec” . 

Assuming the two lines, Ly-a and Ly-B undergo similar variations main- 
taining the same flux ratios, one can deduce the Ly-B flux if Ly-a 
flux is available for that day. Since the value adopted for the Ly-a 

flux for July 24, 1968 is 3.7 x 10^ photons cm”^ sec”^, it follows 

9 -2-1 

that Ly-6 flux is 3,9 x 10 photons cm” sec” . This is within the 
range of values given by Hall et al. [1969] for the Ly-$ flux.' 

For che 977 A line and others in the Lyman continuum, fluxes reported 
by Hinteveggev [1970] have been adopted. The absorption cross sections 
and ionization yields employed here are those given by Ohshio et al. 

[1966], The optical depths are controlled by absorption, and hence 
its variations are due to changes of sec x only. The ion-pair production 
rates are calculated using equations similar to equations (4.1) and (4.2) 
for solar zenith angles 18° and 60°. These results, illustrated in 
Figure 4.3, show that between 95 and 110 km, EUV ionization of 0^ is the 
main source of ion-pair production. 



ION-PAIR PRODUCTION RATE (cm' 3 sec" 1 ) 

Figure 4.3 The icn-pair production rates due to Lyman-3 (102? A) and C III (977 A) ionization 
of C>2 at zenith angles 18° (noon) and 60°. 
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4.2 Sola? X-Rays 

The penetration of solar X-rays into the atmosphere depends mainly on 
the absorption cross sections offered by the major constituents 0 2 and N 2 
to X-rays, The absorptit cross sections for these gases are obtained 
from experimentally determined mass absoiption coefficients, p/p. The 

f 

t 

| values of this quantity measured by several workers were assessed and the 

| 

i weighted means were tabulated by Henke and Elgin [1970], The absorption 

| cross sections for 0 2 and N 2 deduced from these tabulations are shown in 

Figure 4.4. An important feature of these plots is the existence of a 
sharp discontinuity at 23.3 A in the case of 0 2 , end et 30.9 A in the 
■ case of N 2 * This critical wavelength, A K is associated with the ejection 

of a K electron from the atom, and sometimes referred to as the K absorp- 
! tion edge. These absorption cross sections are represented by the following 

empirical formulas in each of the three given wavelength intervals. 

i A < 5 A 


a(0 2 ) = 1.67 x 10- 22 A 3 

C4.4a) 

o(N 2 ) = 9.69 x 10" 23 A 3 

C4.4b) 


5 A < A <\ R 


o(C 2 ‘! = 1.70 x 10 -22 

•z _24 4 

A° - 3.35 x 10 A 

(4.5a) 

a(N 2 ) = 1.02 x 10" 22 

A 3 - 1.72 x 10~ 24 A 4 

C4.5b) 
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i. 
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\ R < X < 100 A 


a(0 2 ) = 2.65 x 10" 23 X 2 * 46 

(4.6a) 

c(N 2 ) = 1.97 x 10~ 23 X 2,40 

(4.6b) 


where A is in Angstroms. The values for a(0) are taken as half of 
a(0 2 ). These expressions are similar to those reported by Slider [1969]. 

The presence of the discontinuity in the absorption spectra of 0 2 
and N 2 divides the spectrum below 100 A into two groups, one group lying 
below A^. and the other lying above A^. Referring to Figure 4.4 one 
sees that the cross sections just above X^ have the same low values as 
the cross sections near 10 A. This similarity in the cross sections makes 
it possible for both of these wavelength regions to penetrate to equal 
depths. The wavelengths just above X^. and near 10 A reach altitudes 
between 90 and 100 km, while the wavelengths below about 8 A whose cross 
sections are much lower reach the ZX-region altitudes. The ionization 
effects of these two wavelength bands are therefore considered separately. 

4.2.1 X-r>ays helm 10 A. X-rays in this band of the solar spectrum 
are produced in the solar corona which is a hot plasma with temperatures 
over a million degrees Kelvin. Several ionic processes take place in this 
region giving rise to different components of the spectrum. The free-free 
transitions (bremsstrahlung-type emissions) and free-bound transitions 
(radiative recombination emissions) cause the continuum while bound electron 
transitions between different energy levels in the ions give rise to the 
line spectrum. Both the theory [Elwert, 1961; MandeVstam , 1965] and 
observations [Rugge and Walker, 1968; Evans and Pounds, 1968] show that 
the region below about 15 A is mainly a continuum spectrum, except for a 
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few lines scattered towards the high end of the band. The relative impor- 
tance between the line spectrum and the continuum spectrum depends on the 
coronal temperature, the line spectrum becoming more prominent at high 
temperatures such as those prevailing during solar flares. Under non- 
solar flare conditions the line spectrum below 10 ft can be easily neglected. 
Photographs of the solar disc taken with rocket-borne pin-hole cameras 
and X-ray sensitive plates show that these X-rays are emitted from highly 
localized regions spread over a small fraction of the total disc area 
[j Friedman j 1963]. Observations with selective spectral sensitivities 
have shown that the emitting regions become more localized with the 
decrease in the wavelength. Theoretical interpretation of the observa- 
tions requires a higher coronal temperature to explain the emission of 
shorter wavelengths while longer wavelengths can be accounted for in 
terms -*f lower coronal temperatures. 

The solar spectrum below 8 ft has been regularly monitored over the 
past two decades by satellite-borne detectors. These detectors respond 
to two wavelength bands, 0.5-3 A and 1-8 ft [Kreptin^ 1961, 1965], 

In order to convert the detector currents to flux values it has been 
necessary to assume spectral distributions over the measured bands . 

The procedure adopted at the inception of this series of measurements is 
to ass um e a black-body distribution given by Planck's radiation formula 
which expresses the energy per unit wavelength interval as 

, oJ\T -1 

= K 1 A (e - 1) , (4.7) 


i 
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where 

X. is a constant, 

a 2 is the second radiation constant, 1.439 cm K, and 

T is the equivalent black -body temperature. 

Fixed temperatures of 2 x 10^ K and 10 x 10^ K have been used in the 
conversion of 1-8 A and 0.5-3 A data, respectively. This assumes that 
the X-ray flux variations are due solely to variations of the nature of 
the emitting regions of the corona, rather than due to any temperature 
variations. Hoivever, it is believed that enhanced emissions are accom- 
panied by hardening of the spectrum caused by an increase in the tempera- 
ture [Culhane et al. 3 1964]. 

This procedure of using a black -body model in the conversion of 
satellite data into flux values is still continued even though more 
realistic models for the energy distribution in the bands are presently 
known [Elwert^ 1961, MccndaZ'stam , 1965]. According to these models, the 
emissions due to bremsstrahlung and radiative recombination could be 
adequately represented by [fforem, 1970] 

oJ\T 

E{\) = Ke 2 e C4.8) 

where T is the electron temperature in the corona, and K is a function 
of the emission measure B which is given by 



V 


where is the electron concentration and the integration is performed 
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over the respective emitting regions. The value of K also depends to 
a smaller extent on X and T . Better agreement has been found with 
spectra measured during enhancements with the above distribution than 
with a black-body type spectrum. 

In order to compute the ion-pair production rates due to these 
X-rays, one has to cons i-uct the spectral distributions from the measured 
broad-band fluxes. The main difficulty here is the estimation of the 
correct coronal temperature and the emission measures corresponding to 
the relevant bands. Horav. [1970], however, has developed a method to 
work out these quantities in terms of satellite-measured fluxes using 
a numerical technique which is somewhat tedious. 

Using a more simple procedure, it is also possible to construct the 
Spectral distribution over the entire range 0-10 A in terms of the 
measured fluxes in the 0.5-3 A and 1-8 A bands. Here, the energy dis- 
tribution per unit wavelength interval is given by an empirical formula 
[Roue et al. s 1970] 


d'(A) = AX n (4.10) 

where A and n are constants , n determines the shape of the spectrum 
while A determines the absolute flux, n is obtained from the ratio of 
the measured fluxes over the two bands using the formula 


_ log(<H8 A)/cK3 A)) . 

log 2.667 


(4.11) 


where <J>(8 A) and <j>(3 A) are the two measured fluxes over the bands 
1-8 A and 0.5-3 A, respectively. The constant A is obtained using 
the flux value over one band. 


^ - Ctt+jQ j(8 A) 
“ 8 n+ ^ - 1 


(4.12) 


The ratio <J>(8 A)/<|>(3 A) is generally of the order of 100 under quiet 
conditions, hut varies by as much as a factor of 5 under extreme condi- 
tions. In two proportional counter spectra measured under conditions of 
low and moderace solar activity [ Pounds ^ 1970], the values corresponding 
to this flux ratio were found to be 1780 and 43, respectively. The 

(1-8 A) flux deduced from the solar active day spectrum, smoothed and 

-3 -2 -1 . 

shown in Figure 4.5, has a value 1.14 x 10 erg cm sec . With this 
flux value and the <|>(8 A)/H3 &) ratio of 43, the spectrum was re- 
constructed using equation (4.10). These results are also shown in 
Figure 4.5. The excellent agreement between these values and the measured 
spectrum shows thi.t equation (4.10) is a valid representation of the X-ray 
spectrum below 10 A. 

These X-rays ionize all constituents present in the atmosphere. In 
calculating their ion-pair production rates, the ionizing effects on N 2 , 
0^, and 0 have been considered, with the production function given by 


f2 

4(1-10 A) - 7 . [X] i j n i (A)o^ (X)J o (X) e -T<:x ’ s:) 

X 1 


(4.13) 


3 

where i denotes the constituents N 2 , C> 2 , and 0( P) • n^U) is the ioni- 

til 

zation yield of the i constituent, and is the energy distribu- 

-2 -1 

tion function given by equation (4.10) expressed in photons cm sec . 
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Figure 4.5 The spectral distributions of X-rays below 10 K correspond! n 
to different flux values in the 1-8 A and 0.5-3 A bands. The 
broken curve shows the spectrum measured by Pounds [1970] . 
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Equations (4.4) and (4.5) give the values of cr^ over the wavelengths 
from 1 to 10 A. 

The optical depth factor is governed by absorption due to all con- 
stituents, and is given by 

CO 

t(x,s) = secy l ot (X) f W^(s) dz . (4.14) 

i •' 

z 

The energies of the incident photons in the 1-10 5 X-ray band 
vary from 12.4 keV to 1.24 keV, which are far in excess of the energy 
required to produce one ion -pair in air. Therefore, each photon is 
capable of producing more than one ion-pair. The first ion-pair is 
produced by the incident photon while the rest are produced by the 
secondary electrons through impact ionization. Photons in the soft 
X-ray bands were found to spend, on an average, about 35 eV in producing 
one high energetic photon, or the ionization yield, can be expressed as 
[Saltier, 1969] 


12400 _ 354 

n ” 35 X(5) " 1(5) 


(4.15) 


This formula, however, does not take into account the efficiency of 
different species in producing secondary electrons through impact ioniza- 
tion, which depends on the impact ionization cross section of each species. 

Ionization yields have also been tabulated by Cfhshio et at. [1966] 
over the entire range of X-ray and EUV wavelengths, using values adopted 
from the results of several workers. In the X-ray wavelengths these 
values show the n(N 2 )/n(0 2 ) ratio is about 1.13, indicating that N 2 is 
more efficient in producing secondary electrons than 0 2> 


157 


More recently, Schlegel [1971] has calculated the ionization yields 
of N 2 , 0 2> and 0 in the upper atmosphere using a Monte Carlo simulation 
of a model ionosphere. According to his results n is a function of both 
the altitude and the solar flux. The height averaged ionization yields 
for the range 130-300 km, when compared with values reported by Ohshio 
et al. [1966] , show that the two sets of values are in good agreement in 
the case of N 2 and 0. However, a wide discrepancy exists with the 0 2 
values, the Monte Carlo values being lower by a factor nearly 10 for a 
photon of wavelength 10 A. Schlegel attributes this discrepancy as due 
to the wide difference in the 0 2 and 0 densities in the height range con- 
sidered. However, even at 120 km, where the C> 2 and 0 densities do not 
differ very much, Monte Carlo values for 0,, are about a factor of two less 
than the corresponding values given by Ohshio et al. for the wavelength 
interval 10-15 A. The implications of this study on the ionization 
rates near 100 km have not been assessed so far. Since the contribution 
by X-rays is only a fraction of the total ionization at altitudes of 
interest it is expected that a reduction in the 0 2 yield will not cause 
any serious error in the present calculations. In view of this, and also 
for ease in computations, ionization yields tabulated by Ohshio et al. 

[1966] have been adopted in this study. 

In the present study. X-ray fluxes as measured on July 24, 1968 
were used, so that the results could be compared with the electron density 
profiles measured on that day. Figure 4.6 shows the fluxes in the bands 
0.5-3 A, 1-8 A, and 8-20 A as measured with instruments on board the 
SOLRAD 9 satellite. T, e two rocket measurements at x = 18° and x = 60^ were 
taken at 1700 UT and 2136 UT, respectively. It is seen from the recordings 
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Figure 4.6 X-ray fluxes in the bands O.S-^3 R, 1-8 ft, and 8-20 ft as measured 
on July 24, 1968 with instruments on board SOLRAD 9 satellite 
(courtesy NOAA Environmental Data Service) „ 
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that these times are free of any solar flares or any other enhancements. 

Further, the X-ray flux in the 1-8 A band appears to be identical at both 

-2 -1 

instances, being equal to 0.7 raillierg cm sec . The flux in the 

0.5-3 A band on this day, however, is below the threshold value for detec- 

-2 -1 

tion by the instruments, i.e., below 0.01 millierg cm sec . Hence, the 

calculations were done using two values for this quantity; an upper limit 

-5 -2 -1 -5 -2 -1 

of 1 x 10 erg cm sec and a lower limit of 0.5 x 10 erg cm sec . 

These two values correspond to flux ratios of 70 and 140, respectively. 

The spectra corresponding to these flux ratios and 1-8 A flux of 0.7 
-2 -1 

millierg cm sec are also shown in Figure 4.5. 

It should also be noted here that the above satellite measured 
fluxes could be in error to a considerable extent. One source of such 
error is the assumption of a gray-body spectrum for the solar spectrum 
in the data processing procedure. According to Kreplin [1970], the 
published fluxes deduced on the basis of this assumption are unlikely to 
differ by more than a factor of 10 from the absolute flux values for 
wavelengths below 20 A. The other source of error is due to the uncer- 
tainty in the effective passband of the detector windows. Recently, 

Wende [1971] investigated the consistency between data from several 
experiments measuring X-ray fluxes below 20 A and found that both the 
passbands and the conversion factors of the relevant bands have to be 
modified in order to bring agreement among different sets of results. 
Unfortunately, these corrections have been carried out for the early 
experiments only, and the extent of this error as applicable to SOLARD 9 
data cannot be assessed. 

The production rates q (1-10 A) calculated using the above mentioned 
spectral distributions for solar zenith angles 18° and 60° are illustrated 
in Figure 4.7. 


ALTITUDE (km) 


X= 60 *^' 

> 


( I -8 A) =7.0 xIO 4 erg cm’ 2 sec 
|(0-3A)= I.OxlO 5 erg cm" 2 sec 

<£«{ 0 -3 A) = 0.5* 10 5 erg cm 2 sec 


ION-PAIR PRODUCTION RATE (crn^sec" 1 ) 

Figure 4.7 The ion-pair production rates due to X-rays below 10 A at zenith angles 18° (noon) 
and 60°, calculated for the two 1-10 A spectra shown in Figure 4,5 (curves 1 ind 2) 
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4.2.2 SO - IOC 8 X-rays. This region of the spectrum, unlike the 
shorter wavelength region, is composed mostly of emission lines as evi- 
denced by both theory [Mandel 'stem, 1965] and observations [Mans on, 

1967, 1972; Freeman and Jones, 1970], These wavelengths come from 
regions of the corona having low temperatures and are not subject to as 
high a variability as the shorter wavelengths. Most of the emission 
lines between 16 and 100 A are due to transitions occurring in heavy 
elements such as oxygen, carbon, sulphur, and other heavy elements 
including silicon, magnesium and iron. Of these lines, a major contri- 
bution comes from the C IV line at 33.7 A as it lies just above the 
critical wavelength X . Since the absorption cross section of this 
line is low, it is capable of penetrating down to heights close to 90 km. 

In order to compute the ion-pair production rates due to these lines, it 
is necessaiy to find out the distribution of energy among the different 
lines. 

Manson [1972], using a rocket -borne telemetering monochromator, 

obtained a highly resolved spectrum between 33 and 128 A on August 8, 

1967. The 10.7 cm flux on this day was 143 units. The 33.7 A line 

-3 _2 -l 

intensity measured was not more than 2 x 10 erg cm sec . Using 

parachute-recovered photographic plates. Freeman and Jones [1970] 

obtained well resolved spectra on two days. The measurement on March 20, 

1968 (10.7 cm flux = 131 units) which extended from 14 A to about 400 A 

-2 -1 

recorded an intensity of 9 millierg cm sec for the 33.7 A line. The 

other measurement on November 20, 1969 (10.7 cm flux = 189 units) which . 

-2 

extended from IS A to 76 A, recorded an intensity of 7 millierg cm 
sec for the 33.7 A line. On the other hand, the flux intensities of 
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33.7 A measured t>y Argo et at. [1970] using a rocket-borne spectrometer 

under a variety of solar conditions during the period 1964-1968 gave 

-2 -2 2 -1 

values lying in the range 2 x 10 - 6 x 10 erg cm sec These 

values are several factors larger than the values of Freeman and Jones, 
and more than an order of magnitude larger compared to Manson's value. 

Since Manson's spectrum begins only at 33 A, one could argue that any 
end-errors could have affected his measurement at 33.7 A, though he claims 
an accuracy of ± 25% over the entire range of measurements. These measure- 
ments also show that the variability of the other lines between these dif- 
ferent days is much less than that of the 33.7 A line, and that this line 
could be specially sensitive to changes in the solar activity [Manson 3 
1972]. 

Broad-band measurements in the range 44-60 A are also made regu- 
larly by satellite-borne detectors. However, in view of the many assump- 
tions made in the processing of data, these measurements can only be 
taken as an index of solar activity rather than as absolute flux values 
[Kreplin, 1970]. Further, it has been recently reported that the 44- 
60 A detector is also sensitive to radiation in the 8-20 A band, and 
that the published data could be in error by several factors [Kreplin 
et at. 3 1973]. Unfortunately, of the days for which rocket measurements 
are available, 44-60 A flux data are available only one day, viz. 

August 8, 1967. Hence a proper comparison between the rocket-measured 
and the satellite-measured flux cannot be made. 

-2 -1 

The 44-60 A flux measured on August 8, 1967 is 0.306 erg cm sec , 

-2 -1 

while that measured on July 24, 1968 is 0.28 erg cm sec . Considering 
the uncertainties involved, these two measurements show identical condi- 
tions as far as soft X-ray emissions are concerned. Hence, the resolved 
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spectrum obtained on August 8, 1967 by Manson [1972] can be assumed to 
be applicable on July 24, 1968 as well. Therefore, in this study Manson's 
data are adopted for solar flux above 30 ft except for the 33.7 ft line. 

For this line a flux of 7 mill.ierg cm" ^ sec -1 has been adopted, consider- 
ing the variability of the measured data. 

X-rays beyond about 50 ft are significant in the ionization of heights 
above 100 km only and not really important in Z?-region studies. Neverthe- 
less, as mentioned in Section 4.1.3, contributions from all wavelengths 
have been included in order to calculate the ion-pair production rates 
up to 120 km. The results of these calculations are shown in Figure 4.8 
for the solar zenith angles 18° and 60°. 

4.3 High-Energetic Particles 

The sources included in this section are the precipitating electrons 
originating from the outer radiation belts surrounding the earth, and the 
galactic cosmic rays. Both of these sources are considered to have flux 
incident on the atmosphere independent of the solar zenith angle, though 
observations indicate the existence of a small day/night asymmetry. 

4.3.1 Precipitating electrons . The effect of these electrons in 
the earth's atmosphere is most prominent in the polar regions. The ioni- 
zation produced by them in these regions has been investigated in detail 
by several workers [Rees 3 1963]. It was not until satellite detectors 
began to monitor these particles over the earth's atmosphere that they 
were recognized as a potential source of ionization in the D region. 

The ion-pair production by precipitating electrons is naturally a function 
of their initial energies and the energy deposition rate in the atmosphere. 
From laboratory measurements of energy-range relationship, it is found 


no 



Figure 4.8 The ion-pair production rates due to X-rays in the bands 30-44 ft 
and 40-100 ft, calculated for zenith angles 18° (noon) and 60°. 


164 



, -.-A. 


■ -4 


j . i i 


165 


that only the electrons with energies greater than about 40 keV could 
penetrate to Z?-region altitudes, and energies over 200 keV are required 
for penetration down to 70 km. 

The energy deposition rates of these electrons have been worked out 
by Rees [1963] for different distributions of the flux over the hemisphere. 
Using these results, the calculation of the ion-pair production rate can 
be easily determined. The ion-pair production rate due to monoenergetic 
electrons of unit flux is given by 


q(PE) = 


E jn X (g/fl) l a i t X k (a) 
e *0 l a^X]. (ff) 


(4. 16] 


where 

is the initial energy of the incident electron, 
e is the energy required to produce one electron- ion 
pair in air (35 eV] , 

R is the maximum penetration depth of the electron, 

S is the penetration depth at an altitude z 3 

is the maximum range of the electrons, 

X is the normalized energy deposition function, and 

a £ is the ionization efficiency of the i constituent, 

Rees [1963] calculated the normalized energy deposition function for an 
isotropic distribution of the incident electrons over the pitch angles 
from 0° to 80°. His results can be expressed by the empirical formulas. 


X(S , /i?) = 1.1 (S/R) 3 - 3.3 G S/R) + 2.2, 1 > 3/R > 0.1 


1.9 


S/R < 0.1 


(4.17) 
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The penetration depth R of an electron with initial energy E^ n is 

t 

given by the formula 


R = 4.57 x 10' 6 eW 5 (4.18) 

-2 

where E. is in keV and 1? is in g. cm . When the incident electrons 

'in 

have a distribution of energy, the corresponding production rates are 

obtained by integrating over the desired energy range. 

Usually, the detectors measuring the precipitating electrons 

respond to the total flux above a certain minimum threshold energy . 

Therefore, the energy spectrum is expressed in terms of this total flux 

N (E . ) . Measurements carried out in the past [Potemra and Z rnuda, 
mm' r 

1970] have shown that the total flux of precipitating electrons can be 
expressed by a power exponential form 

m>E) - » 0 c>W ce/W y 


where 

N(>E) is the flux of electrons with energies greater than E a and 
'i is a constant representing the hardness of the spectrum. 

It is also sometimes possible to express H?(>£) by a simple exponen- 
tial form 


E“E 

H>E) = (>i W ex P ( ” 


(4.20) 


where E is the e-folding energy, a parameter representing the hardness 


of the spectrum. 
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Because of the scarcity of the observations and the wide scatter of 
i\ie available data, it is difficult to determine unambiguously the correct 
form of a general expression for the precipitating electron spectrum that 
can be applied in instances such as ionospheric calculations. In general, 
the exponential spectrum shows smaller fluxes of electrons with higher 
energies. In the present calculations, a power exponential spectrum is 
preferred as it would give higher ionization rates at lower altitudes 
than a simple exponential spectrum with the same N . According to 
O’Brie- 1 [1964], satellite measurements of the average flux of precipi- 
tating electrons with energies greater than 40 keV were f vund to vary 
from 10 to 10 electrons cm” sec” ster” at midlati ,udes. On the 

other hand, rocket measurements of midlatitude electron precipitation 

-2 -1 -1 

show low fluxes in the range 4-7 particles cm sec ster for 
electrons with energies greater than 40 keV under magnetic quiet condi- 
tions [O’Brien et at.* 1965; Gough and Collin 3 1973], Potemva and 
Zmuda [1970] who compared the C-region nighttime ionization rates due 
to scattered Ly-a radiation and precipitating electron found that the 
ionization rates due to precipitating electrons with #^(>40 keV) equal 
to 300 and y = 3 just exceed those due to scattered Ly-a radiation. 

In the present study, it is necessary to look for an 0 2 + source 
between 70 and 80 km that would partly account for the cluster ion forma- 
tion. In order to investigate the precipitating electrons as a potential 

source of O 0 + ions in this height range, ion-pair production rates have 

-2 -1 -1 

been calculated using total fluxes of 300 electrons cm sec ster 
and y varying from 2 to 4. The results are shorn in Figure 4.9. 

4.3.2 Galactic cosmic rays. Cosmic rays penetrating down to D- 
region altitude? are composed mostly of protons possessing energies 



ION-PAIR PRODUCTION RATE (cm" 3 sec'') 

Figure 4.9 The ion-pair production rates due to galactic cosmic rays, and precipitating electrons 
having energies greater than 40 keV (N(> 40 keV) = 300), and hardness indices of 2, 3, 
and 4. 
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greater than about 10 MeV, and the galactic proton flux spectrum in this 
energy range peaks around 300 MeV, The proton spectrum for energies 
above this value is generally expressed in a power exponential form, 
similar to that given in the last section for precipitating electrons. 

Velinov [1968] has worked out the ion-pair production rates due 
to galactic cosmic rays using laboratory derived energy deposition func- 
tions for relativistic particles, and using assumed energy distributions 
for the incoming particles. He showed that this ion-pair production 
rate has a geomagnetic latitude (X^) dependence of cos~ \ m when the 
rigidity spectrum is of the power exponential form with the exponent 
equal to 2.5, arid the geomagnetic threshold of rigidity is given by 

J? =14.9 cos 4 X . In order to calculate the ion-pair production rates 
am 

theoretically, one has to know correctly the energy distribution func- 
tions, the composition of the incident particle flux, etc. In spite of 
several years monitoring of the incident cosmic ray intensities, no 
proper picture has emerged so far t.iiich will represent these fluxes 
correctly. 

In view of these difficulties it has been custormary to obtain the 
ion-pair production rate using the semi-empirical formula 

q(GCH) = q Q n(z)/n Q (4.21) 


where 


q Is the ion-pair production rate at a reference height and 
is the number density of air at height s^. 

The values of q are usually deduced from direct measurements. 
Using balloon observations, Neher and Anderson [1962] found that the 
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cosmic ray ionization rates at a low altitude increases by a factor of 
10 between latitudes 0° and 60° during solar minimum, and by a factor 
of five during solar maximum. Using these data Slider [1969] has tabu- 
lated factors for determining q(GCR) for different latitudes and solar 
epochs, from which a value of 1.0 x 10 (sec ) is chosen as appro- 
priate for the conditions under investigation. Thus, the ion-pair produ- 
ction rates are given by 


qQSCK) = 1.0 x 10" 17 n(a) 


(4.22) 


The height distribution of this production function is illustrated 
in Figure 4.9. 

4.4 Discussion 

The major sources responsible for the daytime ionization of the 
region below 110 km are shown in Figure 4.10. According to these curves 
the region between 110 and 104 km is ionized at noontime mainly by the 
977 A solar line. Ly-g line ionizes the region immediately below, 
between 96 and 104 km. X-rays between 30 and 40 A, consisting mainly 
of the C IV line at 33.7 A, are responsible for ionizing the next layer, 
between 89 and 96 km. The production function due to X-rays between 40 
and 100 A is also shown here. Though it is never a dominant source at 
these heights, it still makes a significant contribution between 100 
and 110 km. 

The ionization of nitric oxide by the solar Ly-a line becomes the 
major source below 89 km, extending down to about 65 km. Below this 
altitude, galactic cosmic rays become the major source, both at daytime 
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1. PRECIPITATING ELECTRONS 7=2 

2 . X-RAYS I-IOA 

3. LYMAN ALPHA 1216 A 

4. X-RAYS 30-40 A 

5. CIII 977 A 

6. X-RAYS 40-I00A 

7. LYMAN BETA I026A 

8. GALACTIC COSMIC RAYS 


ION- PAIR PRODUCTION RATE (cnf^sec 1 ) 

Figure 4.10 The ion-pair production rates in the D region due to major ionization sources at daytime 
(X * 18°). 
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and nigti^time. Ionization due to precipitating electrons also become 
important here, if one considers a hard spectrum for these particles. 

The rest of the sources considered, viz. EUV ionization of ) 

* y 

and the X-rays below 10 A, both remain as minor sources at all heights 

and zenith angles. For X-rays, two spectral distributions were assumed, 

-5 -2 -1-1 

one with a 0-3 A flux of 1 x 10 erg cm" sec" A , and the other 
having half this value. Both of these spectral distributions give ioni- 
zation rates below 90 km which vary in the range 3-4, corresponding 
to a variation of x from 18° to 60°. 

The solar zenith angle variation of the total ion-pair production 
rates in the D and lower E region are given in Figure 4.11. As men- 
tioned above, the sources ionizing the region between 90 and 110 km 
consist of several components. The solar zenith angle variation of 
each of these sources becomes appreciable in the bottom part of the 
individual production curves. This is because the optical depth factor 
is near zero in the upper portion, and any change in x does not affect 
the attenuation of the incident radiation significantly. Since each sf 
these production curves tail off at different heights, the sum total of 
the variation between two x values due to all of these components remains 
approximately the same. Between 18° and 60°, it is about a factor 
of 1/3 in the height interval 90-100 km. At 90 km, howevor, the con- 
tribution to this difference from the EUV ionization is only 1/10. The 
remainder comes from the Ly-a ionization of NO. 

The decrease of the Ly-a production function with increase in x 
is very small at heights above 90 km. At 90 km it is about a factor of 
1.2. As the optical depth factor increases below this height, the varia- 
tion of the production function also increases. At 70 km, the Ly-a 



ALTITUDE (km) 



ION-PAIR PRODUCTION RATE (cm-3sec-') 

Figure 4.11 The variation of the D region ion-pair production function with the solar 
zenith angle, shown with contributions from two different precipitating 
electron spectra. 
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ionization rate drops by a factor of 1/15 between 18° and 60°. Since 
there is no diurnal variation in [NO] between 70 ard 90 tan as obtained 
in the diurnal calculations, this variation is due solely to the attenua- 
tion of the incident Ly-a intensity. 

Below 70 km, however, the variation due to Ly-a ionization is masked 
by those due to precipitating electrons and galactic cosmic rays. Since 
both of these sources do not have a diurnal variation, their contribution 
becomes dominant at large x angles when the Ly-a contribution drops 
rapidly. Therefore, the diurnal variation exhibited by the total ion- 
pair production curves at a given height in this range depends on the 
relative contributions from each of these sources at that height. 

In the case of the precipitating electrons, the relative contributions 
they make at different heights depend on the spectral characteristics of 
the incident particle flux. In this study, two spectral distributions 
given by y = 2 and y = 3 have been used for the precipitating electron 
flux. In the total ion-pair production curves shown in Figure 4.11, 
the solid lines refer to the y = 2 spectral distribution while the 
broken lines refer to the y = 3 spectral distribution Figure 4.9. 

Since these production rates remain in the background as fixed values ^ 

without varying with x> they control the ionization production at large 
X angles and hence the reduced variation shown by the solid curves. 

The ionization oue to y = 3 precipitating electron.* is much less 
than the Ly-a ionization, even at large x angles, in the height interval 
56-70 km. Hence, the broken lines show a much larger variation with 
increasing x* At 70 km, this is nearly a factor of 7 corresponding to 



a 18°-60 < ’ increase in x* 
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Below 65 1cm, the galactic cosmic xay ionization becomes dominant. 
Hence both sets of curves, corresponding to all values of x» converge 
towards the q (GCR) values shown in Figure 4.9. Thus the total production 
function does not exhibit a significant diurnal variation at these 
altitudes. 
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5. 0-REGION ION CHEMISTRY 

Over the past few years, many experiments have been carried out 
using rocket-borne instruments to measure the 0-region ion composition 
as well as the electron densities. Several attempts have been made to 
interpret these measurements through model calculations, but without 
much success. 

Some of the problems encountered in explaining the 0 region ion- 
composition measurements have already been mentioned in Chapter 1. These 
occur mainly due to our poor understanding of the ion reaction schemes, 
and also due to lack of accurate data on various factors such as reaction 
rate constants and minor neutral constituent densities. 

In this chapter, the 0-region ion chemistry is discussed, with a 
view to develop a unified model for the 0-region ion distribution. In 
Section 1, the previous measurements of the ion composition and electron - 
density profiles are discussed. The relevant ionic reaction schemes and 
the available data on reaction rate constants are described in the sub- 
sequent sections. 

5.1 Measurement of Chained Species Distribution 

5.1.1 Positive-ion composition. The techniques available for 
the in-situ measurement of the 0 region positive- ion composition are 
summarized in Chapter J . Using a quadrupole mass spectrometer described 
there, Navcisi and Bailey [1965] carried out the first detailed measure- 
ment of the positive-ion composition in the 0 region. Contrary to the 
previous theories that this portion of the ionosphere should comprise 
0* and N0 H [Nicolet and Aikin 3 1960], these measurements revealed 
for the first time that the dominant ions below about 82 km are hydrated 
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clusters of the type H + *(H 2 0)^ where n = 1, 2, ... (Figure 1.4). The main 
feature of these profiles is the sudden disappearance of these cluster ions 
HgO + and H + *(H 2 0) 2 above about 82 km, and the appearance of the molecular 
ions 0 2 + and NO* above this height. Below about 75 km, ions with mass 
numbers greater than 45 were also observed, while seme metal ions were also 
observed in small quantities above 85 km. 

The unexpected nature of these results naturally caused much specu- 
lation regarding their validity. There was doubt that these measurements 
represent the ambient ion distributions in the D region. The possibility 
of rocket contaminants giving rise to hydrated cluster ions was always 
there. Further, in the event these cluster ions are really present, the 
possibility of their breaking up into small fragments during the process 
of collection by the probes also cause much doubt about the validity of 
these measurements. Napcisi [1966] however maintained that these clusters 
H^O* and H + • (H 2 <3) 2 were indeed true constituents of the D region rather 
than the products of rocket contamination. According to Navc'Csi [1966] , 
the absolute values of these measurements are correct to within a factor 
of 4, while the relative values are correct to within a factor of 2. The 
possibility of ion fragmentation however was not ruled out. 

Subsequently a large number of positive -ion measurements have been 
carried out as summarized by Narc-isi [1970]. Many of these measurements 
correspond to special events such as eclipses, PCA events at high lati- 
tudes, sunrise and sunset, and not to the ambient daytime D region at mid- 
latitudes. Nevertheless, almost all of these profiles exhibit the same 
features mentioned above. In addition, some profiles reveal the presence 
of N0 + *(H 2 0) in appreciate quantities in the D region. These measure- 
ments are also subject to large errors, factors of two or more, particu- 
larly in the cluster ion region below 80 km. 
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Goldberg and Aikin [1971] measured the D region positive-ion 

composition at solar zenith angles 28° and 53° near the geomagnetic 

equator. They too observed similar ionic distributions with the hydrated 

cluster H 3 0 + , H + *(H 2 0) and H + *(H 2 0] 2 dominating below about 82-83 km. 

They estimated the errors to be within ± 20% on major species, and 100% 

on the minor species. However, it is their opinion that the results are 

mostly qualitative in nature. They noticed that the heavy cluster ions 

tend to break up under shock effects created by the rocket trajectory, 

leaving only the more stable H 3 0 + and H + *(H 2 0] behind. 

Another feature noticed in these ion -composition profiles is the 

presence of 0^ down to 65 km in appreciable quantities. On the upleg, 

-3 

concentrations around 100 cm were detected between 65 and 85 km while 

on the down leg somewhat reduced values were detected. The computed 

profiles however showed steady decline of 0 2 + concentrations with values 

_ 3 

less than 1 ion cm below 75 km. 

Positive- ion composition measurements at high latitudes were carried 
out by Johannessen and Rrankottsky [1972] during daytime (x = 55°). They 
observed the cluster-molecular ion transititn around 85 km. Unlike in 
the previous measurements, H + * (H 2 0) 3 was found to dominate the cluster- 
ion region in a narrow height range around 85 km. Above and below 
this height range, N0 + -H 2 0 and H + *(H 2 0) 2 were dominant among the cluster 
ions, respectively. They also observed high abundance of H + *(H 2 0) 4 near 
85 km. The appearance of these heavy clusters at this altitude range 
was attributed to the very low temperatures prevailing there during the 
time of measurements. These observations were made on both the upleg and 
down leg trajectories, and in both cases the fragmentation of the heavy 
clusters was considered to be small in this altitude range. 


I - A- 


179 

A similar set o£ positive-ion composition profiles was also reported 

by Krankowsky et at. [1972] who carried out measurements at high latitudes 

under nighttime conditions during different seasons of the year. They also 

observed cluster ions H + *(H 9 0) where n = 1, 2, 3 dominating below 82-83 

km, with the relative abundance of the heavier clusters increased on the 

downleg trajectory. This is attributed to the reduced shock effects on 

the downleg, resulting in less fragmentation of the heavier clusters on 

this trajectory. They also observed increased amounts of NO *H 2 0 and 
+ 

H • CH^O) ^ in a flight made on a summer night. This again is attributed 
to the very low temperatues which prevail in the high latitude summer 
nights affecting the reaction rates so as to increase the yield of these 
cluster ions. 

More recently, Johannessen and Krankousky [1974] measured the 
daytime positive- ion composition at t mid- latitude station during summer. 
Their results are also similar to those obtained by previous workers. 

The molecular ion to cluster ion transition was observe! at 85 km. 

N0 + *H 2 0 ions were detected in appreciable quantities between 85 and 
95 km. A distinct feature observed is the high abundance of metal ions 
above 85 km. 

In spite of the many positive-ion measurements that have been 
carried out in the recent past, only a very few have been reported as 
relating to ambient conditions prevailing in the mid-latitude, daytime, 
summer D region. Fortunately, when the measurements are made under quiet 
conditions , the positive-ion composition measurements made near the equator 
and at high latitudes do not differ significantly, except that the relative 
abundance of the heavier clusters tends to increase around the high- 
latitude mesopause. 
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5.1.2 Negative-ion measurements. The model calculations that have 
been carried out in the past to determine the negative -ion distribution 
in the daytime D region have shown that negative ions become important 
only below about 75 Ion [Reid> 1970; Thomas et at, 3 1973]. However , the 
high ambient pressures prevailing in this region and the low ionization 
level, make it difficult to carry out in-situ measurements of the day- 
time negative-ion composition. On the other hand, the negative ion 
density above 75 km increases in the nighttime due to attachment of the 
free electrons to the neutral particles, making it possible to carry out 
rocket measurements in this altitude range. 

Measurements made by Naroisi et at, [1971] at high latitudes 
during nighttime under quiescent conditions indicated the presence of 
molecular ions 0 2 ~ and Cl~, and two others which were identified as L r .«j 
and NOj" between 78 and 80 km. Above 80 km, heavy clusters of the type 
CO^^O and NO^" *H 2 0 were also observed. A second flight made under 
similar conditions showed the presence of heavier ions above 80 km. No 
records were obtained below 75 km in both flights. 

Amctd et at. [1971] also measured the nighttime negative-ion 
composition at a high latitude station, but under enhanced ionization 
during a weak aurora. Their results showed the presence of CO^ , Cl , 

HCOj" and NO^ - between 72 and 77 km. They did not detect any hydrated 
clusters, though many of the other ions observed were of the same type 
as observed by Naroisi et at. [1971], 

The only daytime negative-ion composition results available are those 
obtained during a PCA event at high latitudes, and during totality of an 
eclipse at mid-latitudes [Naroisi et at. » 1972a, 1972b]. The PCA event 


measurements showed the presence of the molecular ions 0 2 ~, N0 3 and 
hydrated clusters of the type NO^'CH^O)^ with n = 0-5 dominated the 
D region below 88 km. 

Here again, the problem of constructing a suitable experimental D 
region negative -ion distribution representing the daytime quiet conditions 
at mid- latitudes is encountered. Most of the nighttime measurements give 
ion densities above 75 km, and are therefore not very helpful in con- 
structing a daytime model for the mid and lower D regions. The two 
daytime measurements give values below 70 km, but these have been obtained 
under abnormal conditions . It is evident from these measurements that 
hydrated clusters NO ~*(H 9 0) are present in the daytime D region. 

It was noted that in many experiments the results are expressed in 
terras of the count rate only, rather than in terms of the ion densities. 
This is because the conversion of the count rate to absolute ion densi- 
ties poses many problems. However, Arnold et at . [1971], expressed 
their results in terms of negative ion densities making use of the 
simultaneous measurements of the total positive-ion density and electron 
density. 

5.1.3 Eleatron-density measurements. Some methods available for 
the measurement of D region electron densities are described in Chapter 1. 
These methods have been recently reviewed by Seahrisi [1974] , who con- 
cludes that the electron- density profiles derived from the rocket measure- 
ments of differential absorption, Faraday rotation and dc probe current 
have the greatest accuracy and the best height resolution. In fact, 
among various in-situ measurements of 0-region parameters, accurate 
measurements are possible only In the case of electron densities 
[Meahtly 3 1974], 
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An extensive series of D region electron- density profiles have 
been obtained using the above technique by Mechtty et at. [1972a] since 

1967. Based on these measurements Meehtly et al. [1972a] have prepared 
two reference profiles for the mid- latitude daytime (x = 60°) I'-region 
electron densities, corresponding to solar maximum and solar minimum 
periods. 

In order to study the diurnal variation of the electron densities 
in the D region, a set of four profiles has been obtained on July 24, 

1968, at solar-zenith angles 90°, 84°, 18° (noon) and 60° (Figure 1.3). 
As mentioned before, the present study was undertaken with a view to 
interpret the diurnal variation of the O-region ionization, as ex- 
hibited by these electron-density profiles, in particular those cor- 
responding to 18° and 60° solar zenith angles. Some of the main 
features of these two profiles are summarized below. 

(a) Between 90 and 100 km, the two profiles have approximately the 
same gradient, with the values differing by a constant factor 
of about 3. 

(b) Between 85 and 90 km, the profiles have sharp gradients, particu- 
larly the noon profile whose number density increases from about 

3-3 4-3 

1.8 x 10 cm to about 2.1 x 10 cm , i.e. more than an order 
of magnitude. At 85 km, the difference is only a factor of 1.5. 

(c) The profiles have a slow rise between 75 and 85 km> with the 

difference increasing to a factor nearly 3 at 75 km. The noon 

3-3 3 

profile values increase from 1.2 x 10 cm to about 1.8 x 10 

_ 3 

cm in this height interval. 

(d) The difference between the two profiles decreases below 75 km, 
and around 66 km, they cross over. Below this level, the noon 
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profile has a very sharp gradient with the values dropping 

-3 -3 

from about 150 cm at 66 km to nearly 15 cm at 65 km. 

(e) Below 65 km, the drop in the afternoon profile is rather gradual, 
while the noon profile exhibits a slow variation. This dif- 
ference in the height variations makes them to cross again just 
below 60 km. 

The X-ray measurements taken by the SOLRAD 9 satellite on this day 
did not show any signs of solar flares (Figure 4.6) while the 10.7 cm flux 
recorded on this day was 148 units. This indicates that the measure- 
ments were taken on an undisturbed day during a high solar activity 
period. 

5.2 Toni c-Reaetion Schemes 

5.2.1 Positive-io n reactions. Prior to detection of the hydrated 
cluster ions in the lower D region, the entire D region was thought to 
comprise only 0 2 * and NO* ions [Nieolet and Aikin 3 I960]. As de- 
scribed in Chapter 4, the X-rays and EUV radiation below 1027 A ionize : 

0 2 while X-rays and EUV radiation below 911 A as well as energetic 
particles ionize both 0 2 and N 2 - N 2 + thus formed quickly undergoes a 
charge transfer reaction to yield 0 2 + . Therefore, the abundance of N 2 + 

in the D region during undisturbed periods becomes negligible. The 

+ + ' 

0 2 formed in this manner yields NO through direct the charge transfer i 

reaction 

O/ + NO •+■ NO* + 0 2 ! 

and also through the charge rearrangement reaction I 

0 2 + + N + NO* + 0 | 

* 

| 
I 
5 

1 

f! 
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Of course, NO* is also produced directly through Ly-a ionization 

•Ja -fi 

of NO, In addition, NO is produced from N 2 through the charge re- 
arrangement reaction 

N + + 0 ->■ NO* + N 

and through the charge transfer reaction 

N 2 + + NO + NO* + N 2 

The first of these two reactions is not so important as far as NO* 

production is concerned but it plays a significant role in determining 

4 

the yield of atomic nitrogen, both the ground state N[ S) and the 
2 

metastable N( D) species. 

Above 100 km the production of O' is also important. This ion 
yields NO* and 0 2 * through reactions with N 2 and 0 2 , respectively. 

Unlike in the case of the molecular ions, 0 does not undergo recombina- 
tion with electrons. . Its loss rate is determined solely by these ion- 
neutral reactions. 

This scheme of positive- ion reactions, however, does not take 
into account the presence of the hydrated cluster ions in the mid and 
lower portions of the D region. The present reaction schemes which 
include the cluster ions are based on the scheme suggested by Fehsenfeld 
and Ferguson [1969]. 

According to this scheme, the clustering reaction chain starts with 
0 2 * as the precursor ion with the 3-body association reaction 

0~* + 0- + M + 0 .* + M 
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Though it appeared to be successful initially, it failed when the 0 2 + 
production rates due to 0,( 1 A ) ionization were reduced as shown by 
Huffman et at. [1971], It was later shown that the yielj of the hydrated 
cluster ion 0 2 »H 2 0 through the reaction 

°4 + + H 2° + °2 + * H 2° + °2 
is further reduced due to the reaction 

°4 + + 0 + ° 2 + + °3 * 

which reconverts 0^ into 0 2 , 

The rate constant for this reaction measured and reported by 

Fehsenfetd and Ferguson [1972b] indicates that if the [0]/ [l^O] ratio is 

sufficiently high, it can arrest the formation of the hydrated cluster 

ions. The introduction of this intercepting or short-circuiting reaction 

leads to the hypothesis that the sudden disappearance of the cluster ion 

population above the 82-85 km level could be caused by the sharp increase 
3 

in the 0( P) concentration observed in this height range [Ferguson, 1971]. 

Following the work of Lineberger and Puokett [1969], Ferguson 
[197X3 suggested a new series of reactions where the water molecules 
clustered to N0 + through direct hydration 

N0 + ‘CH 2 0] n + H 2 0 + M + N+ *( H 2°)„ +1 + M 

« = 0 , 1 , 2 , 

This direct hydration scheme, however, does not convert the NO + to its 
hydrate fast enough. An alternate path, faster than this, was found 
by Dunkin et at. [1971]. According to them, the association reaction 
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N0 + + C0 2 + M <*■ N0 + *C0 2 + M 
is followed by the fast switching reaction 

no + *co 2 + h 2 o + no + *h 2 o + co 2 

Drnkin et at. [1971] demonstrated the possibility of a similar 
association reaction involving N 2 , but he could measure only an upper 
limit for the rate constant which was found to be very low. However, 
Eeimevl et at. [1972] later established that its rate constant is much 
higher, and that it could form a faster path for the production of 
NO «C0 2 , according to 

N0 + + N 2 + M -s- NO + »N 2 + M 

NO + -N 2 + C0 2 -* N0 + *C0 2 + N 2 

The N0 + *C0 2 thus formed can undergo hydration to form N0 +, H ? 0 
shown above. 

-f 

The reactions forming the multiple clusters of NO *H 2 0 beginning 
from this ion were given by Linebergen and. Puoket [1969]. However, 
according to Heimerl [Thomas, 1974] the association of NO *H 2 0 
with N 2 and subsequent switching with H 2 0 could form a more efficient 
mechanism to produce the multiple clusters of N0 + *H 2 0. The reaction 
sequence here is 

no + *h 2 o + n 2 + m -*• no + -h 2 o*n 2 + M 

% 

N0 + -H 2 0 + N 2 + H ? 0 N0 + *CH 2 0) 2 + N 2 . 
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Unfortunately, no reaction rate data are available for these 
reactions, and hence are not included in this study. 

According to Lineberger and Puckett [1969] , the multiple clusters 

4 * 

NO • (H 0 Q) with n >_ 3, can get converted to H (H„0) through the binary 
reactions 

NO + *(H 2 0) w + H 2 0 H + -(H 2 0) n + HN0 2 

n >_ 3 

This conversion reaction does not take place for n < 3 in view 

-J* 

of reaction energetics. Therefore, the hydration of NO ultimately 
leads to the products 

no + + no + *co 2 + no + *h 2 o -»■ no + *ch 2 o) 2 -*■ h + *ch 2 o) 3 

*4 

and the higher hydronium clusters. The absence of both H^O and 
H • (H 2 0 ) 2 in this sequence creates a problem, because these are the 
two observed ions which dominate the middle portion of the D region where 
the primary ionization source is NO. 

To avoid this difficulty, Burke [1970] suggested the reactions 

N0 + ’H 2 0 + H ^ H + *H 2 0 + NO 
N0 + -(H 2 0) 2 + H + H + -(H 2 0) 2 + NO . 

Laboratory measurements of the rate constant of the first of these 

-11 3 -1 

two reactions gave values less than 1 x 10 cm sec [ Fergusem , 1971] 

which is too slow to make this path acceptable. 
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It was shown recently by Heimerl et at. [1972] that a fast path 
for the conversion of NO into the corresponding hydroniura ions 

could be achieved thorugh the reactions. 

no + -h 2 o + OH h + -h 2 o + no 2 
no + *ho 2 + ho 2 -> h + *h 2 o + no 3 . 

They showed that when established values of OH and H0 2 are used, 

-9 

the rate constants of these two reactions have to be about 5 x 10 
3 -1 

cm sec in order to become effective in this conversion process. 

In this study the value adopted for the rate constants of both reactions 
is 1 x 10 -9 cm 3 sec”*. Once H + *H 2 0 is formed Heimerl et at. [1972] 
suggested that its multiple hydrates are formed more rapidly through 
association with N 2 and subsequent switching with H 2 0 rather than through 
direct clustering, as mentioned earlier. With assumed data on reaction 
rates, this reaction path has been included in our study. 

The complete positive-ion reaction scheme adopted in this study, 
along with the reaction rate constants used is given in Table 5.1. The 
more important reaction paths are shown in Figure 5.1. 

5.2.2 Negative-icm reactions. The three-body attachment of a 
free electron to a neutral oxygen molecule produces a negative ion. 

The 0 2 ~ ion so formed has a relatively low electron affinity. As a 
result a series of charge transfer and charge rearrangement reactions 
take place giving rise to a complex reaction scheme. The various nega- 
tive ions thus formed include 0^”, 0^”, CO^”, C0^”, NC> 2 , N0 S and a 
large number of hydrated cluster ions. Associative detachment to 
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Rate Constant 

n -1* 
cm sec 


1.0 (-12) (300/T) ** 

Fergus on et at. [1965] 

2,2(-lJ) (300/T) 

Ferguson et at. [1969] 

1.4C-10) 

Fehsenfeld et at. [1970] 


4. 7 C- 11) (300/T) 

Ferguson et at. [1969] 

3. 3(-i0) 

Fehsenfeld et at. [1970] 

1 . 8 (- 10 ) 

Ferguson [1967] 


6. 3(-10) 

Fehsenfeld et at. [1970] 

2. 4 (-30) (300/T) 2 
Good et al. [1970b] 

2.6(-28) (300/T) 2 
Fehsenfeld et al. [1971a] 

2. 2 (-9) 

Fehsenfeld et at. [1971a] 


3.0C-10) 

Fehsenfeld and Ferguson 
[1972b] 

1.9(-9) 


Fehsenfeld et al. [1971a] 




reactions . 
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(Table 5.1 continued) 

R iis 0 2 + *h 2 o + h 2 o — ^ h 3 o + OH + o 2 


4 14 H 3 0 + .0H + H 2 0 — H 3 0 + .H 2 0 + OH 


R 115 H 3 0+ + H 2° + M — *■ H 3 0+ * H 2° + M 


r H6 H 3 0 + + N 2 + M — > H/.N, + M 


R 117 H 3 0+ * N 2 + H 2 0 -^ H 3 0+ ' H 2 ,:, + N 2 


R 


R 


R 


R 


llg H 3°* H 2° + H 2 0 + M — > H 3 0*(H 2 0) 2 + M 

119 H 3 0t t H 2 0 ) 2 + H 2° + M — " H 3° tCH 2 0 ^3 + M 

120 H 3 0t f H 2 0) 3 + H 2 0 - M — > H 3 0t(H 2 0) 4 * M 

121 h 3 °^ ( H 2°) 2 + M — 4- H 3 0tH 2 0 + H 2 0 + M 


0.3(-9) 

Fehsenfeld et al. [1971a] 
3.2(-9) 

Fehsenfeld et al. [1971a] 

3. 46 (-27) (300/T) 2 

Good et at* [1970a, 1970b] 

1.4(-30) (300/T) 2 
Heimerl et al. [1972] 

1 . 0(-9) 

Heimerl et al. [1972] 

2. 24 (-27) (300/T) 2 

Good et al. [1970a, 1970b] 

2. 32 (-29) (300/T) 2 

Good et al, [1970a, 1970b] 

0.9 (-27) (300/T) 2 
Good et al. [1970a] 

7,0(-26) 

Good et al, [1970b] 


R 122 

h 3 o*(h 2 o) 3 + m — >■ h 3 o^(h 2 o) 2 + h 2 o + m 

7.0(-18) 

Good et al, [1970b] 

R 123 

NO* + N 2 + M — NO* N 2 + m 

3.5(-31) (300/T) 2 
Heimerl et al. [1972] 

R 124 

N0 + + C0 2 + M — >- N0 + *C0 2 + M 

2. 4(-29) (200/T) 2 
Dunkin et al. [1971] 

R 125 

N0 + + H 2 0 + M *■ NO + *H 2 0 + M 

1.29 (-28) (300/T) 2 
Howard et al. [1971] 

R 126 

no + *n 2 + cq 2 — 4 - no + *co 2 + n 2 

1.0 (-9) 

Heimerl et al. [1972] 

R 127 

no + *n 2 + h 2 o — v no + *h 2 o + n 2 

l.Q(-9) 

Heimerl et al. [1372] 

R 128 

N0 + *C0 2 + H 2 0 — 4 . N0 + *H 2 0 + C0 2 

1.0(-9) 

Dunkin et al. [1971] 

R 129 

N0 + *H 2 0 + H 2 0 + M — 4 - N0 + *(H 2 0) 2 + M 

1.12 (-27) (300/T) 2 
Howard et al. [1971] 
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(Table 5.1 continued) 


130 

N0*(H O) 0 

* 

H-0 + M — * NO * (H o 0) _ + 

M 

1.30 (-27) (300/T) 2 

4 4 


4 do 


Howard ct al. [1971] 

131 

NO** (ILO) , 

+ 

M — ► N0*H-0 + H-0 + M 


1.5 (-14) 

4 4 


4 4 


Dunkin et at. [1971] 

132 

NO* (HO) 

+ 

M — NO* (H-0) - + H-0 + 

M 

1 . 3(-12) 

4 O 


4 4 4 


Dunk-in et al . [1971] 

133 

KO * (H-0) » 

+ 

H„0 — > H,0* (H-0) „ + HNO- 

7.0C-11) 

4 O 


4 O Z. Z. 


Howard et al. [1971] 

134 

N0 + *H 2 0 + 

H 

— >- H 3 0 + + NO 


5.0 (-12) 
Estimate 

135 

N0 + *H 2 0 + 

OH — »- H^O* + N0 2 


1.0(-9) 

Heimerl et al. [1972] 

136 

NO* *H-0 + 

H0„ H„0* + NO- 


1.0(-9) 

4 


4 0 O 


Heimerl et al. [1972] 





N> + 



Figure 5.1 More important positive-ion reaction paths in the D region. 
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neutrals, mutual recombination with positive ions and photodetachment 
during the daytime are the main loss processes of these negative ions. 

The development of the D region negative-ion chemistry has been 
based largely on the laboratory measurements of the reaction rate 
constants and model calculations, and not so much on actual ion composi- 
tion measurements. 

The reaction chain from 0 2 takes place in two different paths. 

In one, C> 3 ~ is formed through charge transfer to C> 3 , and in the other, 

0 4 " is formed through a 3-body association reaction with 0 2 « In the 
first path, 0 3 ~ yields C0 3 ~ through a charge rearrangement reaction 
with C0 2 , which in turn yields N0 2 ” through a similar reaction with NO. 

In the other path, charge rearrangement reactions with the same con- 
stituents yield C0 4 " and N0 3 successively. 

While the above sequences gives the main reaction paths, a large 
number of reactions involving atomic oxygen and other minor constituents 
also take place causing the negative-ion reaction scheme to be a complex 
one. The present reaction scheme which includes these molecular negative 
ions was developed by Reid [1970]. Among the more important secondary 
reactions are those involving OC^P). These cause reformation of 0 2 
from 0 3 _ and C0 3 ~, conversion of 0 4 ” to 0 3 and CO^ to C0 3 , respec- 
tively. Hence, the negative-ion distribution depends to a large extent 
on the relative abundance of 0 3 and 0. 

The complete reaction scheme used in this study is shown in 
Table 5.2. 

Among these molecular ions, the high electron affinity of N0 3 

$ 

makes it a terminal ion. The possible mechanisms of its destruction 


Table 5.2 


194 


Rate Constant 

n -1* 
cm sec 


1.4 (-29) (3O0/T)e" 1,2 ^ RT ** 
Phelps [1969] 

1.0(-31) 

Phelps [1969] 

2.5(-10) 

Phelps [1969] 

2 . 0 (- 10 ) 

Phelps [1969] 

3.0(-10) 

Fehsenfeld et al. [1967] 

8.0 (- 10 ) 

Fehsenfeld and Fergus on 
[1968] 

l.O(-lO) 

Lelevier and Brans comb 
[1968] 

4.0(-10) 

Fehsenfeld et al. [1967] 
l.O(-ll) 

Fehsenfeld et al . [1967] 
8.0(-ll) 

Fehsenfeld et al. [1967] 
9.0(-12) 

Fehsenfeld et al. [167] 

8.0(-ll) 

Ferguson [1969] 


1 . 8 (- 11 ) 

Fehsenfeld and Ferguson 
[1968] , 
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[Table 5.2 continued) 


R N0 2 ' + H — *■ OH - + NO 3.0 (-10) 

Fehsenfeld and Ferguson 
[1972a] 


R 165 

0 2 - + 0 2 + M — *■ 0 4 - + M 

1.0C-29) 

Estimate 


R 166 

° 4 + 0 *" °2 + °3 

4.0C-10) 

Fehsenfeld et al. 

[1969] 

R 167 

o 4 “ + co 2 — s- co 4 " + o 2 

4.3C-10) 

Fehsenfeld et al. 

[1969] 

R 168 

0 4 _ + NO — *■ N0 3 ” + 0 2 

2.5(-10) 

Fehsenfeld et al. 

[1969] 


CO “ + 0 — CO„" + 0- 

1.5C-10) 



Fehsenfeld et at. [1969] 

R 170 C0 4 ‘ + NO > N0 3 ” + C0 2 


4.8(-ll) 

Fehsenfeld et at. [1969] 
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are the mutual neutralization through reactions with positive ions, and 
photodetachment during daytime. This latter process yields free electrons. 
Therefore, the relative abundance of the D region electrons and negative 
ions during the daytime is governed to a large extent by the photodetach- 
ment of the terminal ions [Thomas et at. * 1973], 

Recent laboratory work and daytime negative -ion composition measure- 
ments have indicated that the D region negative molecular ions can undergo 
hydration reactions forming water cluster ions such as 0^ 

CC> 2 ~*H 7 0 , and NO^-I^O, and their multiple hydrates. 

Fortunately, the rate constants for the non-terminal ions indicate 
that the hydration processes are slower than the charge transfer and 
charge rearrangement reactions mentioned above. Even for the terminal 
ion N0 3 ~, the losses due to hydration become important only if the 
photodetachment rate is less than about 10 sec [Thomas et al. t 
1973]. 

In the event the photodetachment rate of NO^" is indeed less than 
-3 »1 

10 sec , the negative terminal ion will be either NO^ ‘I^O 
or one of its multiple clusters. The abundance of this terminal ion 
will then be determined by its recombination rate with the positive 
ions and its photodetachment rate. None of these quantities is presently 
known. If one is to include these ions in the reaction scheme, appro- 
priate values for these quantities have to be assumed. Since this is 
true for NO^ - also, it makes little difference whether the terminal 
ion used is NO^ - or NOj~ • with n > 1 . 

There is, however, a difference between the two types of ions, 
which could affect the electron -density distribution. That is the 
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mutual neutralization coefficient with the positive ions. Considering 
the fact t v at the positive cluster ions have recombination coefficients 
of the positive molecular ions and that one could expect the mutual 
neutralization rates of positive cluster ions with either the negative 
molecular ions or the negative cluster ions, the hydrated cluster ions 
are totally excluded from the present negative-ion reaction scheme. 

The more important reaction paths in the adopted reaction scheme 
are shown in Figure 5.2. 

5.3 Reaction Rate Constants 

The accuracy of model calculations depends largely on the availa- 
bility of correct reaction rate constants. Unfortunately, the measure- 
ments of rate constants in laboratories are subject to various sources 
of error, and the best measurements are accurate only to within ± 30%. 
This is true for the main products, and for secondaries the errors are 
much larger, sometimes as much as a factor of 2. 

The rate constants of the molecular ion reactions shown in 

Table 5.1 are those measured by Ferguson et at. [1969] and Fehsenfeld 

et at. [1970], The accuracy of most of these values have been estimated 

2 4 

to be ± 3 The only uncertainty is the ratio of the N( D)/N( S) 
yield in the reaction between and 0, The values used in this study 
are 0.75. 

The rate constants for the reactions leading to H (f^CO^ ^ rom 
are taken from the results of Fehsenfetd et at. [1971a] and Good et at. 
[1970a, b]. These values are in good agreement with the original values 
measured by Fehsenfeld and Ferguson [1969]. 
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The reactions leading to the production o£ N0 +, CH 2 0)^ with n = 1, 2, 
3, and the subsequent conversion to H + (H 2 0) 3 were investigated by 
Fehsenfetd et at, [1971b] and Howard et at, [1971] . The measurement 
of the relevant rate constants was carried out by both groups at 295 K. 
For the third body, both groups used He, Ar, and H 2> but the latter 
group used 0 2 in addition. Considering the possible errors involved, 
the agreement between the corresponding values obtained by both groups 
is very good. The values used in this study are those of Howard et at, 
[1971]. 

The rate constants of the 3-body reactions generally vary with the 
type of the third body. A close examination of the results of the above 
workers shows that is more efficient than 0 2 as a third body, the 
ratio of the corresponding rate constants being 3:2 [Fehsenfetd et at, 3 
1971a; Howard et at. s 1971]. It is also noted that the efficiency of He 
as a third body is rather poor, while Ar gives rate constants which are 
close to corresponding values obtained with 0 2 - Measurements with NO 
show that its efficiency is similar to that of N 2 in the 3-body reactions 
[Howard et at, a 1971]. In view of the wide difference in the efficien- 
cies of N 2 and 0 2 as the third body, a weighted average was taken for 
tabulating the values in Table 5.1, wherever values with either N 2 or 
0 2 are available. In other instances, the available values were adopted 
without any correction for the third-body efficiency. 

The rate constants of the 3-body reactions are also known to depend 
on the temperature, with the values increasing with decrease of the 
temperature. Ferguson [1971] estimated the rate constant of the 0 2 
to ^2 + *®2 c ^ uster ^ n 8 reaction to increase from 2.8 x 10 cm sec 
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2Q 3 -1 

to 1 x 10 cm sec between 300 K and D region temperatures. 

Dunk-in et at. [1971] found that the rate constant of the N0 + to 
N0 + *C02 clustering reaction to increase by a factor of 2.5 between 
the temperatures 300 K and 200 K. On the basis of these observations 
the temperature variation of the 3-body reactions may be expressed as 
(T/300) -n where n > 2. However, taking a more conservative figure, 
this expression is taken as (T/300) ^ for the rates of 3-body clustering 
reactions . 

One of the most important parameters that determine the electron 
densities of the ionosphere is the recombination coefficient, a, of the 
positive ions. In the region where the negative-ion concentrations are 
small, the electron losses take place solely due to recombination reac- 
tions with the positive ions. 

Laboratory measurements as we 1 1 as ionospheric observations have 
indicated the wide difference that exists between recombination coef- 
ficients of the molecular ions and those of the cluster ions. Among 
the molecular ions N 2 + has a recombination coefficient of 2.6 x 10 7 

7 — 1 

cm sec - which is temperature independent. The other two molecular 
ions, on the other hand, have temperature dependent recombination coef- 
ficients. Measurements made by Kaener and Biondi [1968] have shown a 
1/T temperature dependence for a (0 2 ) , with the values given by 

a(0 2 + ] = 2.2 x 10" 7 (300/T) cm 3 sec -1 . 

The recombination coefficient of NO*, cc(N0 + ), has a 1/T -1 * 5 temperature 
dependence, and the coefficients are given by [Weller and Biondi , 1968] 
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a(NO + ) = 4.1 x 10“ 7 (300/T) 1 * 5 cm 3 sec 


The accuracies of these laboratory measurements are quite good, 

being of the order of ± 10%. They are also in good agreement with the 

results of theoretical calculations which for 0 2 yields a value of 

2.8 x 10 -7 cm 3 sec -1 at 300 K [Chan, 1968]. 

Laboratory measurements of the recombination coefficients of the 

hydrated cluster ions have been obtained only recently. The values by 

Leu et at, [1973] indicated that a increased from about 1 x 10 ^ to 

about 10 x 10 cm 0 sec - , corresponding to an increase of the cluster 

ion mass number from 19 to 109. They also concluded that it was un- 

-5 3 1 

likely that a could have a value greater than 1 x 10 cm sec" even 

for heavier clusters. In normalizing the data which were obtained at 

different temperatures to a standard temperature of 300 K, they have 

- 1/2 

assumed a temperature dependence of T . 

For many other cluster ions, no direct measurements of a generally 
increased with the size of ion. Assuming this relation to hold true for 
other ions as well, values of a were estimated on the basis of the 
number of clustered molecules. These, together with the values obtained 
for other ions, are shown in Table 5.3 with the appropriate temperature 
dependence. 

The value for the mutual neutralization coefficient used is an 

average value based on the measurements done on different combinations 

of positive and negative ions. Hirsh and Eisner [1972] obtained 
-7 -1 3 -1 

1.7 x 10 and 0.34 x 10 cm sec for this quantity in the reactions 


Reaction 


Table 5. 
Recombination cc 


N 2 + a — + N + N 

0 2 + + a — + O + O 

NO + 0 —v N + 0 

O4 + s + 0 2 + 0 2 

jj 

i + 

I 0 2 *H 2 0 + e — + 0 2 + H 2 0 

I 

H^O *0H + e — *- Neutrals 

1 

I + 

j H^O + e — *■ Neutrals 

1 

I 

[ 

\ H 3 0 + *H 2 0 + e — + Neutrals 

Hj0 + *(H 2 0) 2 + e — + Neutrals 

Hj0 + * (H 2 0) 3 + — + Neutrals 
f + 

f H„0 *N_ + e — *• Neutrals 

fj 3 2 

I 

k N0 + *N 2 + e — + Neutrals 

p: 

r, , ■ j 

N0 + *C0 2 + e — + Neutrals 
4 N0 + *H 2 0 + e — > Neutrals 

* ~Yl 

Read a(-n) as a x 10 . 
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(Table 5.3 continued) 


NO * (H 2 O) j + e — *• Neutrals 


5.0 (-6) (300/T) 
Estimate 


NO • (H 2 O) j + e — -*■ Neutrals 


6,0 (-6) (300/T) 
Estimate 


X + Y 


Neutrals 


2.0(-7) 

Estimate 


II. 
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between NO* and NO.,", and NO* and N0 3 ~, respectively. However, in the 
region where the effect of the negative ions becomes important, the 
dominant positive ions arc the hydrated clusters of H^O . Judging by 
their recombination coefficients, one could assume that the mutual 
neutralization coefficients of these ions are higler than those cor- 
responding to the molecular ions. Hence, a conservative value of 
-7 3 -1 

2 x 10 cm sec is adopted for this coefficient. 

The rate constants given in Table 5.2 for the negative-ion reactions 
are those reported by Fehsenfeld et at. [1967] , Fehsenfeld and Ferguson 
[1968], Phelps [1969] and Fehsenfeld et at. [1969], Most of these are 
fast binary reactions which are independent of the temperature. Excep- 
tions are the three-body attachment reaction between 0 2 and free 
electrons, and the three-body association reaction between and 
Og, which yields The former has a temperature dependence 

given by [Phelps* 1969], and its reaction rate is also 

known. On the other hand, the rate constant of the latter is not so 

7 *] £ 

well known. Paok and Phelps [1970] obtained 4 x 10" cm sec" 
at 300 K for this quantity. According to Fehsenfeld et at. [1969], 

-30 

however, the rate constant for this reaction is likely to exceed 10 
cm sec" at Z>-region temperatures. Thomas et al. [1973], in their 
model calculations, used a value of 2 x 10“^ cm^ sec * for this rate 
constant. 

Rate constants for many of the clustering reactions of negative ions 
with water vapor have been reported in the literature during the past few 
years. Since these hydrated ions are not included in the negative ion 
reaction scheme adopted, values of these are not shown in Table 5,2. 


The photodetachment rates used in the present computations are 
given in Table 5.4. The values corresponding to and 0^” are those 
measured experimentally and reported in the DASA Reaction Rate Handbook. 
The values used for NOj” and NO^"* are those used by Thomas et at, 

[1973], In view of the many uncertainties in both the negative- ion 
chemistry and the experimental data on the daytime negative- ion composi- 
tion, no effort was made to determine the effect of the photo detachment 
rates on the negative-ion composition and the electron density in the 
lower D region. 
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Table 5.4 

Photo detachment rates. 


Reaction 

Photodetachment 

-1 

sec 

0 ^ + hv — ► 0 2 + e 

0.33 

0 S " + h\> — *■ 0 3 + e 

0.66 

C0 3 " + hv — C0 3 + e 

0.04 

no 2 ~ + h\> — *■ no 2 + s 

0.04 

NOg + hv * NOg + e 

0.04 


Reference: Thomas et al . [1973] 
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6. ION COMPOSITION AND ELECTRON-DENSITY PROFILES 

Several photochemical calculations have been carried out in the 
recent past in order to interpret the behavior and the structure of 
the D region as revealed in ground-based and rocket measurements. 

Among those who investigated the daytime D-region positive ion and 
electron densities are Goldberg and Aikin [1971], Hunt [1971a, 1973], 
Keneshea and Swider [1972] and Roue et al .. [1974] . Goldberg and Atkin 
[1971] made use of neutral constituent concentrations, adopted from 
values obtained by others, to determine the distribution of individual 
positive-ion species. The model calculations carried out by Keneshea 
and Swider [1972] and Hunt [1973], on the other hand, are more compre- 
hensive. They calculate the relevant neutral constituent concentrations, 
positive-ion densities, and the negative-ion densities in the same 
computer code. Rowe et al. [1974] developed a simplified model of the 
D region ion composition, subject to several constraints, to study the 
behavior of this region under different conditions. 

In this chapter, a numerical model for the D region ion composi- 
tion is developed concurrently with the minor neutral constituent model 
described in Chapter 3. In the next section, details of the numerical 
method for calculating the height distribution of several ionic species 
are presented. This is followed by a discussion of the computed ion 
composition and electron-density profiles. 

6.1 Ion-Composition Model 

The ion-composition model developed in our study incorporates 17 
positive ions and 7 negative ions. The difference between the densities 
of positive ions and negative ions gives the electron density. 
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The positive ions included are: 
molecular ions: N 2 , 0 2 and NO , and 

cluster ions: ° 2 + *°2» °2 + ’ H 2 0, H 3°* 0H ’ H * H 2 °> H*(H 2 0) 2 » 

Ht(H 2 0) 3 , H*CH 2 0) 4 , H 3 0tN 2> NO^N 2 , 

N0^C0 2 , N0*H 2 0, NO‘ I: CH 2 0) 2 and N0*(H 2 0) 3 

The metallic ions present between 80 and 90 km are not included 
in view of their negligible abundance under normal conditions. 

The negative ions included in the study are the following: 

molecular ions: ° 2 » °3 J C0 3 » N0 2 » N0 3 and 

cluster ions: 0 2 *®2 an< * ^^2 *°2 

For reasons mentioned before, the negative hydrated cluster ions 
are not incorporated into the model. 

6.1.1 Continuity equations. The set of continuity equations 
for the charged species can be written in the form 

^p-=Q-P[X] ( 6 . 1 ) 


where 

Q is the production term for a given ionic species, X, and 
P is the loss coefficient. 

The chemical time constants of the reactions involved are generally 
very small. As such the transport terms are not included in the con- 
tinuity equations for the ionic species. 
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In general, the production term for a positive ion X^ can be 


written in the form 


5 - *m *1 kA lX ■ 


(. 6 . 2 ) 


where 


q m is the primary ionization rate, and 


1 k A [X + ] is the total contribution due to charge transfer or 


rearrangement reactions producing X from other 


positive ions , 


In equation (6.2), q appears only for the molecular ions 0 


N/, 0 2 + , and NO*. 


The loss coefficient for a positive ion generally appears as 


P(0 = I k B + «»l>] + ] 


(6.3) 


where 


£ k B is the total contribution due to reactions converting 


X m into other positive ions, 


a [e] is the term due to dissociative recombination with 
m 


electrons and. 


is the mutual neturalization coefficient, and 


is the total negative-ion density. 


In the case of a negative ion X^ the production is generally 


due to charge transfer or rearrangement processes, except for 0^ and 
0~ where they are formed by 3-body attachment processes. Therefore 
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Q(0 2 ~) = 303 

(6.4a) 

e(0 4 ~) = 3* 0 2 “ , 

6.4b) 

Q(X m ~) = l k A[X~] n min 

(6.4c) 


where 

3 is the attachment rate, 

2 k is the total contribution due to charge transfer or 

rearrangement reactions producing . 

The loss coefficient for X " is in general, 

m 

P(J ”) = 1 k B + a,[/] + <5 (6.5) 

Tn It 

where 

\k B is the total contribution due to reactions converting 
X m ~ into other negative ions, 

[N*] is the total positive -ion density, and 
fi is the electron detachment rate. 

The continuity equation for electrons is 

- o e - - m i6.6) 

where 

Q is the total ion-pair production rate and 
[e] is the electron density. 
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It is possible to write a continuity equation for the total posi- 
tive-ion density, in the form 


■}* 

- jfe - ~ = Q a - «[e] [N 4 ] - o. [N - ] [N + ] 


(6.7) 


where 


a is the average electron recombination coefficient for 
positive ions in the D region. 

The charge neutrality condition yields the relation 


[N ] = [IT] + [e] 


( 6 . 8 ) 


which can be written in the form 


[N 4 ] = (1+X)[e) 


(6.9) 


where 


A = [N“]/[e] 


( 6 . 10 ) 


Eliminating [N ] and [N"] in equation (6.7), 


- (a+Xa^) (1+A) fe] 2 


( 6 . 11 ) 


An effective recombination coefficient for the electrons, which takes 
into account v>e presence of the negative ions , is defined by 


“ef'f = ( 1+A ^ (a+Aou) 






( 6 . 12 ) 
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Equation (6.12) becomes 


djVj 

dt 


Q l 


a eff 



(6.13) 


Under steady-state conditions , 


0 . = 
e 


a * „ [e 1 ' 
err L J 


(6.14) 


The effective recombination coefficient is a convenient parameter 
that can be used in calculating the electron density using known values 
of the ion-pair production rates. Alternatively, it can be deduced 
from the simultaneous measurement of the electron density and the pro- 
duction functions. 

The condition of steady-state conditions is generally valid around 
noon up to zenith angles of 60° at least. The time derivative factor 
d[N ]/dt is appreciable only at large zenith angles. However, for the 
sake of completeness, this factor is retained in our calculations. 

6.1.2 Methods of computation. The set of differential equations 
written for each of the ionic species is converted to a set of difference 
equations in the form given in equation (3.2S). The neutral constituent 
concentrations required for substitution in the production terms and loss 
coefficients are taken from concurrent solutions of the netural constituent 
model described in Chapter 3. The calculations involving the ionic con- 
centrations were incorporated into a separate subroutine which is intro- 
duced at sunrise on the final day of the netural chemistry model run. 

Initially, prior to sunrise, the values of [X - ] and [e] are all 
set to zero. The ionic concentrations were obtained at the same time 
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intervals as used in the neutral chemistry program. At each of these 
time steps, the optical depth factor was calculated separately and the 
corresponding flux levels were determined. 

In Chapter 3, results of the neutral chemistry calculations were 

presented for two different models of the eddy diffusion coefficient. 

The high D q model yielded NO profiles having mesospheric values larger 

than those given by the low D model. Preliminary calculations indi- 

cated that the high D model generally gave electron densities that are 

£> 

too high compared to the observed values. Therefore, the detailed cal- 
culations were carried out using the neutral chemistry results obtained 
with the low D q model only. 

At each time step, the solution was obtained by iterating the concen- 
trations until the charge neutrality condition was achieved. Generally, 
this is achieved in a few iterations at altitudes above about 75 km, 
where the electron density is used as the controlling parameter. Below 
this height, however, where the abundance of the negative ions begins 
to grow, a large number of iterations was required before the charge 
neutrality condition was reached. At 65 km and below, computations 
were terminated before this condition was achieved in view of the 
excessive computer time taken. 

The exact values of several parameters that go into the calcula- 
tions are not known. Therefore, it was necessary to carry out several 
runs using different values for such parameters, keeping the rest of 
the conditions unaltered. These parameters include the branching 
ratios of the reactions R 103 and the recombination reaction of N0 + , 
the rate constants for the reactions and and the produc- 

tion rates. 


I 


i 
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In Chapter 4, ion-pair production rates were calculated corre- 
sponding to different input conditions. For X-rays below 3 A, the 
exact flux intensity <f>(3 A) was not known. Therefore, two values of 
$(3 A) were used to calculate the productions functions. These are 

1 x 10" 5 erg cm -2 sec" 1 and 0.5 x 10 -5 erg cm" 2 sec -1 , respectively 
(Figure 4.7]. Also for the precipitating electrons, the exact nature 
of the spectral distribution was not known. Therefore, contributions 
due to two hardness indices, viz. y - 2 and y = 3 were used in the 
calculations (Figure 4.9]. Of these, y = 2 gives the higher pro- 
duction rate, and this was used in conjunction with the high X-ray flux 
(1 x 10 - ^ erg cm" 2 sec" 1 ] in the program. The results obtained with 

this combination are described under 'high 0^ + production'. The other 

-5 -2 -1 

two parameters, y = 3 and <J>(3 A] = 0.5 x 10 erg cm sec were 
used together and the corresponding results are described under 'low 

0 2 production'. 

6.2 Positive-Ion Composition 

The computed positive-ion profiles are shown in Figures 6.1 - 
6.3, corresponding to x = 18° (noon), x = 60° (morning) and x - 60° 
(evening), respectively. Each figure has two sets of curves; one for 
the high 0^ production case (solid lines) and the other for the low 
C> 2 production case (broken lines) . 

6.2.1 Molecular-ion concentrations. Among the molecular ions, 
only 0 2 + and N0 + are shown in the figures. The concentrations of N 2 
and 0 are too low to show here. In all of these three cases, NO 
and 0 + profiles cross over; the noon profiles at 93 km, and the other 
two at 97 km. N0 + dominates below these height levels, and 0 2 + above. 


ALTITUDE (km) 



POSITIVE ION DENSITY (cm -3 ) 


Figure 6.1 The height distributions of positive -ion species. These results correspond 
to low eddy diffusion and summer noon (x = 60°) conditions. 
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results correspond to 
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The relative densities of these two ions appear to depend on the zenith 
angle. 

The behavior of N0 + and 0 2 + is generally described in terms of the 

4 * 4 * ’ 

ratio of their number densities, i.e. [NO ]/[0 2 ^ * This parameter can 
be easily measured because it does not require absolute calibration of 
the equipment. Further, this ratio can be used to calculate the 

* 4 * 

average electron recombination rate in the height range where NO and 

4 * 

O 2 dominate. 

In all three sets of profiles, NO and O 2 dominate above about 
85 km. Below this height, the hydrated cluster ions dominate. In the 
noon profiles, the ratio [NO ] / J ^ as a raax i mum value of 4.4 at 87 km, 
which drops down to unity at 93 km. Above this height the ratio drops 
through unity to 0.3 at 100 km. At x = 60° this ratio is a factor 
about 8 at 88 km, reduces to unity at 97 km, and above this height 
reduces to about 0.2 at 100 km. These results indicate that the point 
of unity ratio gets lifted as the zenith angle increases. Under the 

4 * 4 

same conditions, the difference between the NO and 0 2 profiles widen 
near the mesopause. 

Similar results have been obtained in model calculations done by 
Keneshea and reported by Naroisi [1970]. According to these calcula- 
tions, the [N0 + ]/[0 2 ] ratio has values around 10 at 60° zenith angle 
and at heights between 85 and 90 km. It approaches unity at heights 
near 100 km. For smaller zenith angles, the ratio approaches unity 
at lower heights. These results are similar to the observations made 
by Nc&oisi [1970] . 

The 0* concentration at these altitudes is due to direct ioniza- 
tion of 0 2 by X-rays and EUV radiation, as well as due to conversion of 
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4 * *}* *!• 4 » 

0 and N 2 into 0 2 through reactions R 1Q2 and R 104 * Part of the 0 2 thus 
produced is converted into NO through reaction R 1Q7 » The contribution 
made by reactions Rj 01 > Rjqj, ^105* ant * R 106 convert ^- n S 0+ and N 2 + 

4 * 

directly into NO is not significant below 100 km. Above about 90 km, 
the ionization caused by X-rays and EUV radiation results in the 0 2 
production, and the N0 + production due to Ly-ct ionization of NO is 
rather negligible. The formation of NO in this altitude range is mainly 
due to conversion of this 0 2 through reaction R 107 * Hence, the ratio 

4 * 4 * 

[NO ]/[0 2 ’] represents the relative amounts of NO available in this 
region* 

4 - + 

At altitudes below about 85 km, both NO and 0 2 drop rapidly 
causing the hydrated cluster ions to dominate. The effect of high and 

4 * 

low 0 2 production due to variation in the X-ray flux below 3 A and the 
precipitating electron flux appears in this region. Under noon condi- 
tions, this variation is only a factor of 1.5 - 3,5, the values increasing 
as the height decreases. The decrease of both 0 2 * and N0 + in the T) 

\ 

region is really an exponential variation, being approximately linear 
on the semi- logarithmic scale. The 0 2 + profile follows closely the 0 2 + 
production curve. However, in the case of the NO* profile, the values 

4 * 

drop faster than the decrease in the NO production rate, particularly 
above 70 km. This is caused by the conversion of NO into its hydrated 
cluster ions. 

4 * 4 * 

The NO and 0 2 profiles obtained are generally in agreement with 
the profiles measured by Nccroisi [1970] and Johannessen and Krarikousky 
[1972]. 

In certain observations [Goldberg and Aikin 3 1971], it was noted that 
the o/ profile remained more or less constant at a fixed value below 
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about 75 km, instead of decreasing rapidly. One reason for assuming a 
high 0 2 + production rate, effective in this height range, is to verify 
whether this increase in the 0 2 + concentration is caused by an increase 
in the 0 2 * production rate. The computations show that in spite of an 
increase in the 0 2 + production rate, 0 2 * continues to drop. This indi- 
cates that one has to look for a decrease in the loss coefficient of 
0 + in order to explain such a constant value of 0 2 + . However, it has 
to be noted that 0 2 + is only a minor ion at these altitudes, and there- 
fore, its measurement cannot be considered as very reliable. 

6.2.2 Cluster-ion ooneentratione . The noon profile as well as 

the two x - 60° profiles of the positive-ion composition shows that 

the hydrated cluster ions dominate the region below about 86 km. However, 

contrary to the observations, the major cluster ion according to these 

profiles is N0 + *H 2 0 rather than H*»(H 2 0) 2 or H + *(H 2 0) 3 . The height 

range from about 88-89 km to about 72-74 km is dominated by this NO 

-3 

hydrated cluster ion with a peak value of 860 cm near 83 km. 

* 1 * 

According to the given reaction scheme the formation of NO *H 2 0 
from NO* ions takes place with a time constant of about 13 secs at 80 km. 
This is achieved through the formation pf the intermediate cluster 
N0 + *N 2 and subsequent switching to N0 + *H 2 0. Its fastest loss path 
is due to electron recombination, which has a time constant of about 
300 secs at 80 km. The other loss reactions have time constants much 
greater than this value. However, one possibility for achieving a faster 
loss rate is through the formation of the cluster N0 + *H 2 0*N 2 or 
N0 + *H 2 0*0 2 and subsequent switching with H 2 0 to yield N0 + *(H 2 0) 2 . 

This ion in turn can produce N0 + *(H 2 0) 3 in a similar association - 
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r 

i- 

!: 
r- 
\< 

1: 

j; 

l switching process [Arnold and Krankotisky 3 1974], However, these reactions 

#■' . - 

I are not included in our calculations because of the lack of evidence as to 

l 

V- f 

£ their existence and rate constants of the reactions involved. 

1 

| The two paths suggested by Eeimevl et at . [1972] which involve OH 

and H0„ have time constants about 3000-6000 sec. These high values are 

i due to low concentrations of OH and H0„ obtained in the low Z? model 

! 2 e 

| 

calculations. If the concentrations of these two constituents were suf- 

l ) 

; ficiently high, they would have resulted in the formation of H^O through 

V. + 

| reactions R 135 and R 13 gj which could subsequently yield H 

The reaction involving H atonic also has a time constant greater than 1000 
1 sec, corresponding to a rate constant of 5 x 10~ A cm" sec (Reaction 

\ R 112 5 * 

\ ^ . 3 1 

| The multiple hydrated cluster of NO , i.e. NO * (H 2 0 ) 2 an ^ NO • (^0) 3 

s’ ; 

l have negligible concentrations according to the computations. This is 

!; mainly because of the low production rates for these ions. The faster 

:: 4. . 

i path of producing H -(H 2 0) 2 is given by reactions R 108 , R 110’ R ll2’ and 

| R 114 J proceeds according to 

£ 

? '• > - 

; 0 2 + — * 0 2 + > 0 2 — * 0 2 + *H 2 0 — * H 3 0 + «0H — * H + *(H 2 0) 2 ll 

j: } |j 

l + + ! 

s; The alternate path of producing H *(H 2 0) 2 from 0 2 *H 2 0 is j 

I' * 4 : ; \ 

| through the formation of H^O . However, the conversion of H^O into : 

H + • (H 2 0) 2 through the conventional 3-body reaction is found to have j 

i a time constant of about 600 sec. On the other hand, the association/ ] 

l l 

V ' 

switching reactions involving the major neutral particles can convert J •?] 

j; . ) . : j 

j; HjO into H - (H 2 0) 2 within about 5 secs. ' l\ 


l S 



The third path of producing this cluster ion is through the forma- 
tion of NO + *H 2 0: 

OH, HO H O, M 

NO *H 2 0 ^ H *(H 2 0) v H *(H 2 0) 2 


■ 

Thi path is again limited by the conversion rate of NO *H 2 0 into 
H + *CH 2 0), which is over 1000 sec., as mentioned before. 

The existing reaction schemes, therefore, give time constants which 

•f + 

are too long for converting NO into H A 5 s result, the 

concentration of H *(H 2 0) 2 reaches only a low value. It reaches 

_ 7 

a peak value of 390 Jons cm at 81 km. It remains as the second major 
cluster ion between about 74 and 85 km in the noon profile. This height 
range is little affected by the zenith angle variation. 

Above 88 km, H *(H 2 0) concentration drops rapidly, becoming 

*f + 

less important than NO *C0 2 ana H^O . One reason roz this rapid 

3 

drop is the increase in 0( P) concentration which occurs towards the 
high altitudes. 

Towards the lower heights, particularly below 72.5 km, its abundance 
becomes less than that of H + *(H 2 0) 3 and H + *(H 2 0) 4 . The effect of 
0 2 + production rates become more prominent in the x - 60° profiles than 
in the noon profile. 

The higher clusters H + *(H 2 0) 3 and H + *(H 2 0) 4 dominate below 
about 71 km in the noon profiles. The conversion of H •(H 2 0) 2 into 

+ 4 * 

H • (H 2 0) ^ and H * (H 2 0j ^ becomes more efficient with decrease in 
the altitude due to corresponding increase in the [H 2 0] and neutral 
particle density. Since the electron density decreases rapidly below 
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about 75 km, the loss of these heavy ions due to electron recombina- 
tion also gets reduced. At the same time, the negative-ion concentra- 
tions increase below about 70 km. These ions undergo mutual neutrali- 
zation with the positive ions. Therefore, the negative ions become 
increasingly important in controlling the positive ions at these 
altitudes . 

The next higher hydronium cluster H + * (I^O) ^ appears below 82 km, 

_3 

reaches a peak value of about 200 cm at 69 km, and then drops to 
about 35 at 60 km. It also has two profiles corresponding to high and 
low 0 2 + production rates , between about 63 and 81 km. Its abundance 
becomes equal to that of H • (H 2 0) ? near 72.5 km in the noon profile and 
near 72.5 km in the x = 60° profiles. 

The cluster ion H + *(H 2 0) 5 appears to be the dominant D region ion 

below 72 km. Its variation with 0 2 + production varies in the height 

-3 

range from 72 to 74 km only. Its peak value is 7.2 cm which occurs 

at 65 km. This is the terminal positive ion in the hydronium cluster 

series. The ion density in this region is determined by the ion-pair 

production rate which is due to cosmic rays and hence the profiles do 

not show much variation towards the lower heights. 

Another cluster ion that is present in significant amounts is 

NO *C0 2 . It has a peak value of 260 cm around 88 km, which reduces 
•2. 

to 230 cm at x = 60°. This i°n is an intermediate one in the forma- 
tion of N0 + *H 2 0 from N0 + . At higher altitudes, H 2 0 concentration 
drops, and hence the conversion of this ion into N0 + *H 2 0 becomes less 
efficient which results in increasing the concentration of NO *C0 2> 

The distribution of the heavy cluster ions H' ‘(HgO), and 
H + *CH 2 0) 4 agree well with the measured profiles discussed in 
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Section 5.1.1. These are the ions dominating the lower D region. 

However, as mentioned earlier, a major discrepancy occurs regarding 
the dominant ion in the mid-P region. The observations indicate the 
presence of H + *(H 2 0) 2 and H + *(H 2 0) 3 in large quantities in this 
height interval, while N0 + *H 2 0 has been found -;o be only a minor ion. 

The computed distributions, however, show the opposite. N0 + *H 2 0 is 
the dominant ion while H + *(H 2 0) with n - 1, 2, 3 are only minor ions 
in the height range 72-85 km. This demonstrates a serious weakness in 
the existing theories concerning the P region positive-ion chemistry. 

Some clues as to a possible circumvention of this situation have been 
suggested by Arnold et al. [1974]. According to them the rapid forma- 
tion of heavy hydrated clusters could be made possible through (N 2 +0 2 ) 
and C0 2 association/switching reactions. The inclusion of such a 
reaction scheme into the current computer codes is not however feasible 
due to the absence of data on the rate constants of these reactions. 

6.3 Negative-Ion Composition 

The ion composition model yields the density profiles of the nega- 
tive ions 0 2 ", 0 3 ", C0 3 ", N0 2 ", NC> 3 _ , 0^" and CO^ - . All of these ions 

-3 

except 0 3 reach density levels greater than 0.1 ions cm during the 
daytime. Figures 6.4 - 6.6 show the density profiles of these ions 
corresponding to x = 18° (noon) , x - 60° (morning) and x = 60° 

(evening), respectively. The difference in the values obtained for 
the two cases, high f?(0 2 ) and low q(0 2 ), is insignificant for the 
negative-ion densities. Therefore, only the results corresponding to the 

.j. 

high 0 2 production case are given here. Also shown in these figures 
are the corresponding electron- density profiles. 


ALTITUDE (km) 




ALTITUDE (km) 
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NEGATIVE ION DENSITY (cm** 3 ) 

Figure 6.5 The computed height distributions of positive -ion species. These results 
correspond to low eddy diffusion and summer morning (x = 60°) conditions. 
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As mentioned before, the negative-ion model adopted in our study 
is a very approximate one. In reality, it is believed that these ions 
exist clustered to water molecules or possibly other neutral molecules 
as well. However, for the purpose of computing the electron-density 
profiles below 75 km, the present model seems to be adequate. 

In the noon profile (Figure 6.4), the total negative-ion density 
and the electron density become equal at 70 km, i.e. X = 1. Above 
70 km, the negative-ion density drops rapidly so that at 75 km, it 
reduces to 1/27, and at 80 km, it is even below 1/600. The negative- 
ion density is an order of magnitude less than the electron density at 
74 km. Therefore, the negative ions could be completely neglected in 
D region studies above this height. 

Below 70 km, however, the negative ions dominate over the electrons. 
Between 65 and 66 km, the negative-ion density reaches values 10 times 
the electron density. At 63 km, X is nearly 100, and at 60 km, it is 
about 100. 

In the other two negative- ion profiles (Figures 6.5 and 6.6) X 
behaves somewhat similarly. In the morning profile, X = 1 just below 
70 km, while in the evening profile this is close to 72 km. X becomes 
1/10 and 10 near 74 km and 66 km, respectively, in the morning profile. In 
the evening profile this happens near "5 and 66 km, respectively. There- 
fore, the overall behavior of the negative ions relative to the electron 
densities is almost identical in the three profiles. 

The absolute values of the negative-ion density, on the other hand, 
appears to vary with the zenith angle. In the noon profile, X = 1 when 
the negative-ion density is about 360 ions cm , while in the morning and 
evening profiles this occurs when the ion density is 170 and 320 ions 
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cm -3 , respectively. Further, the peak density of negative ions is 700 

7 T 

ions cm at x 5 18°, 360 ions cm” at x = 60° (morning) and 540 ions 
cm at x = 60° (evening) , respectively. 

One should, however, remember that these values are subject to several 
uncertainties. Among these are: firstly, the production function; 

secondly, the photodetachment rates, and finally the mutual recombination 
coefficient of the major negative ions. 

The ion- composition profiles show that N0 3 , CO^ , and CO^ are the 
major negative ions below 75 km. In the noon profile, N0^ dominates 
below 66 km, while C0 3 ” dominates above this height. CO^ is a secondary 
ion with values about a factor 2-3 below the major ion species. The 
next set of ions includes 0 2 “ and N0 2 ~. These have concentrations which 
are about a factor 20 - 30 below the major ion concentrations. Between 
these two ions N0 2 ” dominates below about 67 km, while above this height, 

0 2 " dominates. 

In the negative-ion chemistry, the primary ion produced is 0 2 , and 
the rest are all produced subsequently through charge exchange or rearrange- 
ment reactions. The species that have low loss rates can attain high 
concentration levels. In the case of NO^ , the loss mechanisms are only 
photodetachment and mutual neutralization with the positive ions. Hence 
it is considered a terminal ion, though in reaiity, further reactions to 
form hydrated cluster ions may be possible. Thus N0 3 could attain a 
high concentration. 

In the case of C0 3 ” and CO^”, however, large concentrations are 
possible because of the relatively high abundance of C0 2 present in this 
region. As discussed in Section 5.2.2, the negative-ion reactions proceed 
in two paths, one via 0 3 ” and the other via 0 4 “. Both of these react with 
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COj to yield CO^ and CO^*", respectively. The time constants of these 
reactions are much shorter than those of the other reactions. This results 
in these ions attaining a high abundance in the lower D region. 

The x = 60° profiles also show a similar distribution of the negative 
ions. However, the relative abundance of the different ions appears to 
have changed. In the morning profiles, the absolute concentrations of the 
three major ions have gone down by a factor of about 1.5 - 2 below 70 km, 
while maintaining the relative distribution somewhat identical. Such a 
reduction in the overall negative-ion density can be understood because 
of a similar reduction in the total ion-pair production between noon and 
X = 60°. 

In addition to the photo detachment, collisional detachment processes 

also control the relative abundance of the negative ions and electrons. 

Of such processes, the reaction R. S3 where 0( P) reacts with 0 2 ~ plays 

an important role in controlling the 0 2 level. Between x = 60° in the 

3 

morning and noon, [0( P)] in the height range from 60 to 70 km does not 
change very much (Figure 3.12). However, between noon and x = 60° in the 
afternoon [0( P)] drops through a factor about 1.5 in the same height 
interval. Such a drop in [0( P) ] can result in enhanced negative-ion 
concentrations. This is what is shown in the evening set of profiles. 

3 

In addition to changes in [0( P)], changes in [NO] and fNO^] also contri- 
bute towards the disimilarity in the negative-ion concentrations in the 
morning and afternoon profiles. 

As discussed in Section 5,1.2, no proper measurements of the absolute 
negative- ion densities carried out under quiet daytime conditions are 
available for comparison purposes. Nevertheless, the available 
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measurements indicate the presence of CO^ , NOg and their hydrated clusters 
as the major negative ions, in agreement with the computed results. 

6.4 Electron Density 

The computed electron- density profiles are shorn separately in 
Figure 6.7. These correspond to 18° and 60° zenith angles, the solid lines 
corresponding to the high 0 2 + production and the broken lines corresponding 
to the low 0/ production. The electron-density profiles measured at these 
zenith angles at Wallops Island on July 24, 1968 aTe also shown here for 
comparison. 

6.4.1 Comparison with measured profiles. In discussing the 
positive- ion composition in Section 6.2, it was noted that there is a 
discrepancy between the type of hydrated cluster ion that is dominant in 
the D region as computed and as observed. Fortunately, the dissociative 
recombination coefficient of both these ions has the same value, and 
therefore the electron densities obtained is independent of the actual 
type of the positive ion. As such, the comparison of the computed 
profiles with those measured will still be valid here. 

The general agreement between the two sets of profiles is good. 

However, certain discrepancies do occur. In the notn profile, the com- 
puted values are significantly less than the measured values above 88 km. 
This difference is about a factor of 2.5 above 90 km. The two profiles 
have equal values at 88 km. In this region of the ionosphere, the 
dominant positive ions are 0 2 + and NO*. The metallic ion abundances are 
also negligible. Hence, the effective electron recombination rate can 
easily be determined once a(0 2 + ) and a (N0 + ) are known Therefore, the 
missing factor has to be in the production function. 


LTITUDE (km) 


HIGH q(Oi) 

LOW q (d£) 

MEASURED ( Mechtly and Smith, 1970) 


ELECTRON DENSITY (cm** 3 ) 


Figure 6.7 The computed electron-density profiles corresponding to 18° and 60° solar 
zenith angles. The measured profiles are also shown for comparison. 
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Above 90 km, it is mainly the O 2 molecules that contribute to the 
ion-pair production, which is ionized by X-rays in the 30-40 A band (mainly 
by the 33.7 A line) and Ly-g (1016 A) radiation. Though C> 2 undergoes small 
changes in its density, it is unlikely that the change is large enough to 
cause the observed discrepancy. Therefore, the only other factor left for 
consideration is the incident flux of the 33.7 A and 1215 A radiation 
lines. The values used in this study are those reported by Hinterreger 
[1970], The Ly-g line intensity has been measured many times and its flux 
variation is fairly well known. On the other hand, the 33.7 A line 
intensity is not so well known. The flux intensity used here is that 
measured by Manson [1972]. This is nearly an order of magnitude less than 
the measurements of Argo et al. [1970]. However, measurements of Argo 
et al. cover only a few lines. Manson' s data covers the entire spectrum 
bettveen 30 and 100 A. Naturally, one would prefer to use Manson’ s data. 
However, being on the extreme edge of the spectrum measured, the 33.7 A 
line may have been subject to a greater error than the rest of the spectrum 
in his measurements, which could have resulted in getting a too low value. 
An increase in the 33.7 A intensity can easily increase the electron 
densities between 90 and 95 km considerably. Such an increase can also 
make the x = 60° profile agree better with the measured profile above 
90 km. However, this has to wait until more accurate measurements of this 
line is made available. 

In the height range between 80 and 88 km, the agreement between the 

noon profile and the measured profile is very good. The electron- density 

_ 3 

profile has values between 1300 and 2000 cm in this region. The domi- 
nant positive ions are the hydrated clusters. Though the computations 
show N0 + *H 2 0 to be the dominant ion here, the actual one could be 
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H + *(HJD) , with n 4. The recombination coefficient of such ions 
-5 7 > -1 

are nearly 1x10 cm' sec [Leu et at. , 1973] . This is particularly 
so in view of the valley in the temperature profile which is present in 
this height interval. 

The x - 60° profile, however, shows a significant deviation from 
the measured profile. Between 80 and 88 km, the computed electron densi- 
ties are about 1.3 - 1.4 times greater than the observed values. This 
implies that either the loss coefficient has increased or the production 
function has decreased over the interval of time corresponding to x “ 18° 
and 60°. 

The loss coefficient depend? on the actual domination present at these 

heights, and the effective recombination coefficient. The formation of 

the dominant ion is controlled by the concentrations of several minor 

3 

neutral constituents such as 1^0 and 0( P) . These have long time constants 
around 80 - 85 km, and it is very unlikely that their concentrations would 
vary significantly over this interval of time. 

The effective recombination coefficient depends to a large extent on 
the recombination coefficient of the dominant cluster ions in this region. 
To bring about an increase of a factor of 1.4 in the electron density, the 
effective recombination coefficient ( a e £f) has to increase by a factor of 
2. The value of “ e ££ can change for two reasons. Firstly, it can change 
due to a change in the composition of the positive ions, when the dominant 
ion changes from one having a low recombination coefficient to one having 
a large recombination coefficient. Alternatively, the dominant ion could 
remain unchanged, but its a may change due to a change in the ambient 
temperature. 
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According to Leu et at, [1973], in a simulated environment, the 

dominant hydrated cluster ion could be made to vary from one having a low 

mass to one having a high mass, by changing the ambient temperature. They 

found that between 300 K and 200 K, the dominant ions changed from 

H + »CH 2 0)j and H + *(H 2 0]^ to H + *(H 2 0) 5 and H + *(H 2 0)g, respectively. 

-5 3 -1 

This means that « e ££ changed from about 4 x 10 cm sec to about 
- 5 3 - 1 

8 x 10 cm sec , corresponding to this 100° C drop in temperature. 

On the other hand, the change in a with temperature for a given ion was 

-0 5 

found to be rather small, given by T ’ only. For a temperature change 
from 300 to 200 K, this gives a difference of only 25%, which is much 
smaller than the required factor of 2. 

During the summer days, the mesopause temperatures are considered to 
reach very low values. THe measurements of the mesosphere temperatures 
by Smith et at. [1970] during the summer of 1968 have shown that the 
mesopause temperature could vary generally between 160 and 200 K. They 
also have observed that on a given summer day, the mesopause temperature 
had changed through 20° C within a time interval of about 10 hours. Further, 
near the mesopause, the gradient in the temperature was observed to be 
quite high. Temperature variations of about 20° C within a distance of 
2 km are not uncommon in the mesosphere. Generally, much larger varia- 
tions are observed in the winter months than in the summer months. 

Therefore, it is conceivable that the low temperatures prevailing at 
the mesopause could have given rise to heavier clusters over a narrow 
height range. As discussed earlier, the formation of the heavier clusters 
take place through 3-body association reactions. These reactions have 
rate constants which are strongly temperature dependent. In our study, 
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a T temperature dependence is assumed for these rate constants. In 
reality it is quite possible that the temperature dependence is even 
stronger than this. Thus, it is apparent that the sharp valley in the 
temperature profile around 85 km could give rise to firstly, an enhance- 
ment in the heavier clusters having large values for a, and secondly, an 
increase in a of each cluster ion around the mesopause. The possibility 
of such a scheme was demonstrated in the D region p.-sitive-ion measure- 
ments made recently by Krankcwsky et at. [1972] who observed high con- 
centrations of heavy clusters such as H + *(H 2 0)^ around 85 km on a summer 
night. 

Next, one could consider the possibility of reducing the production 
function within the short period of time concerned. Below 88 km, it is 
the NO ionization due to Ly-a radiation that is most important. Therefore, 
either the [NO] or the intensity of Ly-a has to decrease by a factor of 
nearly 2. However, NO at this altitude has a life time against chemical 
and diffusion losses over a day. Thus, the expected change to occur in a 
few hours appears unlikely. In the case of the Ly-a line, Timothy and 
Timothy [1970] observed only a maximum of 10-12% variation in the inten- 
sity of this line over short term periods, i.e., periods less than a day. 
Hence the necessary 100% variation again seems unlikely. The other means 
by which the production rate could change is through a variation of the 
optical depth. Such a variation could be brought about by changing the 
0 2 density. Below 88 km, however, this again appears to be difficult to 
realize. 

The missing factor of 1.4 in the electron- density profile corresponding 
to x = 60° and 80-88 km range could therefore be attributed most probably 


237 


to the change in caused by a sharp temperature variation at the 

mesopause. 

The most striking feature in the 70-80 km range of the computed pro- 
files is the relatively low dependence of their values on the solar zenith 
angle. Throughout the height interval of 70-75 km, the measured profiles 
maintain a nu.'>b»r density ratio of 3, while the corresponding ratio in the 
computed profiles is about 1.5 - 1.7. 

It is also noted that all four computed profiles, including the profiles 
corresponding to high and low 0 ^ production rates, lie enveloped by the 
two measured profiles; i.e., at x = 18° the measured values are greater 
than the computed values whereas at x - 60° the measured values are 
smaller than the computed values. Therefore, any adjustment that is to be 
made to the computed profiles to bring thorn into agreement with the measure- 
ments has to be selective regarding the solar zenith angle. 

Based on the same arguments presented earlier, any possibility of the 
minor neutral constituents such as NO and H 2 0 varying through the required 
extent may be ruled out. Being away from the mesopause any abrupt change 
in the temperature profile also may be considered unlikely. 

Towards the lower end of this height interval, of course, the presence 

of the negative ions become important. A relative increase in the negative 

ion density makes increase causing a reduction in the electron density. 

Such an increase in N is indeed possible through the collisional 
- 3 

detachment of 0 2 in its reaction with 0( P] . The present computations 
show that 0( 3 P) below 75 km drops markedly from the noon values. This drop 
could get further enhanced if lower values for the eddy diffusion coef- 
ficient are used. However, above 70 km this mechanism is only marginal 
in reducing the electron densities. 


J 
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Referring to Figure 6*7, it is also noted that the difference 
between tne computed and the measured profiles is not so great when the 
lew <7 (Oj ) results are employed. However, with these results, the values 
in the computed noon profile fall short of the measured electron densities 
which is worse with the low ^(0 ^ ) case than with the high q(0^ ) case. 

An explanation for this behavior could be sought in terms of the produc- 
tion mechanisms. 

Below 75 km, Ly-a is still the dominant ionization source. However, 
the X-rays in the band 2 - 10 A and precipitating electrons also make 
significant contributions, particularly towards the lower region. The 
precipitating electron flux has no diurnal variation. The agreement be- 
tween the measured and computed profiles is obtained when low flux (y = 3) 
for this source is used. In fact making it zero would make the agreement 
at x = 60° still better. 

The X-ray wavelengths less than 3 A cause ionization below 75 km. 

An enhancement in the 0-3 R band flux at noon, with its contribution 

exceeding that of Ly-a, can cause a corresponding increase in the noon 

electron density below 75 km. In Chapter 4 when the appropriate fluxes 

were selected from ths available data, 0-(3 A) was taken as 1 x 10 - ^ erg 
- 2-1 + 

cm sec for the high production case. However, Figure 4.6 shows 
that shortly before 1200 hours on the day the measurements were made, a 
sudden enhancement in the 0-3 A flux had indeed taken place. Such an 
increase in the X-ray flux in the 0-3 A band occurring near noon and absent 
at x = 60° could account for the differences appearing below 75 km in the 
computed and measured profiles. In order to obtain quantitative values 
for this increase, however, more detailed recordings of the 0-3 A band flux 
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are required. As discussed in Section 4.2.1, many uncertainties exist in 
both the measurements and interpretation of the solar X-ray fluxes below 
3 ft. Nevertheless, this mechanism provides a satisfactory explanation 
because of the selective manner in which this enhancement affects the 
electron- density profiles, i.e., affecting only the noon values lying 
belovr 75 km. 

The region below 70 km is controlled mostly by the negative ions. 

The overall agreement between the computed electron densities and the 
measured values could be considered satisfactory in view of the many un- 
certainties present in the negative-ion chemistry. However, the computed 
profiles fail to reproduce a characteristic feature revealed in the 
measured profiles. That is the cross-over of the noon and the x = 60° 
profiles near 65 km. The production function does not show such a cross- 
over at these heights. Hence, this could be safely attributed to changes 
taking place in the negative-ion chemistry. As mentioned in the last 
section, many of the parameters such as the photodetachment rates and the 
mutual neutralization coefficients are not properly known. Hence, at 
this stage, attempts to interpret the detailed behavior of the electron- 
density profiles below 70 km would seem futile. 

6.4.2 Solar zenith angle variation. The computed electron density 
profiles do not exhibit a marked solar zenith angle variation, as revealed 
in Figure 6.8. Here, the computed electron densities are plotted against 
the solar zenith angle, for heights of 75, 80 and 85 km. The solar zenith 
angle variations as observed in partial reflection, and cross -modulation 
experiments are also shown in this figure [Throne, 1969]. These show 
larger variations, particularly for 85 km. One could also compare the 
diurnal variation of absorption measurements. However, these are mostly 
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sensitive to heights above 85 km, and Hence not related to the computed 
variations. 

Figure 6.9 shows the effective recombination coefficient deduced 
from the computed electron density values and the production functions. 

As shown in this figure, there is no significant variation in the 

values deduced for x = 18° and x = 60°. | 

The effective recombination coefficient shorn here also exhibits the 
regions where different loss mechanisms occur. Above 88 km, “gff is 

3 _ n + + 

less than 1 x 10 cm sec . The molecular ions, NO and 0 2 domina.te 
this region. Between 70 and 85 km, has values lying in the range 

£ r 1 

3 x 10" - 1 x 10” cni sec” . This region is associated with the hydrated 

cluster ions. In the intermediate region, i.e. 85-88 km, electron loss 
mechanism changes from cluster-ion recombination to molecular-ion 
recombination. 

Below 70 km, increases rapidly. In the 10 km range from 70 km •> 

to 60 km, “ e £f increases by a factor of nearly 1000. This factor would 

undergo a solar zenith angle variation only if the dominant ion in a par- l'\ 

ticular height range varied systematically from one having a low a. to 
another having a high a., or vice versa. In these computations, however, 

- 4 . 

it is found that no such variation occurs. § 

■] i 

In the 75-85 km region, the minor constituents, particularly NO and 1 

- , • 1 

£ j 

H-0, have long time constants and, therefore, their chemistry does not v 

* i :: 3 

reveal a detectable solar zenith angle variation. Hence, in the production :j 

S' K- 

function, the solar zenith angle dependence appears through the variation 

in the incident flux. This takes place due to changes in the sec \ factor j 

' i ■ 

in the optical depth. At higher altitudes, of course, the 0^ density could 

I h 

% ■■ 3 


EFFECTIVE RECOMBINATION COEFFICIENT (cmV) 

Figure 6.9 The effective recombination coefficient of the electrons, obtained from the electron 
densities computed at 18° and 60°. 
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have a diurnal variation which could change the optical depth and hence 
the attenuation of the downcoming radiation. However* the results of the 
neutral chemistry model show that this is not important below 100 km. 

Therefore, it appears that above 70 km, the only systematic solar 
zenith angle variation of the electron densities occurs through changes 
in the sec x term in the optical depth factor. The concentrations of the 
major ionizable constituents such as NO and 0 ^ remain unchanged. Further, 
there is no systematic change in the effective recombination coefficient. 
For heights below 70 km, systematic changes in the negative-ion chemistry 
which take place during the daytime could cause corresponding variations 
in the electron densities. However, in the absence of proper data on the 
negative-ion chemistry, the behavior of the electron densities in this 
region is only a matter of conjecture. 

It has now become apparent that the systematic variations of the 
computed electron densities fail to account completely for the behavior 
of the measured profiles. Hence it was found necessary to invoke addi- 
tional sources that would cause short term variations in the U-region 
electron density. These include the temperature variations around the 
mesopause and enhancements in the solar X-ray flux in the 0-3 A band. 
Variations of such quantities cannot be predicted in advance, and could 
only be detected through continuous monitoring related parameters. There- 
fore, any predication of the solar zenith angle variation of the D region 
ionization will generally fall short of the actual measured variations. 
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7. SUMMARY AND CONCLUSIONS 

An attempt has been made in this study to investigate the solar zenith 
angle variation of the D-region ionization. In particular, the variation 
occurring between 60° zenith angles corresponding to a sununei, quiet-sun day 
during the high solar epoch has been the main interest. One objective in 
this study is to interpret the measured electron-density profiles available 
at two zenith angles 18° and 60°, on July 24, 1968. These two profiles are 
among a series of electron-density profiles measured using in situ rocket 
experiments which are capable of yielding accurate electron densities in 
the D region. 

The ionization parameters studied include the positive-ion density, 
negative-ion density and the electron density. Prior to constructing a 
model for the D region ion chemistry, it was found necessary to construct 
a model for the neutral constituent chemistry encompassing the D-region 
altitudes. In order to minimize the effects of boundary conditions, the 
neutral chemistry model was extended from 30 km to 120 km. 

Model calculations were carried out to study the diurnal behavior of 
14 minor constituents which include oxygen only, oxygen-hydrogen, and 
oxygen-nitrogen species. Initially, the carbon compounds were also included, 
but these were later removed from the model as they were found to contribute 
only a little to the ion chemistry of the D region. As far as the behavior 
of this region is concerned, the most important minor neutral constituents 
are 0( 3 P), 0 3> H 2 0, NO, N0 2 and perhaps OH and HO,. 

The neutral chemistry model includes both the photochemistry of the 
neutral species and the transport terms to account for the eddy diffusion 
of long-lived constituents. In both of these areas many of the parameters 
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used are quite uncertain. The solar-radiation intensity over certain bands, 
the absorption cross section of 0 2 in the S-R band and the rate constants 
of several reactions are some of the lesser known photochemical parameters 
used in the model. In addition to these, the parameters such as the eddy 
diffusion coefficient a 1 the boundary conditions of the constituents 
under diffusion are also subject to much uncertainty. 

Based on the observations of several workers, as well as on the pre- 
vious work of a similar nature, two models for the eddy diffusion coeffi- 
cient were adopted, which are designated high and low D models. Using 
the photochemical production and loss terms, continuity equations were 
written which turned out to be simple first-order differential equations 

for the short-lived constituents, and second-order partial -differential 

3 

equations for the long-lived constituents. These include 0( P) , H, H 2 , 

H 2 0, NO, and H 2 0. 

This set of coupled differential equations was solved numerically by 
computer simulation. The computer solution was carried out in three parts. 
First, the initial values were obtained by solving the steady-state equa- 
tions. Next, the time-dependent equations were solved keeping the solar 
angle fixed at noon. Finally, the diurnal variation was simulated by 
running the program with the solar zenith angle varying through day and 
night. It was found necessary to carry out the diurnal simulation through 
a minimum of ten days to achieve convergence. 

The results of the diurnal study show that most of the neutral con- 

3 

stituents undergo diurnal variation at some height interval. The 0( P) 
concentrations have a diurnal variation throughout the 30-120 km height 
range with the concentrations vanishing below 80 km at night. In 0 2 (*A^) , 
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| the variations during the daytime, appear below 75 km. Corresponding to 

1 

| these variations, 0^ also undergoes significant changes below about 80 km. 

I Among the 0-H species, both H^O and H 2 do not show any diurnal varia- 

tion. On the other hand, H, OH and H0 2 exhibit large diurnal variations. 
The concentration of H vanishes below 77 km in the night, and during the 
daytime, shows a strong diurnal asymmetry, mostly between 50 and 85 km. 
Both OH and H0 2 show somewhat similar variations during the daytime. How- 
ever, at night the depletion of H0 2 is much greater than that of OH. 

In the 0-N species, N0 2 has a strong diurnal variation while NO has 
a diurnal variation only below 70 km. N 2 0 on the other hand has no 
diurnal variation. 

The ion-pair production rates were calculated using the neutral con- 

! stituent concentrations obtained in Chapter 3. The results of these cal- 

! 1 

culations indicate that the ion-pair production due to 0„( A ) is small 
! * 9 

compared to other sources although it is a significant source of 0 2 ions. 
The major contributions were found to be due to Ly-a in the mid-Z) region 
In the typer D region, Ly-6 and 33.7 A radiation lines were found to be 

important, while in the lower region, galactic cosmic rays were found to 

i 

1 be the major source of ionization. 

i 

The diurnal variation in these major production functions was solely 
< due to the sec x term in the optical depth factor. The concentrations of 

L the ionizable constituents NO and 0- remain unaltered. The small varia- 

f! 2 

f: 

i- tion of NO below 70 km does not contribute much because of the rapid 

tv. 

attentuation of Ly-a radiation below 70 km. 

To investigate the variation of the P-region ionization with the 
solar zenith angle, a model for the charged species comprising positive 
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ions, negative ions and the electrons was developed. This model was run 
concurrently with the neutral chemistry model on the last day of its diurnal 
simulation. 

The positive-ion composition results of these computations agreed 
generally with the observed distributions, except for the discrepancy 
observed in the major ion in the mid-U region. The computations show 
NO *^0 to be the major ion whereas direct measurements indicate hydrated 
clusters of the type H + *(H 2 0 ) k with n - 3. The negative-ion composition 
results are generally in agreement with the results of the previous workers. 
Comparisons with the observations also show some similarity in the results, 
though such a direct comparison with the observations is not really valid 
because of the several approximations made in the negative ion reaction 
scheme. 

A comparison of th^ electron c 'nsities show good agreement with the 
observations. However, there are regions where small discrepancies appear. 
Further, the systematic variation of the electron density with the solar 
zenith angle as computed is only a fraction of the observed variation. 

This is because the only term that can be included in any model calcula- 
tion which contributes to the solar zenith angle variation is the sec x 
term in the optical depth factor. However, in reality short term varia- 
tions such as temperature changes near the mesopause. X-ray flux changes 
in the 0-3 A band could cause variation in the electron densities within 
the same day. Such variations, unfortunately, do not follow a given 
pattern to be included even empirically in a computer code. 

Therefore, in conclusion, it might be mentioned that the observed 
solar zenith angle variation in the ^-region electron densities may be 


i 

i 
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considered as due to two components. Firstly, a regular variation caused 
by the sec x term in the optical depth factor, and secondly, an irregular 
variation caused by such effects as temperature variations and X-ray flux 
variations. 

A large number of input parameters go into the numerical modeling of 
the D region, as demons < ..ate d in this work. Unfortunately, many of these 
remain uncertain. In the neutral chemistry model, more reliable data on 
such quantities as the solar fluxes above 1000 A, absorption cross sections 
of the constituent gases, rate constants of fast neutral reactions, eddy 
diffusion coefficient and its height distribution are urgently required. 
Further, actual measurements of the distribution, of minor neutral species 
such as NO, 0( P) and H^O are required to verify the validity of the model 
calculation results . 

The ion chemistry model also requires data with better reliability 
for satisfactory interpretation of the measured ion composition and electron- 
density profiles. Among these are X-ray fluxes, energetic particle fluxes, 
ionization yield of 0 2 and N 2 due to these high energy sources, and rate 
constants, particularly the electron recombination coefficients. 

Also, one need not emphasize the importance of conducting coordinated 
experiments to measure the ion composition and electron- density profiles 
of the D region. Such measurements are necessary on quiet days in addi- 
tion to those on special occasions so that a clear picture of the normal 
Z?-region ionization could be obtained. It is hoped that in the near future 
these data will become available to D-region investigators. 
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APPENDIX I. NUMERICAL SOLUTION OF PARTIAL DIFFERENTIAL EQUATIONS 


The finite difference scheme adopted in the conversion of the set 
of partial differential equations (equation (3.16)] into the corre- 
sponding difference equations (equation (3.21)), and the method of 
solution of the resultant equations are given in this appendix. 

Using the spatial and time divisions given in equations (3.19) and 
(3.20), the spatial derivatives can be expressed in the form 


a y m Y m - y m . 
n ~ n+l w-1 

3s ~ 2 As 


(A.l) 


s 2 y m r m , 

n ~ n+1 


2Y m + Y 
n n-i 


9s 


(As)‘ 


(A. 2) 


where Y m is the constituent concentration at altitude s and time t . 
n n m 

The time derivative is written in the semi-implicit form 


^ j ^ Y m y ^ 

n - n ~ x n 


9t 


At 


(A. 3) 


When equatio, s (A.l) - (A. 3) are substituted in equation (3.16), it 
becomes 
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2As 


In this equation, the coefficients a and b are functions of s only, 
while the coefficients o and d are functions of both s and t. The coef- 
ficient o includes the loss coefficient and therefore its value at time 
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t m is evaluated using the concentrations Similarly, d, which 

is the production term is also evaluated using concentrations obtained 

at time t . . 

m-1 

Rearranging, equation (A. 4) becomes 



This simplifies to 


(A. 5) 
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(A. 6) 


which is the same as equation (3,21). 

At the beginning of the computations, the coefficients g and 6 are 
evaluated using the initial values of {y 

rb 

For the N height intervals, and corresponding to a given instant, 

(iV- 2) equation of the above type can be written. The resulting system 
of equations could therefore be solved using the two boundary conditions. 

In general, a boundary condition can be written in the form 


u 



VY = w 


(A. 7) 


where u 3 V and w are evaluated at the two boundaries. 
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When written in the difference form and rearranged, this becomes 


u 

/i? ♦ 

As 1 


u 




A f 1 = W 

As o o 


(A* 8) 


% + 


V \ jn 

z } N ~ s N-l 


- w 


N 


(A. 9) 


for the lower and upper boundaries, respectively. 

Following the method given by Richtmyer [1957], a recurrence formula 

is written for the variable Y m in terms of two view variables, r and s 3 

n 

so that 






s 


n 


CA. 10) 


Eliminating terms including and 7^^ from the equation (A. 6) 

using the above formula, one gets 


( — + S yr . ) Y^ + (ys , - — S - 6) = 0 
'•r 1 n-l J n n - 1 r n 

n n 


(A. 11) 


(In writing this equation, the suffix in the coefficients are left out 
for simplicity.) 

This expression is true for all values of Y ™ , so that one may write 


— + 5 + Sr . = 0 
n-1 


(A. 12) 
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This yeilds 
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n 8 + y 2 * 


(A. 14) 


n - 1 


and 


i - y s 


s = 


n - 1 


n My? 


(A. 15) 


n-1 


The starting values of r and a are determined using the lower 

n n 

boundary conditions. 


u 


jr> = 


o u — V As 
o o 


(A. 16) 


W A 3 

3 = “ T 

0 As 


(A. 17) 


Thus, starting from v and s , all the values of r and s could 
° o o n n 

be readily evaluated. 

Next, using the upper boundary condition (equation (A. 9)), one 

gets 


y _ w n Ag * S J7-1 u jy 

N ” v N As + (l-*^) u N 


(A. 18) 


Therefore, once Y^ is known, the other values of Y n could be obtained 
using the above equation, in the descending order of n. 


§. 
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